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We review our current understanding of the critical dynamics of magnets 
above and below the transition temperature with focus on the effects due to 
the dipole-dipole interaction present in all real magnets. Significant progress 



o 
in . 

ON , in our understanding of real ferromagnets in the vicinity of the critical point 



has been made in the last decade through improved experimental techniques 
and theoretical advances in taking into account realistic spin-spin interactions. 



a . 

O , We start our review with a discussion of the theoretical results for the criti- 

cs 

cal dynamics based on recent renormalization group, mode coupling and spin 
wave theories. A detailed comparison is made of the theory with experimen- 
tal results obtained by different measuring techniques, such as neutron scat- 
tering, hyperfine interaction, muon-spin-resonance, electron-spin-resonance, 
and magnetic relaxation, in various materials. Furthermore we discuss the 
effects of dipolar interaction on the critical dynamics of three-dimensional 
isotropic antiferromagnets and uniaxial ferromagnets. Special attention is 
also paid to a discussion of the consequences of dipolar anisotropics on the 
existence of magnetic order and the spin-wave spectrum in two-dimensional 
ferromagnets and antiferromagnets. We close our review with a formulation 
of critical dynamics in terms of nonlinear Langevin equations. 
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I. INTRODUCTION 



Remarkable progress has been achieved in our qualitative and quantitative understand- 
ing of the critical dynamics of magnetic materials. This is partly due to the advances and 
new developments in experimental techniques such as neutron scattering, electron spin res- 
onance and hyperfme interaction probes. Simultaneously, on the theoretical side, important 
developments were provided by dynamical scaling theory, mode coupling theory and the 
renormalization group theory. The theoretical progress has been mainly promoted by in- 
cluding effects of magnetic interactions, such as dipole-dipole and spin-orbit interaction, on 
top of the exchange interaction. These interactions, which are present in all real magnetic 
materials, were found to change the critical dynamics quite drastically. Around 1986 the 
experimental situation was quite puzzling showing both quite excellent agreement but also 
large discrepancies with the theories existing at that time. Motivated by this seemingly 
contradictory situation anewed theoretical interest on the subject of critical dynamics was 
promoted. 

In this review we shall mostly be interested in the effects of the dipole-dipole interaction 
on the critical dynamics of isotropic ferromagnetic materials in three dimensions. The static 
properties are reviewed only in as much as they are needed for the dynamics or if they 
emerge naturally together with the dynamical results. We will also discuss other dipolar 
systems, such as dipolar antiferromagnets and uniaxial ferromagnets in three dimensions, 
as well as the consequences of dipolar interactions in two-dimensional magnets. We do not 
consider pure dipolar systems such as the nuclear magnets Cu, Ag, and Au. We aim at 
giving an up-to-date account of the situation in the field, where we try to give attention 
to the theoretical as well as the experimental progress. In order to make the discussion as 
self-contained as possible we give brief discussions of the various theoretical concepts where 
needed. 

Ferromagnets have always played a special role in the field of critical phenomena. Several 
simplified models, such as the Ising |[L26|1 , Heisenberg ||106|| , and Hubbard ||120|| model have 
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been developed in order to understand the fundamental questions of the statistical behavior 
of magnetically ordered systems. The attraction of these models resides in the simplicity 
of their form and complexity of behavior they are capable to describe. Over the last two 
decades various approaches to study the above models have been intensively developed, 
which lead to a profound understanding of the statistical mechanics of magnetic systems 
on the basis of the above simplified models. However, one has to keep in mind that for 
a more realistic description of magnetic materials it is not sufficient to consider solely the 
electrostatic interaction between the electrons in the partly filled shells leading in conjunction 
with the Pauli principle to the exchange interaction. There are magnetic interactions, such as 
the mutual interaction of the electron spins (dipole-dipole interaction) and the interaction 
of the spin and orbital moments of electrons (spin-orbit interaction) which may lead to 
magneto- cryst alline anisot ropies . 

The exchange interaction is responsible for the phenomenon of magnetic ordering itself 
(at least above two dimensions). The characteristic feature of the exchange forces is their 
isotropic and short-range nature. They do not impose any definite orientation of the mag- 
netic moments with respect to the crystallographic axis. Magneto-crystalline anisotropies 
are due to relativistic or magnetic forces (spin-spin dipole, quadrupole, etc.; spin-orbital, 
orbital-orbital). Microscopic models of various magnetic interactions have shown that, in the 
majority of cases, the spin-orbital interaction is the basic one responsible for the magneto- 
crystalline anisotropy (see e.g. Ref. | 256J200| ). It relates the directions of the spin magnetic 



moments of atoms through their orbital states with the crystallographic axes. Irrespective 
of the microscopic nature of the magnetic anisotropy forces, their macroscopic manifestation 
in a crystal seems to be determined mainly by the type of symmetry of the lattice. A special 
role among the magnetic forces is played by the dipole-dipole interaction between the mag- 
netic moments of the electrons, because it is of long range and in cubic crystals leads to an 
anisotropy with respect to the wave vector but to leading order not to the crystallographic 
axes. 

A Hamiltonian which takes into account exchange interaction as well as dipole-dipole 



interaction between the magnetic moments of the electrons has first been given and studied 



by Holstein and Primakoff ||115|| . The dipolar interaction is usually two to three orders 
of magnitude weaker than the exchange interaction, and it can therefore (in many cases) 
be considered to be a small perturbation. However, because of its long-range nature, the 
effects of the dipole-dipole interaction become significant at least very close to the critical 
temperature T c and for small wave vectors q, as will become clear in this review. Note, that in 
particular it leads to the appearance of the demagnetization factors. Using renormalization 
group theory it has been shown by Fisher and Aharony |74,1,||] that the dipolar interaction 



is a relevant perturbation with respect to the Heisenberg model. They find that the short- 
range Heisenberg fixed point of the renormalization group is unstable against perturbations 
resulting from the dipolar interaction, and the asymptotic critical behavior is characterized 
by a new dipolar critical fixed point. Furthermore, the long-range nature of the dipolar 
interaction reflects itself in the Fourier transform containing a contribution of the singular 
direction dependent form q a q l3 /q 2 . An immediate consequence is the fact that longitudinal 
fluctuations are reduced in comparison to the two transverse ones and do not have a divergent 
susceptibility any more at T c . Here, by longitudinal and transverse we refer to the direction 
of the wave vector q. A consequence of this reduction of the number of effective order 
parameter components is a change in the static critical indices. Hence in a system with 
dominating exchange interaction there is a crossover from Heisenberg critical behavior to 
dipolar critical behavior. For instance the effective temperature dependent critical exponent 
of the static susceptibility 7 approaches the Heisenberg value then goes through a minimum 
until finally it ends up at the dipolar limiting value. 

In general the long-range dipolar interaction is of importance whenever fluctuations 
get large. This is the case in the vicinity of critical points and in systems of reduced 
dimensionality. In the vicinity of critical points longitudinal fluctuations are suppressed and 
rotational invariance is destroyed. This leads to modified static critical behavior and to 
drastic changes in the dynamics. In systems of reduced dimensionality which on the basis 
of short-range interactions would not have a phase transition at a finite temperature due 



to the large fluctuations destroying the order parameter, the dipolar interaction suppresses 
these fluctuations and thereby allows a finite order parameter. With the detection of high-T c 
superconductors and their fascinating magnetic properties, the study of mechanisms which 
lead to phase transitions in such quasi-two-dimensional systems (inter-plane interaction, 
anisotropy, dipolar interaction, etc.) is of prime importance. 

In this paper we review the critical dynamics of magnetic systems. As a reference and 
a simpler situation to start with we treat in chapter 2 first isotropic ferromagnets with- 
out dipolar interaction. This allows us to introduce the main theoretical concepts such 
as dynamic scaling, mode coupling and dynamic renormalization group theory in a quite 
elementary and hopefully pedagogical way. 

In chapter 3 we describe the theoretical results on the dynamics of dipolar ferromagnets 
with an emphasis on the mode coupling theory for the paramagnetic phase. A detailed 
analysis of the consequences of the dipolar interaction on the functional form of the dynamic 
scaling laws, critical exponents and the line shape and line width crossover will be given. 
For the ferromagnetic phase we give some results based on spin-wave theory and we also 
comment on some recent theoretical approaches, which go beyond the linearized spin-wave 
theory. 

In chapter 4 the theoretical results are then compared with the findings from a variety 
of experimental techniques. They include neutron scattering, electron spin resonance and 
magnetic relaxation, hyperfine techniques and muon spin resonance experiments. 

Whereas the main part of the review concentrates on dipolar effects in isotropic ferro- 
magnets, chapter 5 concerns dipolar effects in other magnetic systems. These include three- 
dimensional isotropic antiferromagnets, bulk uniaxial ferromagnets and two-dimensional 
systems. In the latter case the dipolar interaction leads to long-range order, which would 
not be possible for the isotropic Heisenberg model in two dimensions. 

Finally, in chapter 6, we discuss alternative derivations of the mode coupling theory 
via the generalized Langevin equations of Zwanzig and Mori. These methods allow for 
a systematic derivation of the mode coupling theory in the framework of a diagrammatic 
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analysis. Also, some systematic improvements are possible. We conclude with a summary 
and an outlook in chapter 7. Some technical details and important conceptual background 
material is collected in the appendices. 
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II. ISOTROPIC FERROMAGNETS 



Our main concern in this review is the immediate vicinity of the critical point and the 
influence of the dipolar interaction on the critical dynamics. As a reference and a simpler 
situation to start with we treat first isotropic ferromagnets without dipolar interaction. This 
allows us to introduce the basic concepts and results of the different theoretical approaches 
and compare them with the experimental situation. 



A. Dynamical Scaling and Hydrodynamics 

In systems where a continuous symmetry is broken, hydrodynamics together with dy- 
namical scaling allows one to obtain definite conclusions about the dynamic critical behavior. 

To start with, we remind the reader of the structure of the hydrodynamic modes in 
isotropic ferromagnets. In an isotropic ferromagnet magnetization is conserved giving rise 
to three hydrodynamic equations for the magnetization vector M(x). In the paramagnetic 
phase, i.e. for temperatures above the Curie temperature T c and in zero magnetic field H, 
(T > T c , H = 0) the magnetization obeys a diffusion equation, 

^ = DV 2 M, (2.1) 

with a diffusion constant D. Hence the long wave length excitations are diffusive 

u(q) = -iDq 2 . (2.2) 

In the ferromagnetic phase the magnetization M is finite and the spin fluctuations which 
are perpendicular to the mean magnetization obey spin-wave equations of motion. The 
spin-wave frequency is according to hydrodynamics given by 

^q) = ^-^q A A■ (2.3) 

Here x T {q) is the transverse susceptibility, M the magnetization, and A a damping constant. 
The real part of the frequency is related to static critical quantities. For the imaginary part, 
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the damping, hydrodynamics predicts a decay rate proportional to the fourth power of the 
wave number []. 

Dynamical scaling states that the critical frequency is of the homogeneous form 

= <z*n(?£), (2-4) 

where z is the dynamic critical exponent and £ oc |T — T c \~ u the correlation length with the 



static critical exponent v. Using the scaling behavior of the static quantities ||175|| in the 
hydrodynamic region (M ~ and x T {q) 

~ q 2 £ v ) and the scaling relations between 
the static exponents (7 = v(2 — 77), 2(3 = (d — 2 + T])u) one finds for the spin-wave frequency 
in the hydrodynamic region 

Re M<?,0} = M/ X T (q) oc q ^~ d ^ = ^-^^-d+r,)^ (2 5) 

Here d is the spatial dimension of the system. One thus finds as a result of the dynamic 
scaling law, Eq. ( |2.4|) and the hydrodynamic behavior of the static quantities an exact 
relation for the dynamic exponent 

z= d + 2 -". (2.6) 



Using again Eqs. Q2.5|) and (|2.6|) the damping coefficient A of the spin-waves and the spin 



diffusion coefficient D in the paramagnetic phase are 

A^^-d+v)/^ and _D ~ j:( 2 - d +v)/2 _ (2.7) 

The spin diffusion coefficient D goes to zero at the critical temperature, a phenomenon 
known as critical slowing down. Thus hydrodynamics and dynamical scaling allow one to 
determine the main critical dependencies of the transport coefficients. The picture emerging 
from Eqs. (|2.2|) -(277) is summarized in Fig. 2.1. 



indeed the microscopic theory of the Heisenberg model predicts a spin-wave decay rate of the 
form q A (c 2 (lnq) 2 + c 1 lnq + c ) &33|,|102|i257|1 . 
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The critical frequency is a function of the wave number and the inverse correlation length. 
The hydrodynamic region is given by q <C whereas the non hydrodynamic critical region 
by q S> Of course the critical region is limited to q <C a" 1 and £ ^> a, where a is a 
microscopic length scale (e.g. the lattice spacing). 

We close this subsection by giving the hydrodynamic equations of the low temperature 



phase [|244| . |112 



dt « X T (q) 



M*= -^M* - Aq 4 M* , 



—M y = — M x - Aq^M y , 

dt q x T (g) q q 



d 
dt 



-M*= -T(g)M^ 



the excitations of which were the basis of our discussion. Here we assumed that the mag- 
netization is oriented along the ^-direction. In addition to the transverse equation the 
magnetization component along the order parameter obeys a diffusion equation in three 
dimensions T(q) = Duq 2 ■ 

B. Mode Coupling Theory 

One of the most successful theoretical approaches in critical dynamics are mode coupling 
theories. The first such theory was proposed by Fixman [ff5| |, and then put on a more 



rigorous basis by Kadanoff and Swift ||130|| , and especially Kawasaki ||137| , |138| , |140|| . In this 



section we exemplify the mode coupling theory for the critical dynamics of an isotropic 
ferro magnet, where the spins are coupled only by short range isotropic exchange interaction. 
The Hamiltonian for such a spin system is given by 

H = I J(q)S q -S_ q , (2.8) 

Here we have introduced the notation J* q = v a J ^yr, where v a is the volume of the primitive 
cell of the Bravais lattice. Upon introducing the cube edge length a of the corresponding 



cubic unit cell, the dimensionless quantity b = ya 3 /v a characterizes the lattice structure. 
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The Fourier transforms of the Cartesian components S l (x) of the spin operator are defined 
by 



4 = / rf 3 xe iq ' x ^(x) 



(2.9) 



and J(q) = —Jo + Jq 2 a 2 characterizes the exchange interaction. We are retaining only terms 
up to second order in the wave vector q and have supposed that the exchange interaction 
extends up to the second nearest neighbors. For bcc and fee lattices we have J = J\ + J2, 
where J\ and J2 are, respectively, the values of the exchange parameters between the nearest 
and next-nearest neighbors. For a simple cubic (sc) crystal the relation is J = J\ + 4J2. 
The parameter J does not enter the equations of motion. It is convenient to introduce the 
ladder operators 



(2.10) 



and 5*q instead of the Cartesian components of the spin operator. Then, using the com- 
mutation relation for spin operators one finds for the Hamiltonian ( |2.8| ) the following set of 
equations of motion (we take % = 1) 



d_ 

dt 

d 
dt 



S^-if [ J(k) - J(q - k)] Sj_ k 5£ , 
S± = ±2i I [ J(k) - J(q - k)] S±_ k Sl . 



(2.11) 



(2.12) 



The equations of motion ( |2.11| ) and ( p. 12 ) exhibit explicitly the vanishing of dS^/dt at q = 0, 



i.e., the order parameter itself is a constant of motion. Starting from these microscopic 
equations of motion there are a variety of different schemes for deriving the mode coupling 



19,262 


232 


233 


137 


118 


119 



The quantity of interest is the Kubo relaxation matrix for a set of dynamical variables 
{X a (x)}, which is defined by 



JO 

$ a/3 (q,t) = Him / rfre- £r ([X Q (q,r),X /3 (q,0) t ]) , 

e^O J 



(2.13) 
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where (...) denotes the thermal average. The dynamical variables are normalized such that 

(X a ,X ) :=$ a/3 (q,t = 0) = 5°P, (2.14) 



i.e., we use normalized spin variables X®(t) = 5 , q(t)/yx a (q), where x a (<l) are the static 
susceptibilities. 

One of the most concise ways of deducing mode coupling equations utilizes the projection 
operator technique, originally introduced by Mori ||197|| and Zwanzig ||276|| . The main idea 



is that one can separate the set of dynamical variables into two classes, one slowly and one 
fast varying class. With the aid of this projection operator method the fast variables are 
eliminated and one can derive generalized Langevin equations for the dynamic variables or 
equivalently for the correlation functions | |199|jl98| , |139|| (see also chapter [VT|, and appendix |TJ) . 
The corresponding equations for the Kubo relaxation function is 

^i = iw «rr( q)$ 7/J( qjt )_jf rf r r^(q,t-r)$^(q,r), (2.15) 

and for its half-sided Fourier transform 

$ a/3 (q,cj) = J dt^(q,t)e iujt , (2.16) 
o 

one obtains 

$(q, u) =i[ul + w(q) + iT(q, cu)}' 1 . (2.17) 
The frequency matrix u al3 (q) is given by 

' Q/3 (q) = (X«,X?) = -i([X«,X^) , (2.18) 



ILO 



where we have used the Kubo identity (A, B) = —i([A, B^]) ||157| . The non-linear aspects of 
the spin dynamics are contained in the matrix of the transport coefficients (memory matrix) 
r. As a result of the projection operator technique these can be written in terms of the 
Kubo relaxation matrix [ 138| , 140 | 



r#{q,t) = (6XZ(t),6XZ(t)) (2.19) 
13 



of the non conserved parts of the currents^] 

5X% = X%- iuj a/3 (ci)XP . (2.20) 

The simplest approximation which can be made at this stage is to consider only two mode 
decay processes, which in technical terms amounts to a factorization of the Kubo formulas 
Q2.19| ) after insertion of the equations of motion. Frequently one makes additionally an 
approximation for the line shape (i.e. frequency dependence) of the relaxation matrix, e.g. 
a Lorentzian approximation. In principal, however, one can solve directly the set of self 
consistent equations for the shape functions resulting from the decoupling approximation 



only. For the isotropic ferromagnet this was first achieved by Wegner ||262|| and Hubbard [|118 
in the paramagnetic phase. For most practical purposes an excellent approximation for the 
line width can be obtained from the mode coupling equations simplified by the Lorentzian 
approximation. 



1. Paramagnetic Phase 

In the paramagnetic phase the order parameter is zero < S q >| q= o= implying that the 
frequency matrix uj al3 vanishes. Upon using the above decoupling procedure one obtains the 
following set of coupled integro-differential equations for the Kubo relaxation function 

a$(g,t) 

dt 

and the transport coefficients 



/ drr(q,t-r)$(q,r), (2.21) 
Jo 



r(q, t) = 4k B T I v(k, q) X(q k)x(k) $(k, t)$(q - k, t) , (2.22) 
A y q 



x(q) 



2 Here we have chosen a linear projection operator VX = (X, X a )X a . Then the random forces 
can be written in terms of the projection operator as <5X q = exp[i(l — V)Ct](l — V)Xq, where 
C is the Liouville operator. Neglecting the projection operator in the time development one gets 
Eq. fl2.20| ); see also Appendix B. 

14 



with the vertex function (H = 1) 



u(k, q) = (J(k) - J(q - k)) z = [ 2Jq 2 a 2 



q- k 1 

~ 2 



(2.23) 



Essentially the same equations have been derived by numerous authors [19|, 262, 232 , 233 



||137| , |140|| , ||118| , |119|| using different approaches. The temperature dependence enters the 
equations only implicitly via the correlation length 

'T-T r 



[2.2A) 



The equations for the Kubo relaxation function must in general be solved numerically. 
However, in the critical region one can deduce certain important properties of the solution 
analytically. Upon inserting the static scaling law [ [269|| 

1 



x(q,0 



2Ja 2 



(2.25) 



with the scaling variable x = l/q£ one can show by inspection that the solution of Eqs. ((2.21[) 



and ( |2.22 ) fulfills dynamic scaling 



u;) = (Aq 



g\-l 



(2.26) 



T(q,£,uj) = Aq z 1 (x, v) 



with the scaling variable v = u/Aq z and the non universal constant 



A 



x 5 / 2 l2Jk B T c 



(2.27) 



(2.28) 



b V 4vr 4 ' 

where b is a dimensionless parameter which depends on the crystal structure (see Table |) . 



Note that in Eqs. ( 2.25| )-( 2.27|) we have explicitly incorporated the correlation length £ 
into the list of arguments of the correlation functions in order to indicate the reduction of 
arguments accomplished by the dynamic scaling form. In order to simplify notation this 
temperature dependence is in most of the remaining text not written out explicitly. 

The above mode coupling equations give a dynamic critical exponent z = (5 + 77) /2 
instead of the correct expression z = (5 — rj)/2 |99| , |174| , [T27| , |r^] . This inconsistency of the 
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conventional derivation of the mode coupling equations is fortunately not a very serious 



problem since r\ is very small in the case of 3D ferromagnets, r, ~ 0.05 [ |189|| . In order to be 
consistent one has to take for the scaling functions an Ornstein-Zernike form 

= > ( 2 - 29 ) 

1 + X 1 



neglecting the exponent q. In chapter we will give a derivation of modified mode coupling 
equations based on a path integral formulation of the stochastic equations of motion. There 
we will show how the above inconsistency can be resolved by taking into account certain kinds 
of vertex corrections which are neglected in the conventional derivation of mode coupling 
equations. 

The scaling relations ( |2.26| ) and ( |2.27| ) for the Fourier transformed quantities imply for 
their time dependent counterparts 

*{q,Z,t) = <Kx,T), (2.30) 

T(q^ J t) = (Aq z ) 2 1 (x,r), (2.31) 

with the scaled time variable 

r = Aq z t. (2.32) 

The mode coupling equations for the corresponding scaling functions are 

d(p(x, t) 



o 



dr 
and 

-' roc x(x/p)x{ x/p-) 

-i 'Jo X{x 

The scaled vertex function reads 



dr'-i{x,T -t')<P{x,t'), (2.33) 



j(x,t) = 2vr 2 / dr, / dpp 2 _v(p, r,) XK ^ "JX' ^ </>(x/p, rp z )<j>{x/p^ rp z _) . (2.34) 



v(p,r ] ) = 2(p V -l/2) 2 , (2.35) 



where we have defined p = k/q, p_ = |k — q\/q, and r, = cos(k, q) 
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Before turning to the numerical solution of the mode-coupling equations, let us quote 
some results which can be obtained analytically. For temperatures not to close to T c one can 
infer from Eqs. ( |2.22| ) and ( |2.23| ) that in the limit q — > (hydrodynamic limit) the vertex 
factor v(h, q) and hence the memory kernel T(q,t) becomes small. Hence one could argue 
that the relaxation function <3>(g, t) varies very slowly and the solution of Eq. fl2.21| ) becomes 
an exponential ||118 



$(q,t) = exp{-Dq 2 t} , (2.36) 

where the diffusion constant is given by 

1 r°° 1 
D = lim— / T(q,t)dt = lim — T(q,u = 0) . (2.37) 
g^o q z J q^o q l 

A scaling analysis of the right hand side gives for the temperature dependence of the diffusion 

coefficient D ~ ^~ l l 2 in agreement with the scaling result by Halperin and Hohenberg |99] . 

The above argument leading to the spin diffusion behavior has been questioned by 

Manson | 180 |. Starting from a relaxation function which is of spin diffusion type (see 



Eq. fl2.36| )) Manson ||180|| showed that the memory kernel becomes of the form 

T(q,t) ocr 5 / 2 e- Dg2t/2 (2.38) 

for asymptotic times and small q. Therefrom Manson ||180f| concludes that the spin diffusion 
type of behavior can not be the correct form of the relaxation function at asymptotic times, 
which would raise some questions on the validity of the mode coupling theory (since it 
invalidates the results obtained from a hydrodynamic theory based on the conservation of 
the magnetization). The above argument leading to an exponentially decaying relaxation 
function becomes invalid also close to T c because the static susceptibilities in the expression 
for the memory matrix diverge and therefore the relaxation function no longer varies slowly. 
A shape crossover from a Lorentzian to a different critical shape takes place by approaching 
the critical temperature ||118|| . Nevertheless, a Lorentzian approximation for the line shape 
still gives reasonable an approximation for the line width, since the latter is not so sensitive 
to the precise form of the line shape. 
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In the Lorentzian approximation the mode coupling equations reduce to a single integral 
equation for the line width ri or (g) = Ag z 7 lor (x) 

7l ° r(X) = Jix)L dV L dpP - V ^ } / x ■ 5/2 • ( 2 - 39 ) 

Therefrom one can deduce the asymptotic behavior of the typical line width analytically 

{1 for x << 1 , 
(2.40) 
y/x for x >> 1 , 

implying T(q) ~ q 5 ^ 2 right at T = T c and T(q) ~ g 2 £ -1 / 2 in the hydrodynamic limit g£ << 1, 
i.e the temperature dependence of the diffusion constant is given by D ~ ^~ 1 / 2 ||137|| as we 
have already deduced from scaling arguments. The full scaling function resulting from 
Eq. ( f2.39|) is shown in Fig. 2.2. It is usually called Resibois-Piette scaling function ||23 1 



since Resibois and Piette did the first numerical solution of the mode coupling equations in 
Lorentzian approximation. 

In order to find the complete behavior of the relaxation function $(g, t) one has to solve 
Eqs. ( [2.33D and ( |2.34| ) numerically. This was done by Wegner ||262|| at T c and extended 



to temperatures above T c by Hubbard ||118| , |119|| . The results are shown in Fig. 2.3. It is 



found that for small wave vectors (not too close to the zone boundary) there is a shape 
crossover from a Lorentzian (see Eq. (|2.36f) ) to a more Gaussian like shape by approaching 
the critical temperature. The critical shape at T c is essentially the same as the one obtained 
from RG-theory plj^ll (see also section [II C 1|) 



Recently, this shape crossover has been reexamined by Aberger and Folk |l] and Frey et 
al. |85| in detail with emphasis on constant energy scans. Their results, shown in Figs. 2. 4. a 



and 2.4.b for the scaling function of the spin relaxation function versus time and frequency, 
respectively, confirm the shape crossover from a Lorentzian to a critical shape first found 
by Hubbard | |118|J119|| . In addition, strongly over-damped oscillations in the time-dependent 



spin relaxation function at T c are found. These oscillations, however, do not lead to an 
observable structure in the Fourier transform apart from a flatter decrease at small frequen- 
cies. With increasing temperature these oscillations practically disappear. Presently it is 



not clear, whether these oscillations are an artifact of the mode coupling approximation 
and go away when higher order terms are included. A RG analysis does not show these 
oscillations P,|2llp3|. 



Hubbard ||118|| also discusses the shape function for cubic ferromagnets with nearest- 
neighbor interaction for wave vectors close to the Brillouin zone boundary: J(q) ~ 
J[cos(q x a) + cos(q y a) + cos(q z a)}. He finds that there is a tendency of the shapes to become 
more squarer than a Lorentzian and as the wave vectors come close to the zone boundary one 
observes the formation of small shoulders. Recently Cuccoli et al. []54| have studied the shape 
of the correlation function at the zone boundary for EuO and EuS with a face-centered- 
cubic lattice taking into account nearest- and next-nearest-neighbour exchange interaction. 
The numerical solution of the mode coupling equations give, as in the simple cubic case 



with nearest-neighbor interaction considered by Hubbard ||118|| , inelastic shoulders at the 
zone boundary, but less intense than seen in the experiment [EflpI • 



2. Ferromagnetic Phase 



In this section we review the mode coupling equations for isotropic ferromagnets below 
the Curie point [p45|j8l|] . Assuming that the spontaneous magnetization points along the 
z-axis the frequency matrix is given by 

/0 \ 



u^{q) = u(q) 







-1 
\o +1J 

where a, (3 = z, +, — . The frequency of the transverse modes is 

M 



(2.41) 



x t (q) ' 



(2.42) 



where M = {S% =0 ) denotes the magnetization and x T (<?) the static transverse susceptibility. 
Due to the rotational symmetry of the Hamiltonian the Kubo relaxation matrix Q al3 (q, u) 
is diagonal 
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g"(g,fa;)= \ ' (2-43) 

uj + zl 2Z (g, cjJ 

$±±( g , w ) = ^zx7 \ • (2.44) 

The mode coupling approximation for the transport coefficients 

r zz (g,t) = (^(t),s*(o))/x%) = r( g ,t) , (2.45) 



V^{q,t) = (S±(t) ±iu(q)S±(t),S±(0) ±tuj(q)S±(0))/2x T (q) = A ± (g,t) , (2.46) 

where x L (l) is the longitudinal susceptibility, results in the following set of integral equa- 
tions ||245|| 

,T{„ U\^T 



T(c l ,u J ) = k B T I I t;(k,q) X (q (k) $+ + (q - k, a; - i/)$"(k, i/) , (2.47) 

.a, a x (q) 



A ± (q, a;) = 2A; B T / / u(k, q) ^ k) f L(k) ^(q -k,u- v)<S> zz (k, v) . (2.48) 

Ju Jk x (q) 



Here we have used the notation /„ = /§£. Eqs. (^43|) , ([24l together with (|2T47|) , (^48 



constitute a complete set of self consistent integral equations for the Kubo relaxation func- 
tions $ Qa (g, u), which in principal could be solved numerically. For M = and x L = X T 
Eqs. ( |2.47| ), (|2.48| ) reduce to the mode coupling equations for the paramagnetic phase, 



Eqs. (ggg ). 

The above mode coupling equations have been analyzed in the Lorentzian approximation 
for the relaxation functions 

«) = ' ^(g, w) = - f , (2.49) 

with 

r(g) = r(g,w = 0), A(g) = A + (g) = A~(g)* = A+(g,u;(g)) . (2.50) 

The frequency integrals can now be carried out readily and one finds the following set of 
coupled integral equations for the line widths 
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r(q) = ^/ «(k,q) 



X T (q-k) X r (k) 



w(k) - cj(q - k) + zA(q - k) + £A*(k) 



A(q) = ^ / r(k.q) 



x T (qj ^ 



X T (q-k) X i (k) 



(2.51) 



^2.52) 



w(q) - cu(q - k) + zA(q - k) + ir(k) ' 

As is easily seen T(g) is real, but A(q) in general is complex. The imaginary part of the 
transverse damping function A(q) leads to a shift of the frequency of the transverse spin- 



waves, which however is a negligible correction in comparison to the frequency matrix ( 2.41|) 
as will be seen later. 

In the hydrodynamic regime the Eqs. Q2.51|) , Q2.52 ) can be solved analytically with the 
result 



r(g)oc 



q 



x L (q) 



and A(g) oc q 4 



(2.53) 



where Co and C\ are constants []. With the well known scaling properties of the static 
susceptibilities neglecting the Fisher exponent r\ 

1 



X L ' T (q) 



X L ' T (x) 



2Jq 2 a 2 



(2.54) 



Eq. ( pTITD gives 



uj(q) = Aq z Lo(x) 



[2.55) 



where the dynamical critical exponent z equals 5/2 as in the paramagnetic phase. 

The scaling function for the bare frequency of the transverse modes following from 
Eq. fl2.42|) can be written as 

fy/x for T <T C , 



OJlX) 







for T>T r 



(2.56) 



3 The low-temperature spin-wave theory of Dyson 65] gives in addition a term 
<? 4 (ln(l/^)) P^Wp57| . This will probably come out from a mode coupling theory where 



decays of the transverse mode into three transverse modes are included. 
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Analyzing the scaling properties of the mode coupling equations and combining this with 



the static and dynamic scaling law it was shown |[24(]| , [241|| that the amplitude / for the 
scaling function of the spin-wave frequency is a universal quantity and determined by other 
universal amplitude ratios 



/ 



c\ l/2 ( R 



1/2 / t \ d-2 / j, \ z-2 
T 



(2.57) 



Here c is an arbitrary normalization constant for the scaling functions. If one chooses the 
value of the scaling functions at criticality to be 7(0) = 5.1326, c becomes c = 8tt 4 ||241|| . 



The quantities R c and Rt are universal amplitude ratios as defined in the review article by 



Privman et al. [226], £j is a transverse correlation length below T c [226], and £ + , £_ are 
longitudinal correlation lengths above and below T c , respectively. 

The amplitude / for the spin-wave frequency can be determined from RPA-arguments 
(see e.g. Ref. |[245|| ), which gives / = 7r 3 / 2 . This value has been used in consecutive applica- 



tions of mode coupling theory on magnets [51] below T c . Upon using the known values for 
the static amplitude ratios | |226| | it is found [ |24U| J 



/ = 9.5 ± 1.8. 



(2.58) 



This amplitude can also be determined from the available experimental data in 
Refs. [p9| , p6| , p0| , |34| , |220|| , and references cited therein. The results are summarized in Ta- 
ble ||240|J241|| , where depending on which experiment [p9] , p6|p0] , p4] , p20|| one analyzes one 
gets slightly different values for /. 

Hence Eqs. ( [2.51|) , ( ]2.52|) can be solved by using an dynamic scaling Ansatz 



T(q)=Aq z 1 (x), A(q)=Aq z \(x), 



[2.59) 



where the dynamical scaling functions j(x) and X(x) obey the following set of coupled 
integral equations 

X T (x/p-)x T (x/p) 



7(x) = 2ir i drj dpp_ v(p,rj) 



-1 



X L (x) 
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—p z £o(x/p) + p z _u)(x/p-)+ip z X{x/p)+ip z _X*{x/p^)'' ^ 

X(x) = 2ir z i / dr] dpp_v(p,i]) ^-r 

J-i j % (:rj 

x 1 (2 61) 

uo(x) — p z uo{x/p) + ip z X{x/p) + ip z _^(x/ p-) 

Here we have used the same notation as in section [II B 1| . 

In order to solve those mode coupling equations one has to know the static susceptibilities. 
In the ferromagnetic phase the global continuous rotation symmetry is spontaneously broken. 
Although one of the equivalent directions of the order parameter is selected, no free energy is 
required for an infinitesimal quasi-static rotation of the magnetization vector, which in turn 
leads to a diverging transverse correlation length. This physical effect is mathematically 



expressed by the Goldstone theorem [0, stating that there is exactly one massless mode 
for each generator of the broken-symmetry group. In the context of a ferromagnet below T c 
this implies that the transverse susceptibility is given by 

™ = W*- (2 ' 62) 

The longitudinal correlation functions entering these integral equations has been computed 
by Mazenko [184] to first order in e = 4 — d using Wilson's matching technique 



' - 2Jq 2 a 2 



x L {q) 



9e of-, n r^7, /Vl + 4a; 2 - T 



-x 



1 + Vl + 4x 2 In 



n + 8 \ V 2x 



+x 2 ( n + 8 + (5-*)e 



, . (2.63) 
9 + (n-l)x e /J' v ; 

where n is the number of spin components. The last term in x L (q) results from the presence 
of Goldstone modes below T c , and it implies that also the longitudinal susceptibility diverges 
in the limit q — ► for any temperature below T c . 

The resulting numerical solution of the mode coupling equations ( |2.60| ), ( |2.61| ) have been 
achieved in Ref. [pT| , |240| , p41[| . The results are shown in Fig. 2.5 for three different values of 
the frequency amplitude /. 

One recognizes that the scaling function Im\(x) for the frequency shift of the transverse 
modes is very small compared to £o(x). In the critical region ImX{x) starts at the critical 

23 



point with infinite slope and is negative in the hydrodynamical region. The scaling functions 
for the longitudinal and transverse line widths split off linearly at the critical temperature 
and differ by orders of magnitude in the hydrodynamic region. This linear split-off of the 
longitudinal and transverse widths and the infinite slope of the frequency shift at the critical 
temperature below T c is an immediate consequence of the presence of Goldstone modes below 
T c . This feature can be derived analytically from Eqs. Q2.51|) , Q2.52j ). The sign of the slope of 
the longitudinal line width depends on the magnitude of the amplitude / for the frequency 
of the spin-waves. For values of / close to the RPA-value the slope is positive. If this value 
is increased towards the universal value determined in Ref. ||240|J241|| the slope becomes 



negative and one gets a minimum in the longitudinal scaling function 7. This minimum 
has been observed in a recent experiment by Boni et al. |34| (see below). Above T c the 
scaling function for the line width starts quadratically in agreement with a renormalization 
group calculation by Iro ||123|| , but in contrast to the numerically found infinite slope of 



Hubbard ||118|| . It disagrees also with a computation of Bhattacharjee and Ferrell ||22|| , who 
predict, using Ward identities, a linear dependence on 

The numerical data can be fitted in the limits x ^> 1 (hydrodynamical region) and 1 < 1 
(critical region) by simple approximants as summarized in Table |TJ (note that all functions 
are given in units of the value at criticality 7(0) = ReX(0) = 5.1326) 

In unpolarized neutron scattering experiments on Fe [f|8|], Ni ||194|| , and EuO |[210fl no 
quasi-elastic peak from spin diffusion, as predicted by the mode coupling theory [|3T 



was 



discernible. Only the side-peaks originating from the transverse spin waves were observed. 



This is plausible in the light of the mode coupling results in Ref. [81|. In the hydrodynamical 



region (x = l/q£, 3> 1) the width of the longitudinal peak is much wider than the separation 



of the transverse peaks ||245|| . Moreover, its intensity is smaller than that of the transverse 
magnons, which altogether implies that it may be very difficult to distinguish the longitudinal 
peak from the background. In the critical region the line widths are of the same order of 
magnitude. In this limit however the frequency of the transverse modes tends to zero. Using 
unpolarized neutrons one can only observe a superposition of the peaks. Lacking a theory 
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for the line width in the critical region below T c it was impossible up to recently to resolve 
the longitudinal and transverse peaks. 



The first observation of the longitudinal peak was reported by Mitchell et al. ||195|| using 
polarized neutrons. This study shows in agreement with the theory that the width of the 
quasi-elastic longitudinal peak becomes comparable with the spin-wave energy explaining 
why this peak was not observed by by neutron scattering experiments with unpolarized 
neutrons. However, there are not enough data as yet to permit a quantitative comparison 
with the theoretical predictions. Furthermore, the material is disordered (palladium with 
10°/ o iron) which makes it not an ideal system ||34)| . Very recently, Boni et al. |34| have 
investigated the spin dynamics of a Ni single crystal by means of polarized neutron scatter- 
ing. They observe that the longitudinal fluctuations are quasi-elastic in agreement with our 
theoretical predictions ||245|j81[| and RG calculations ||1 74|| . In Fig. 2.6 we show a quantitative 
comparison of the longitudinal line width, obtained from solving the mode coupling equa- 
tions in Lorentzian approximation p74) , p40|J241|| , with the experiment []34]]. The light dashed 
line represents the result of the mode coupling equations with an amplitude / = 7r 3 / 2 , taken 
from RPA-arguments. The solid line and the dot-dashed line represent the results of solving 
the mode coupling equations with an amplitude of / = 5.1326 x 1.49 and / = 5.1326 x 1.60, 
respectively. The agreement between theory and experiment is quite well for an appropriate 
choice of the universal amplitude. One should especially note, that the scaling function of 
the longitudinal line width shows a minimum in accord with the experimental data. 

Finally we note that the above analysis does not take into account effects from the 
dipole-dipole interaction. Those effects have up to now not been studied quantitatively in 
the ferromagnetic phase, but, one may expect similar effects as above T c , which we are going 
to describe in the next chapter. 



C. Renormalization Group Theory 
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1. Paramagnetic Phase 



Renormalization group calculations of the critical dynamics of ferromagnets start from 
a stochastic equation of motion for the spin density S(x, t) 



(2.64) 



where £(x, t) is a random force with a Gaussian probability distribution with zero mean and 
variance 



(C(x, t)C j (x', t')> = 2Tk B T5 {3) (x - x')5(t - t')5 ij . 

The effective Landau-Ginzburg- Wilson free energy functional is given by 

1 



(2.65) 



H = J d d x - (rS 2 + (VS) 



u 
4!' 



j2\2 



^2.66) 



These equations can be derived ||198| , |139| , |199|| using a Mori-Zwanzig projection opera- 
tor formalism ||197| , |276|| (see also chapter [VI]). An exhaustive discussion of these semi- 



phenomenological equation of motion can be found in the article by Ma and Mazenko [174 



The first term in Eq. (|2.64|) describes the Larmor precession of the spins in the local mag- 
netic field, 5H./5S, and the second term characterizes the damping of the conserved order 
parameter. The precession term of the spins in the local magnetic field plays a major role 



in the dynamics. From the RG analysis in Refs. [|174|| L7| one can infer that its effect can 
be ignored above the upper critical dimension d c = 6, and can be treated by perturbation 
theory in e = 6 — gQ. 



The RG-theory proves |p.74jJ17|| the dynamical scaling hypothesis |70| , |99|| and shows that 
the spin correlation function fulfills the dynamical scaling form near the fixed point 



4 The upper critical dimension for the corresponding static problem (see Hamiltonian in Eq. ( 2.65 )) 
is d c = 4. Since the precessional term is of second order in the equations of motion corresponding 
to third order in the Lagrangian the upper critical dimension is shifted to d c = 6. 
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(x, v) , 



(2.67) 



with the scaling variables x = l/q£ and 



V 



where the characteristic frequency has itself the scaling form 

w c (g,£) = Aq z V(x). 



[2.68) 



(2.69) 



The dynamic critical exponent z is known exactly from RG theory [|174|J1271|T7|j 

d + 2-7] 



^2.70) 



in accord with the general dynamic scaling considerations of section [11 A| . Here r\ is the 
Fisher exponent from the static scaling law 



x{q,0 = q 2+v x(x) . 



The Fourier transform of the spin correlation function can be written in the form 



(2.71) 



C(q,Z,u) = x(q,0^rRe- 



[2.72) 



Xq z —i& + \x(x)IL(x,u))] 1 

where H(x, Cj) is the self-energy of the dynamic susceptibility and we have defined Cj = uj/\q z . 
The asymptotic behavior of the self-energy is known exactly from a RG-analysis [|17j .|6l|l 

One-loop RG-calculations give to order 0(e) 



H(x, Cj) 



[2.73) 



with pTjP^I 



F{0,iu) 



U(x, Cj) = 1- e[F(x, ub) + ]- In 2] + 0(e 2 ) 



| ln(— icD) - |ln2 + 0(c£;- 1 / 2 ) for Cj > 1 
(tt + |)/8 + 3(2 + 21n2-vr)f for Cj < 1 



(2.74) 



(2.75) 
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and 111 741,11 23 



-ilnx- (In 2 + 1)12 - ^oT 2 for x > 1 , 
F(x,0) = <; 2 4 128 ( 2.76) 

-(tt + f)/8- (3tt- f -3 In 4)%- for x < 1 . 

In an e-expansion with respect to the upper critical dimension d c = 6 the 4 -term is irrel- 
evant for 4 < d, hence the static critical behavior is classical. This is no more the case for 
d < 4, a fact which has to be kept in mind if one extends the results of the RG e-expansion to 



e = 3. The explicit form of F(x, iuj) to order 0(e) can be found in Ref. ||123|| . The logarithms 
in the limits (x — > 0, v — > oo) and (x — > oo, v — > 0) are the 0(e)-contributions of the power 
law behavior in Eqs. ( |2.73| ) which are known exactly for any dimension d. Exponentiating 
these logarithms in such a way that the exactly known asymptotic behavior Eqs. fl2.73| ) is 
matched one obtains pl| , |124j , p!25fl the two-parameter interpolation formula 



H(x, u) 



(l + br 



2-e/4 

— aiCo 



e/(8-e) 



(2.77) 



with a = 0.46 (a = |(6 + 61n2 — 3n)) and b = 3.16 for e = 3. This reasoning for expo- 
nentiating the leading logarithms of a first order e-expansion to match the exactly known 
asymptotic results is due to Bhattacharjee and Ferrell One should note, however, that 
this exponentiation procedure is not unique. To test the validity of the exponentiated ex- 
pression one would have to calculate contributions to order 0(e 2 ) which would be quite 
cumbersome a task. 
The shape function 

<p(x,u)=2Re — ttt-^ttt — rr~i — t (2.78) 

shows the crossover from the critical shape at x = to a Lorentzian at x = oo in agreement 



with the mode coupling results by Hubbard ||1 18|| and a more recent reanalysis by Aberger 



and Folk 0. From Eqs. ( E2~72"D one can also determine the half width at half maximum u r 
defined by 

C(q,Z,u c ) = ±C(q,£,0). (2.79) 
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In Fig. 2.7 the scaling function u> c (x) resulting from Eqs. ( |2.72| ) and ( |2.79| ) is compared with 



mode coupling results ||231| ,|Tj|. Whereas the MC-result in Lorentzian approximation ||231 



shows large deviations from the RG-result in the hydrodynamic limit, the complete MC- 
result abandoning the Lorentzian approximation |I| follows closely the RG scaling function. 
One should note, however, that if all scaling functions are rescaled in such a way that they 
coincide in the hydrodynamic limit the differences between the scaling functions appear 
much less pronounced. 

Let us now compare the theoretical predictions with the experiment. In early neutron 
scattering experiments almost all data have been fitted by a Lorentzian shape function. 
Recently, however, with the advances in neutron scattering techniques leading to higher 
intensities and better resolution, deviations of the measured spectra from a Lorentzian have 
been observed. By comparing the RG-result at T c [|6T] , pTj| with constant energy scans on Fe 



|264| , |171| , |1 72|| it has been shown |77j that theory was in accord with the data in the observed 
experimental wave vector and frequency window. As has been demonstrated in Ref. |^] the 
peak positions as well as the peak profile in constant-energy scans of EuO could be explained 
on the basis of RG-theory, taking into account short range exchange interaction only, for 
wave vectors in the range 0.15A -1 < q < 0.3A -1 and energies 0.2meV < uj < OAmeV. 

If a model based on the exchange coupling between neighboring spins is the correct 
description of the critical behavior of real ferromagnets, one would have expected that the 
results from RG and MC theories would become even closer to the experimental data as 
one comes close to the critical temperature or/and for very large wavelength. It came as a 



completely unexpected surprise, when Mezei [|190|1 found in spin echo experiments on EuO 



that the observed shape at even smaller wave vectors q = 0.024A -1 clearly resembled a 
Lorentzian shape in disagreement with the predictions of RG theory for the dynamics of an 
isotropic Heisenberg ferromagnet, which would give a bell shaped decay. 

As we will explain in chapter [TTI| this ultimate crossover to a Lorentzian can be explained 
by taking into account the long-range dipolar interaction. Further evidence of dipolar effects 
have been found in EuS, where it was observed that the peak positions (in constant-E scans) 
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do not scale l31| . 



Concerning the line width the experimental situation is as follows. Right at the critical 
temperature there is almost perfect agreement of the wave vector dependence of the line 
width with T ~ q b ' 2 in EuS |IJ, EuO |57|,|28|p9| , Fe PI MM and Ni In early 

experimental studies on Fe it seemed that the experimental data PUREUO] are in reasonable 



agreement with the theoretically predicted scaling function of Resibois and Piette [|231 



Recent neutron scattering experiments, however, showed large deviations from the Resibois- 
Piette scaling function in Fe |188| , |189| , EuO ||193| . |192| and EuS ||33|| . This puzzling situation 
can only be resolved by additionally taking into account the dipolar interaction, which is 
the subject of chapters |T| and |TV[ 



2. Ferromagnetic Phase 



The critical dynamics below the transition temperature has been studied also by renor- 



malization group methods. Ma and Mazenko [|174|| calculated the transport coefficient for 
the longitudinal magnetization for small wave vectors in an e— expansion (e = 6 — d). Their 
result was 



tv \ r (?) 2 
r (g) = ~u^ q 



(2.80) 



with 



d-6 

T oc q 6 



(2.81) 



With x L (q) 

oc - in d = 3 dimensions this would give T(q) oc in contradiction to the mode 



coupling result for small q (i.e. in the hydrodynamical limit ||231|| ). However, Sasvari's ||237|| 



reanalysis of Ma and Mazenko's [ 174|| exponentiation method showed, that taking into ac- 
count the regular parts of T(q) results in 



f (?) oc q 



d-6 
3 



^2.82) 



30 



This leads to T(q) oc q 2 for d = 3 dimensions, in agreement with the mode coupling result in 
Table |ITI|. The q A — dependence of the transverse transport coefficient in the hydrodynamical 
limit is also confirmed by the renormalization group calculations ||174|| . 

A thorough renormalization group study of the critical dynamics of a Heisenberg ferro- 
magnet below T c is still lacking. Such a study would have to take into account all the pecu- 
liarities resulting from the presence of the Goldstone modes below T c . As a first step towards 
this end, there is a recent study ||250|1 of the critical dynamics of the 0(n)-symmetric relax- 
ational models with either non-conserved (model A) or conserved order parameter (model 
B) below the transition temperature (see also chapter [V|). 
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TABLES 

TABLE I. Crystal structure dependent parameters of cubic Bravais lattices, c counts the num- 
ber of next nearest neighbors to a given lattice site. The parameter b is defined as b = (a 3 /va) 1 ^ 2 , 
and characterizes the lattice structure. 5 is the distance between nearest neighbors ions and v a the 
volume of the Bravais lattice primitive cell, a is the cube edge. 



sc bcc fee 

c 6 8 12 

b 1 y/2 2 

5/a 1 ^3/2 y/2/2 

v a /a 3 1 1/2 1/4 



TABLE II. Experimental values for the spin-wave amplitude b = //5.1326. 





Fe 


Ni 


Co 


EuO 


EuS 


s 


1.5(1) 


1.5(1) 


1.6(2) 


1.3(2) 


1.4(2) 




1.8(1) 


2.1(1) 






1.9(3) 



Data are collected from Refs. [29, 26,30,34, 220 ), and references cited therein 



TABLE III. Asymptotic behavior of the scaling functions below the critical temperature in 
units of the value at criticality 7(0). 



7 (x)/ 7 (0) 


r(ff) 


ReX(x) 
7(0) 


Im\(x) 
7(0) 


ImA(q); ReA(q) 


x»l 1.37 + 0.37|^ 
x<l 1.0 + 0.55a; 


g 2 r l/2 
q 5/2 


0.16x~ 3/2 ln(x) 
1.0 - 1.34a; 


-0.07x~ 3 / 2 ln(x) 
0.77X 1 / 2 
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Figure Captions: 

Figure 2.1: 

The macroscopic domain of wave vector q and correlation length £. In the three shaded 
regions the correlation functions have different characteristic behaviors, (a) Hydrodynamic 
regions: g( < 1, T > T CI and T < T c , (b) critical region: g£ 1, T « T c . There is 
a change-over from under-damped spin-waves to spin diffusion when the temperature is 
raised from below T c to above T c , as schematically indicated in the diagram. 

Figure 2.2: 

The Resibois-Piette scaling function versus x = l/q£, resulting from the numerical solu- 



tion [231] of the mode coupling equations in Lorentzian approximation for a Heisenberg 



ferromagnet in the paramagnetic phase. 
Figure 2.3: 

The universal function <j)(uj/aq 5 ^ 2 , l/q£) for several values of l/g£: A) l/q£ = 0., B) l/q£ 
0.2, C) l/q£ = 1.04. Taken from Ref. JTT81 . The scale a is defined in Ref. p8| . 



Figure 2.4a: 

The spin relaxation function (j)(Aq 5 ^H, l/g£) for several values of l/q£ indicated in the graph. 
Figure 2.4b: 

The scaling function <fi(uj / Aq 5 ^ 2 , l/q£) for the spin relaxation function u) = 

(j)(u I Aq 5 ^ 2 , l/g£)/Ag 5//2 for several values of l/q£ indicated in the graph, showing the shape 
crossover to a Lorentzian shape for l/q£> 1. 

Figure 2.5: 
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Dynamic scaling function for ferromagnets with short range exchange interaction only versus 
below T c for three different amplitudes of the spin- wave frequency: / = 7r 3//2 (solid), 
/ = 7.65 (dashed), and / = 8.20 (point-dashed). Taken from Ref. JM|>B- 



Figure 2.6: 

Comparison of the longitudinal line width with the polarized neutron scattering experiments 
on Ni p4] . All line widths are normalized to 1 at criticality. The result of mode coupling 



theory is shown for two different values of the spin-wave frequency amplitude / = 7r 3 / 2 
(solid), / = 7.65 (dashed), and / = 8.20 (point-dashed). Taken from Ref. |53p]| 



Figure 2.7: 

Comparison of the line width obtained from mode coupling and renormalization group anal- 



ysis. Taken from Ref. [125 
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FIGURES 




Figure 2.1: The macroscopic domain of wave vector q and correlation length £. In the 
three shaded regions the correlation functions have different characteristic behaviors, (a) 
Hydrodynamic region: ^ « 1, T > T c , and T < T c , (b) critical region: ^ > 1, T 
T c . There is a change-over from under-damped spin-waves to spin diffusion when the 
temperature is raised from below T c to above T c , as schematically indicated in the diagram. 
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Figure 2.2: The Resibois-Piette scaling function versus x = resulting from the 



numerical solution [231] of the mode coupling equations in Lorentzian approximation for a 
Heisenberg ferromagnet in the paramagnetic phase. 
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Figure 2.3: The universal function <j)(uj/aq 5 ^ 2 , l/q£) for several values of A) l/q£ = 0., 
B) l/q£ = 0.2, C) l/qi = 1.04. Taken from Ref. |TT8| . The scale a is defined in Ref. |ffTS |. 
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Figure 2. 4. a: The spin relaxation function 0(Ag 5,/2 t, l/g£) for several values of l/q£ indi- 
cated in the graph. 
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Figure 2.4.b: The scaling function 0(cj/Ag 5//2 , for the spin relaxation function 

a;) = 0(c<j/Ag 5 / 2 , l/g£)/Ag 5 / 2 for several values of l/g£ indicated in the graph, showing 
the shape crossover to a Lorentzian shape for l/q£ > 1. 
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Figure 2.5: Dynamic scaling function for ferromagnets with short range exchange inter- 
action only versus l/g£ below T c for three different amplitudes of the spin- wave frequency 
/ = vr 3 / 2 (solid), / = 7.65 (dashed), and / = 8.20 (point-dashed). Taken from Ref. gggggg 
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Figure 2.6: Comparison of the longitudinal line width with the polarized neutron scattering 
experiments on Ni [Q. All line widths are normalized to 1 at criticality. The result of 
mode coupling theory is shown for two different values of the spin-wave frequency amplitude 
/ = 7T 3 / 2 (solid), / = 7.65 (dashed), and / = 8.20 (point-dashed). Taken from Ref. pg^Hf . 
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III. DIPOLAR FERROMAGNETS 



In this chapter we review the static and dynamic critical behavior of dipolar ferromagnets, 
i.e., spin systems with both short-range exchange and long-range dipole-dipole interaction. 
Special emphasis is put on the discussion of the mode coupling theory in the paramagnetic 
phase. 

But, before turning to the detailed analysis we would like to emphasize the following 
characteristic features of the dipole-dipole interaction, which have important implications 
on the critical dynamics. (1) In contrast to the short range exchange interaction the dipolar 
interaction is long-ranged and thus dominates the asymptotic critical behavior of ferromag- 
nets. (2) It introduces an anisotropy of the spin fluctuations longitudinal and transverse to 
the wave vector q. This implies that the longitudinal static susceptibility remains finite for 
q — > and T —>■ T c [|J. (3) The order parameter no longer is conserved as can be inferred 
from the equations of motion. (4) The strength of the dipolar interaction introduces, besides 
the correlation length £, a second length scale g^ 1 , where q D is the so called dipolar wave 
vector defined below. This leads to generalized scaling laws for the relaxation functions and 
the line widths. 



Our starting point is a Hamiltonian for a spin system including both isotropic short-range 
exchange and long-range dipolar interactions 



where Jw denotes the short range exchange interaction, usually restricted to nearest or next 
nearest neighbors, and Af t , is the dipolar interaction tensor given by 



A. The Model Hamiltonian 




(3.1) 




(3.2) 
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Here g^ is the Lande factor, //g the Bohr magneton. As shown by Cohen and Keffer |47 
the lattice sums 

Af = J2 A m^ (3.3) 



can be evaluated by using Ewald's method ||68|| , and one finds for infinite three-dimensional 
cubic lattices g^iH 



« 2 g 2 - MO 2 5 Q/3 + 0(g 4 , (g Q ) 4 , (O V) 2 ) \ , (3-4) 



where t> a is the volume of the primitive unit cell with lattice constant a, and ctj are con- 
stants, which depend on the lattice structure (see Table [TV]). Upon expanding the exchange 
interaction, 

J q = E'^oe iq - X; « Jo - J<?V + 0(? 4 ) , (3.5) 
j 

and keeping only those terms, which are relevant in the sense of renormalizat ion-group 

theory this results in the following effective Hamiltonian for dipolar ferromagnets 

7 a„/r 



q 



-J + Jq 2 a 2 + Jg^—Tr 



SVJ , (3.6) 



where the Fourier-transform of the spin variables is defined by 

^ = ^E^V q ' Xi . (3.7) 

Here we have defined a dimensionless quantity g as the ratio of the dipolar energy (<7z//b) 2 /a 3 
and exchange energy 2 J, multiplied by a factor 47ra 3 /t> a , which depends on the lattice 
structure. 

Ana 3 (g L /i B ) 2 /a 3 Dipolar Energy 

g = oc — — . (3.8) 

v a 2 J Exchange Energy 

Strictly speaking there are dipolar corrections of order 0(q 2 ) to the exchange coupling. But, 
those can be neglected, since the strength of the dipolar interaction is small compared with 
the exchange interaction. 
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The dipolar interaction induces an anisotropy of the static as well as the dynamic spin- 
spin correlation function with respect to the wave vector q. It was shown by Aharony and 
Fisher [T4J, that the static transverse susceptibility diverges with the dipolar critical exponent 
7 as the critical temperature is approached, whereas the longitudinal susceptibility remains 
finite. The finite value is inversely proportional to the strength of the dipolar interaction 
g. The dipolar anisotropy becomes substantial when q 2 + £~ 2 < q 2 D , where the dipolar wave 
vector is defined as g = (q D a) 2 . The strong suppression of the longitudinal fluctuations have 
been observed in EuS and EuO by polarized neutron scattering experiments [| 1 5 3 1 . 

The matrix of the static susceptibility is given by 

x aP (g, q) = x T (g, q) (fi°* - tpj + x ^ g , ff) «V ; (3 9) 

where one often uses the Ornstein-Zernike forms for the longitudinal and transverse suscep- 
tibilities 



X T (9,q)= 2 J._ 2 , (3-10) 



X L ^g) = ga + qJ +e - a » ( 3 - H ) 
with the non-universal static amplitude 

Conventional mode coupling theory does not account for effects of the critical exponent r], 
which will be neglected in the following^. Here £ = £o f T ~ Te ) is the correlation length. 
The static crossover from Heisenberg to dipolar critical behavior is partly contained in £ 
through the effective temperature dependent exponent pl| , |84] , p9| v = 7 e g/2. The full dipolar 



201 


41 


42 


,84 


89 



5 For a refined version of mode coupling theory which allows for a consistent treatment of the 
Fisher exponent r/ see chapter VI. 
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The tensorial structure of the static susceptibility suggests a decomposition of the spin 
operator S q into one longitudinal and two transverse components with respect to the wave 
vector q, i.e., 



S q = 5^ + ^ 1 t 1 (q)+^"t J (qO > 



where the orthonormal set of unit vectors is defined by 



(3.13) 



- q £i, M qxe 3 

Q (qf + qi) 



t 2 (q)=qxt 1 (q) 



(3.14) 



For vanishing components of q the limits are taken in the order of increasing Cartesian 
components. 

From the Hamiltonian, Eq. ( |3.6| ), one deduces the following microscopic Heisenberg equa- 
tions of motion JT5| , |5U | 



d s - = Ja2 i 10l{ q ■ (2k " q) ( V^+^{s?-k, st} + A*{S?- k , s?}) 



dt 



(3.15) 



d oTx 



-Ja 



J 3 k 1 



(2tt) 3 k 
k 2 (faq - k 2 



q(2k-q) 



-fa fa 



1 k 2 q 



q — k | Jkf + k 



rrt2 q-m i ^ ( r qL qL\ J0T2 qT 2 x 

: t°q-k5 °k J 1* o TT 777 V L°q-k5 °k J t°q-k> °k J 



+3 



/c-^ "l - 



{^q-kJ ^k) + 



2 I q- k 

fa(faq - k 2 ) 



q — k I Jk\ + k 



(3.16) 



and for S^ 2 correspondingly, where {, } denotes the anticommutator. The terms proportional 
to g, resulting from the dipolar term in the Hamiltonian remain finite as the wave vector q 
tends to zero, whereas all the other terms vanish in this limit. This is responsible for the 
fact that the dipolar forces lead to a relaxational dynamics in the limit of long wavelengths, 
i.e., the order parameter is no longer a conserved quantity. 



B. Mode Coupling Theory for the Paramagnetic Phase 
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1. General Mode Coupling Equations 



As we have explained in chapter II, mode coupling theory amounts in an factorization 
approximation for the transport coefficients 



r L (g,g,t) 



x L (q,g) 



(S L (q,t),S L (q,0)), 



(3.17) 



r r (g,^,t)=^^(^(q,t),^(q,0)) 

x {q,g) 

= ^£(^(q,t),^(q,0)) 



(3.18) 



The mode coupling equations resulting from considering two-mode decay processes have 
been derived in Refs. JF^j3(| ( a = L, T) and are given by 

ksT f v a d 3 k 



T a (q,g,t) = 2(2 Ja 



2\2 



X 



(3.19) 



Here the fc-integration runs over the first Brillouin zone (BZ). The vertex functions for 
the decay of the mode a into the modes (3 and a are given by |3(J 

1 2 



v TT (k, q, g, rf) = 2q cos 9 



k cos 9 1 

q ~ 2 



(3.20) 



v L T( k ,<l,9,v) = 2g sin 



fc cos 1 g 
~~ g 2 + V 



(3.21) 



v TT (k,q,g,ij) = q sin 



cos 6* 1 

g ~ 2 



i + 



2| q-k 



|2 ' 



(3.22) 



vl T {k,q,g,r}) = q 4 



k cos 9 1 g 
~~q 2 + V 



1 2 



q — k 



-5 sin^ , 



(3.23) 



v LL (k,q,g,r]) = q sin 



/c cos 6 1 1 
9 ~ 2 



2| q- k 



|2 ' 



(3.24) 
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with i] = cos 9. In passing we note that there were certain attempts to develop a mode cou- 
pling theory already twenty years ago, although nobody succeeded to derive the appropriate 
mode coupling equations. A type of mode coupling calculation was used by Huber ||122| 



to determine the uniform spin relaxation for temperatures larger than the dipolar crossover 
temperature. This analysis was extended by Finger ||72|| , who put forward certain scaling es- 
timates and computed the uniform spin relaxation in the strong dipolar region. An attempt 
to construct a mode coupling theory was launched by Borckmans et al. [P7|| , using an incom- 
plete basis and ending up with equations containing undetermined vertices. The theoretical 
and experimental situation was rather controversial and no explanation was available for the 
apparent contradictions ||151| , |191|| . Only in 1986 a complete self-consistent mode coupling 



theory was developed by the authors |79|,|80| and its various properties and consequences were 



studied in detail |8T|,|82|j8^,85l. Some of the results were confirmed numerically by Aberger 



and Folk and by Kalashnikov and Tret'jakov [p.31|J132|| using analytic approximants. 



The mode coupling result for the transport coefficients, Eqs. ( |3.19| ), together with the 
relation 

-$ a (q, g,t) = -J drT a (q, g, t - r)<S> a (q, g, r) (3.25) 
o 

for the Kubo relaxation functions constitute, as in chapter II, a complete set of self consistent 
equations. 

As emphasized before, the dipolar interaction introduces a second length scale g" 1 besides 
the correlation length £. This entails the following extension of the static scaling law for the 
spin susceptibilities 



with the scaling variables 



x a (q^,g) = r q - 2 x a (x,y), (3.26) 



x = — and y = —. (3.27) 
« 1 



Note, that here and in the remainder of this section we have explicitly indicated the depen- 
dence of the susceptibility on the correlation length. In all other parts of this review we 
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suppress this dependence for notational convenience. Since the vertex functions v% a are pro- 
portional to the fourth power of the wave number, v* a oc g 4 and because of the homogeneity 
of Eq. (|3.26| ), the relaxation functions and line widths derived from Eqs. ( |3.19| )-( |3.25| ) obey 
the dynamical scaling laws 

$ Q (gZ, gl 2 , ul z ) = r z $ a (q, £, g, u) , (3.28) 

and 

r*(ql,gl 2 ,ul z ) = l t T"(q,£,g,u), (3.29) 

with z = 5/2 and a scaling parameter I. We emphasize that despite of z assuming the 
isotropic value 5/2, there is a crossover to dipolar critical behavior contained in the functional 
form of the correlation functions, as will become clear below. 

An immediate consequence of Eq. (|3.29|) is the following scaling property of the charac- 
teristic longitudinal and transverse frequencies u)"(q,£,g) 

u:(q,Z,g) = Aq z n a (x,y), (3.30) 

where A is a non universal coefficient. 

Now there are various ways to rewrite the scaling laws Eqs. ( 3.28 ) and ( p.29j ) by appro- 



priate choices of the scaling parameter /. If one sets I = q 1 one finds 



$ Q (g,£,^) = g- 2 $ a (l 4,^), (3.31) 

q£ q l q z 



and 



V a {q^,g,u) = q z Y a {- ^^). (3.32) 



A disadvantage of the representation given in Eqs. fl3.31|JT3~2| ) is that both the crossover of 
the time scales and of the shapes of the correlation functions are intermixed in $ Q . Since 
the time scales for the isotropic and dipolar critical and hydrodynamic behavior differ quite 
drastically, it is more natural to measure frequencies in units of the characteristic frequencies. 
Hence we fix the scaling parameter by the condition 
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and find from Eqs. ( |3.3iy3.32| ) the scaling forms 



and 

T a (q, £, g, to) = Aq z n a (x, y) 7 a (x, y, v a ) , (3.35) 
with the scaling variable for the frequency 

v a = i — 77-7 r • (3.36) 

Aq z Q a (x,y) K ' 

With Eq. ( |3.33| ) one has separated the crossover of the frequency scales and the crossover 
of the shapes of the correlation functions. The former mainly is contained in Q a (x, y) the 
latter in <p a (x, y, u a ). 

There is still some freedom in the choice of the characteristic frequencies u;° in Eq. Q3.33|) ; 
for instance, one could take the half width at half maximum (HWHM) of the frequency 
dependent Kubo functions. This, however, would require to solve Eqs. Q3.19Q and (|3.25|) 
simultaneously for the time scales and the shapes of the correlation functions. Therefore, 
in the following we will use as characteristic frequencies the half widths resulting from the 
Lorentzian approximation for the line shape (see section 3.3). The Lorentzian line widths 
qualitatively obey the same scaling laws as the HWHM and have the same asymptotic (hy- 
drodynamic, dipolar, isotropic) properties. Thus this choice for the characteristic frequencies 
solely is a matter of numerical convenience and does not introduce any approximations in 
the final result. From the final result one can obtain the HWHM and rewrite the scaling 
functions in terms of these new variables. 

Eqs. ( |3.34| ), (|3.35f) and (|3.36|) imply for the Laplace transformed quantities the scaling 



laws 

<5> a (q,Z,g,t) = <f ) a (x,y,T a ), (3.37) 
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and 

T a (q, f , g, t) = [Aq z n a (x, y)f 7 a (x, y, r a ) , (3.38) 

where the scaled time variables r a are given by 

r a = Aq z tt a (x,y)t. (3.39) 

One should note that the characteristic time scales l/[Aq z Q a (x, y)\ are different for the lon- 
gitudinal and transverse modes. This is mainly due to the non critical longitudinal static 
susceptibility implying that the longitudinal characteristic frequency Aq z Q L (x,y) shows no 
critical slowing down asymptotically. In other words for T = T c and q — > the longitudi- 
nal characteristic frequency does not tend to zero, which implies an effective longitudinal 
dynamical critical exponent z^ s = for the wave vector dependence in the limit q — > 0. In 
comparison, the effective exponent for the transverse characteristic frequency at T c shows a 
crossover from zj s = 5/2 to zj s = 2 (see also the following section). This mode coupling 
result disagrees with a calculation based on nonlinear spin wave theory, where zj s = 1 is 
found in the dipolar region | |176 |. 



Inserting Eqs. ( |3.37| , |3.3g| ) together with the static scaling law ( |3.26p into Eqs. ( |3.19j ) and 



( |3.25| ) we find the following coupled integro-differential equations 

/ \ 2 +1 Pcut 



7 a (x,y,r Q ) 



7T 



fi Q (x,3/), 



-i n /V 



and 



X a (x,y) \p p J \p- p- 



d 

(j) a (x, y, r a ) = ~ I dT^ a (x, y, r a - r)(j) a (x, y, r) , (3.41) 



u o 



connecting the scaling functions for the transport coefficients with the scaling functions 
for the Kubo relaxation functions. In Eq. (|3.40|) we introduced the notations p = k/q, 
p- = | q- k | /q, 7} = cos(q, k) and T a/3 (x,y,p) = t^QP (x/p,y/p) /Q a (x,y). 
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The non universal frequency scale resulting from the transformation of Eq. ( |3.19|) in 
Eq. ( ggDp is 



A 



(3.42) 



b V 4vr 4 

The apparent critical dynamic exponent contained in Eq. ( |3.40|) equals z = 5/2. However, 
as noted before, the crossover to dipolar behavior is contained in the scaling functions for 
the transport coefficients 7 a (x, y, r a ), the scaling functions for the Kubo relaxation functions 
4> a (x,y,r a ) and the scaling functions of the characteristic frequencies fl a (x, y). 
The scaled vertex functions appearing in Eq. ( |3.40|) read |8(| 



v 



2 V 2 5 a > L + (l - V 2 ) (V' T + ^ 5 a ' T 



pv-- 2 



v 



LT 



2(l- v 2 5 a > L ) - (l - r7 2 ) ( 1 



1 

7- 



6 



PV 



1 y 2 
- + — 

2 2 



(3.43) 



(3.44) 



which are related to the vertex functions v% a of Eqs. ( |3.20| - |3~2"4"|) by Vg a = q A v ( p (T - For both 
longitudinal and transverse modes, the dipolar interaction enters only in vertices involving 
decays into a longitudinal and a transverse mode, since the dipolar interaction enters the 
Hamiltonian only through the longitudinal modes. 

Because the /c-integration is restricted to the Brillouin zone the p-integration of Eq. ( |3.40| ) 
contains the cutoff 



Pcut 



Qbz Qbz 



-y- 



(3.45) 



where qsz denotes the boundary of the first Brillouin zone. All other material dependent 
parameters are contained in the frequency scale A. The cutoff is important for small times, 
because the integrand of Eq. ( |3.40| ) is of order 1 for t = and p ^> 1. Hence for small times 
also wave vectors near the zone boundary contribute to the relaxation mechanism. 



As explained before, in the numerical solution of the MC-equations one has taken p2|j85 
for the characteristic frequencies the line widths resulting from the Lorentzian approximation 
of the MC equations, i.e., 
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W(x,y)=T&(x,y). (3-46) 

Using as input the solution of the mode coupling equations in the Lorentzian approximation 
one can solve the complete set of MC-equations for different values of qsz- Because there 
are three scaling variables (x, y and u a ) it is impossible to present here all numerical results. 
Instead, in chapter IV, we will concentrate on a limited number of temperatures and wave 
vectors motivated by the available experiments on the substances of primary interest EuS, 
EuO and Fe. 



2. Lorentzian Approximation 
For later reference we want to close this section with quoting the results from the so 



called Lorentzian approximation |Hyj8T|. The results of the numerical solution of the full 
mode coupling equations will be discussed in chapter IV in conjunction with the experimental 
data. 

If the transport coefficients vary only slowly with uj we may approximate the relaxation 
functions by Lorentzians, i.e., we replace the transport coefficients by their values at u = 

rf or (g, f , g) = T L (q, £,g,u = 0), T^q, f , g) = T T (q, £,g,u = 0). (3.47) 

Despite a Lorentzian being not the correct shape of the correlation function for all values 
of the scaling variables, the resulting line widths obtained in the Lorentzian approximation 
already capture most of the crossover in the time scale. 
The Lorentzian line widths obey the scaling law 

K r (q,^g) = Aq z lZr(^y)- (3-48) 

From Eq. ( |3.20|) it is then easily inferred that the scaling functions of the Lorentzian line 
widths 7ior( a; ) y) are determined by the coupled integral equations 

+1 oo 

V"(.r. //) 



7l ° r(x,2/) = yHx v) I d7] I ^p-^^lAyip^) 



-j o P u 
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As summarized in Table 3.1 the mode coupling equations ( p.49| ) can be solved analytically in 
the dipolar (D) and isotropic (I) critical (C) and hydrodynamic (H) limiting regions. These 
are defined by DC: y > 1, x < 1; IC: y < 1, x < 1; DH: y > x, x > 1; IH: y < x, x > 1. 

Concerning the critical dynamical exponent one finds for the longitudinal line width a 
crossover from z = 5/2 in the isotropic critical region to z = in the dipolar critical region, 
whereas for the transverse line width the crossover is from z = 5/2to;z = 2. The precise 
position of this crossover can only be determined numerically. 

The numerical solution |79| , |80|| of Eq. ( |3.49|) shows that the dynamic crossover for the 
transverse width is shifted to smaller wave vectors by almost one order of magnitude with 
respect to the static crossover, whereas the crossover for the longitudinal width occurs at the 
static crossover. For the numerical solution of the mode coupling equations it is convenient 
to introduce polar coordinates 

r = \J x 2 + y 2 , and p = arctan(y/a;) . 

The transverse and longitudinal scaling functions 7 T (r, ip) and 7 L (r, (p) are shown in Figs. 3.1 
and 3.2 as a function of the radial scaling variables r and ip. A different representation of 
the results can be given by plotting the linewidth versus the single scaling variable x for 
several values of p. This is shown in Figs. 3.3-5, where we have drawn the two-parameter 
scaling functions 'y^ix, y) in units of the value at criticality 70 = 7ior T ^(0> 0) — 5.1326. The 
physical content of the two parameter scaling surfaces is illustrated best by considering cuts 
for fixed q D and various temperatures. In Figs. 3.3,4 the scaling functions versus x = 1/qC, 
are displayed for different values of <p = arctan(g D ^) = Nn/20 with N = 0,1,..., 9. For 
tp — 0, corresponding to vanishing dipolar coupling g, the scaling functions coincide with 
the Resibois-Piette scaling function. If the strength of the dipolar interaction q D is finite, 
the curves approach the Resibois-Piette scaling function for small values of the scaling 
variable x and deviate therefrom with increasing x. For a given material, q D is fixed and the 
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parametrisation by <p corresponds to a parametrisation in terms of the reduced temperature 
(T-T c )/T c . 

To examine the dipolar crossover precisely at the Curie point, Fig. 3.5 displays the scaling 
function for the transverse and longitudinal width at T = T c versus the wave number, i.e., 
y^ 1 = q/q D . This graph shows that the crossover from the isotropic Heisenberg to dipolar 
critical dynamics in the transverse line width occurs at a wave number, which is almost 
one order of magnitude smaller than the static crossover wave vector q D . The crossover 
of the longitudinal width, from z = 2. 5 to z = 0, is more pronounced and occurs in the 
intermediate vicinity of q D . The reason for the different location of the dynamic crossover is 
mainly due to the fact that it is primarily the longitudinal static susceptibility which shows 
a crossover due to the dipolar interaction. Since the change in the static critical exponents is 
numerically small the transverse static susceptibility is nearly the same as for ferromagnets 
without dipolar interaction ||. Hence the crossover in the transverse width is purely a 
dynamical crossover, whereas the crossover of the longitudinal width being proportional to 
the inverse static longitudinal susceptibility is enhanced by the static crossover. In order to 
substantiate these arguments we have plotted in the inset of Fig. 3.5 the scaling functions 
of the Onsager coefficients SCl a a t the critical temperature versus q/q D , showing only the 
dynamical crossover. 



3. Selected Results of the Complete Mode Coupling Equations 

All the scaling functions for the dipolar ferromagnet depend on the three scaling variables 
x — l/<?£> V = Qd/Q an d v a = uj/ Aq z Q a (x,y). Therefore it is impossible to present all the 
results obtained from the numerical solution of the complete mode coupling equations. In 
this section we intend to review the most important features of the shape and line width 



crossover reported in Refs. p2| , p5| . Further results and discussion will be presented in the 
next chapter in conjunction with a comparison with experimental data. 

The results from the Lorentzian approximation show that the dipolar line width crossover 
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in the transverse line width sets in at a wave vector almost one order of magnitude smaller 
than the dipolar wave vector q D . In order to get information about the line shape, one 
has to dismiss this approximation and solve the complete set of mode coupling equations, 
Eqs. ( |3.40D and ( |3.41| ). This has been achieved in Ref. p2| , p5| , and one finds the following 



crossover scenario. First of all, the crossover in the line shape sets in at wave vectors of the 
order of the dipolar wave vector. 

Let us first consider the case of temperatures very close to T = T c . Fig. 3.6 and 3.7 show 
the transverse and longitudinal scaling functions <f) a (r,ip,T a ) versus the scaling variables r 
and r a for ip = 1.49. Referring to EuO, characterized by the (non universal) parameters q D = 
0.147A -1 , T c = 69. IK and £ = 1.57A, this corresponds to a temperature T ps (1 + 0.003)T c . 
The line shapes of the longitudinal and transverse relaxation function agree in the isotropic 
Heisenberg limit, i.e., for r — > corresponding to large values of the wave vector q (q ^> q D ). 
In this limit the dipolar interaction becomes negligible and the shape is of the Hubbard- 
Wegner type as discussed in section II. Upon increasing the value of the scaling variable 
r the line shape of the transverse and longitudinal relaxation function become drastically 
different. Whereas the transverse relaxation function shows a nearly exponential decay, 
pronounced over-damped oscillations show up for the transverse relaxation function. The 
shape crossover is also shown in Fig. 3.8, where the transverse and longitudinal relaxation 
functions are plotted versus the scaling variable r a for three different values of the scaling 
variable r = ^/x 2 + y 2 (r = 0, r = 1, and r = 10). 

For temperatures well separated from T c , i.e. q D ^ <C 1, the dipolar interaction becomes 
negligible. Hence, the difference in the shape crossover of the longitudinal and transverse 
relaxation function diminishes with decreasing q D £. For q D £ <C 1 the shape crossover as 
a function of r corresponds to the crossover from the critical (Hubbard- Wegner) shape to 
the hydrodynamic shape as discussed in section II (Note that for q D ^ <C 1, the scaling 
variable r reduces to x = l/q£,.). Figs. 3.9 and 3.10 show the shape crossover for q D £ = 3.52 
(if = 1.294) for the transverse and longitudinal relaxation function, respectively. Figs. 3.11 
and 3.12 display the corresponding crossover for q D ^ = 1.46 (ip = 0.97). In the latter 
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case the shape crossover of the transverse and longitudinal relaxation function are already 
quite similar and almost identical to the shape crossover found for the isotropic Heisenberg 
ferromagnet without dipolar interaction. The various crossover scenarios can also be read 
off from Fig. 3.8. 

Finally let us add a comment on how the line shape crossover affects the line width 



crossover. It has been shown f82|,|85j (see Section IV) that there are only slight changes 
of the line width crossover, when the line width is determined from the full solution of 
the mode coupling equations, as compared to the widths obtained from the Lorentzian 
approximation. Roughly speaking, the overall effect of taking into account the correct line 
shape is approximately a shift by a constant factor. For details we refer the reader to section 
IV. 

C. Spin Wave Theory in the Ferromagnetic Phase 

Holstein and Primakoff [ |1 1 5|1 have investigated the dynamics of dipolar ferromagnets far 
below T c using linear spin-wave theory. Upon neglecting fluctuations in the longitudinal 
component, Sf ~ S, they get the following equations of motion for the transverse spin 
fluctuations (see also Ref. ||168|| ) 



ip+^A^ + B^, (3.50) 



ij f S k = AiSZ^ + B k S+ , (3.51) 



in terms of the raising and lowering operators <Sj = Sf ± iSf. From these equations of 
motion the spin-wave dispersion follows (see also Ref. ||141|| ) 



e k = V^k+I^kl V^k-|£k|. (3.52) 

The coefficients A k and B k are given in terms of the exchange interaction J k = 
2SJ2m ie' k ' (xrXm) and the dipolar interaction tensor A^f ( see Eq. (pO)) ) 

A k = g L fJ. B Ho + (Jo - Jv) + (2£ + E v ) , (3.53) 
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Eu= SA^ , 

5 k = -S {A™ - Al y - 2iA^) = SyJ(A%° - A vy ) 2 + 4^V 2 ^ k 



(3.54) 
(3.55) 



For crystals with cubic symmetry the dipolar tensor A%r is given by A%. ~ 
\{9L^B) 2 ^f (S a/3 — ^pr-) near the zone center ka <C 1, plus small structure dependent terms 
proportional to k 2 (see Section III. 1) ) . Note however that the dipolar tensor becomes severely 



structure dependent (requiring numerical evaluation, see Ref. |^7j) as k approaches the zone 
boundary. 

Hence, for ka ^ 1 one finds for crystals with cubic symmetry (for more details we refer 



to the review by Keffer [1411) 



A k = g L ^ B H + (J - Jk) + 7,9lVbM [4vr - 47r(k z /k) 2 



4:7V 



(k x ) 2 -i(ky) 2 
k 2 



27igLjjBM sin 2 6kexp(— 2i(fk) 



(3.56) 
(3.57) 



where 0k and ip^ define the orientation of the wave vector k with respect to the z-axis; Mq 
is the saturation magnetization. Note, that the latter equations are strictly valid for long 
thin samples only. Otherwise, one has to take into account demagnetization effects, which 
amount in an additional magnetic field Hq — > Hq — N z Mq, where N z is a demagnetization 
factor ( jUTj, see also Ref. (TJ). 

Using the formalism of Holstein and Primakoff [115] the influence of the dipole-dipole 
interaction on the inelastic scattering of neutrons has been investigated by Elliott and 
Lowde (see also fll70fl ). 

The neutron scattering cross section is proportional to the correlation function corre- 
sponding to the spin fluctuations (see e.g. Ref. ||166| , |50| ). If one considers N identical 
atoms with fixed positions one finds for the magnetic partial differential cross section for 
unpolarized neutrons in forward direction 

d 2 a 



TTroinQ)| 2 ^£e* Q(x, - x,/) / ' ' ^-'^V/(f),S-(0)) 



k' N 

dfldE' ~~kli 

= jjr 2 \F(Q)\ 2 S T (Q^) 



IV 



+°° dt 
oo 2tx 



(3.58) 
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where Sj 1 ^) is the component of the spin density perpendicular to the momentum transfer 
(scattering vector) Q = k — k', and hu is the neutron energy transfer. The length r = 
— 5.391/m is analogous to the nuclear scattering length (see also chapter IV), and F(Q) 
the magnetic form factor. We have also defined the transverse magnetic scattering function 
S t (Q,uj) by 



6 al3 



Q a Q l3 \ /■+<» dt 

~Q 2 



)r+oo rlf 
/ ^e"^(^ q (t)^(0)) 



(3.59) 



where Q = q + r with q the wave vector of the magnon and r a reciprocal lattice vector. 

The equations of motion, Eqs. |3.50| and |3.51| , are diagonalized in terms of the Bose 
operators a and a) that satisfy 



[a q , aL] = 8 m < 



(3.60) 



with 



(3.61) 



A convenient choice for the coefficients u q and u q is [141, 168|| f| 



<=(2^iV)(A q + e q )/2eq 

Uq = -W q S*/(Aq + e q ) . 



(3.62) 
(3.63) 



The contributions from single-spin-wave events to the dynamic structure factor was 



calculated by Elliot and Lowde f67| . |17(f 



s t (q,cu) = !k + i 



(Q 



z\2 x 



Q 2 



(Q 



\B r , 



Q 2 



■ C0S[2((^ Q - <^ q )] 



+ Cq) , 

(3.64) 



Holstein and Primakoff [115] have ignored the phase relationship between u q and t> q , without 
which an incorrect expression for the scattering cross section between spin-waves and neutrons 
would be obtained. This error has been corrected by Lowde [37|. 170,141|. 
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where n q is the occupation number of the magnon oscillator with wave vector q. The 
second term in the dynamic structure factor is due to the dipolar interaction. For zero 
dipolar interaction or for wave vectors larger than the dipolar wave vector q D the dynamic 
structure factor reduces to a simple angular distribution proportional to 1 + (Q z /Q) 2 - When 
dipolar effects become of importance, the angular dependence of the scattered intensity gets 
quite complicated. 

More recently Lovesey et al. |168|j254j] have extended Lowde's analysis to scattering of 



polarized neutrons, and a discussion of the longitudinal cross section (SfSi,), which contains 
two-spin-wave scattering events. Furthermore, they analyze the static susceptibilities in the 
framework of linear spin-wave theory, correcting work by Toh and Gehring [|252|| who used 



the incorrect expressions for the coefficients w q and t> q from Ref. ||115f| . 



As a first step beyond the spin-wave theories described above, Toperverg, and 



Yashenkin ||253|| have used the perturbation approach developed by Vaks, Larkin and 
Pikin ||257|| , to investigate the frequency dependence of the uniform transverse and lon- 
gitudinal susceptibilities []. The applicability of their results is mainly restricted to low and 
intermediate temperatures. Their perturbation approach for the dipolar interaction breaks 
down not only close to the critical temperature, but also at any temperature for low fre- 
quencies. For this parameter regime a non perturbative approach like mode coupling theory 
is needed to account for the strong fluctuations. A first attempt towards such a theory 



has been made some time ago by Raghavan and Huber ||227|| . There are, however, several 



shortcomings in their approach. First of all, the static susceptibilities used in their anal- 
ysis do not account for the coexistence anomalies and dipolar crossover properly. Instead 
the longitudinal susceptibility is taken of Ornstein-Zernike form xl oc l/(<? 2 + £ _2 )> which 
neglects dipolar crossover effects as well as the by now well known coexistence anomaly 



7 An excellent review on the theoretical and experimental work prior to 1984 was given by Maleev 
[17S], with a particular emphasis on Green's function techniques. 
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Xl oc 1/q in the limit q — > 0. The expression used for the transverse susceptibility is valid 
in certain limits only. Furthermore, in the presence of dipolar interaction, for a general an- 
gle between the wave vector and the direction of the spontaneous magnetization, there are 
three and not just two non degenerate eigenvalues for the static susceptibility matrix. This 
fact has been neglected completely. Also, they didn't evaluate the complete functional form 
of the relaxation functions, but used instead a parametrisation, which is strictly valid in 
the hydrodynamic regime only, and calculated the corresponding parameters. Nevertheless, 
the theory seems to give a quite reasonable description of the data obtained by a neutron 
scattering study on EuO [j207|| in the range g£ < 1 and q D C, < 1 for not too small values of 
the wave vector q. Beyond this range the approximate treatment of the dipolar interaction 
in Ref. ||227|| breaks down. 



More recently Lovesey ||169|| reported on an approximative mode coupling approach for 



dipolar ferromagnets below T c . Similar to Raghavan and Huber the spin wave dispersion 
relation used by Lovesey is applicable for not too small wave vectors only. The analysis in 
Ref. ||169|] is restricted to the exchange region q D ^ < 1. This is due to the assumption made in 



Ref. [ |169| ] , that the only non vanishing relaxation kernels are those for the spin fluctuations 
longitudinal and transverse to the direction of the spontaneous magnetization. Neglecting 
off-diagonal matrix elements excludes the applicability of the theory to the dipolar region. 
Note, that in the dipolar region close to T c the memory function and the relaxation functions 
become - similar to the situation above T c - diagonal in terms of the spin fluctuations 
longitudinal and transverse with respect to the wave vector and not to the spontaneous 



magnetization. Furthermore, the analysis in Ref. [16S] is restricted to very small wave 



vectors. In this limit, however, the expressions for the static susceptibility and the spin- 



wave dispersion relation used in Ref. ||169|| are not valid, since they do not account for the 
presence of the subtle combined effects of Goldstone modes and dipolar anisotropy. Also, 
no self-consistent solution but only a first iteration of the mode coupling equations, based 
on neglecting the damping, is performed. In summary, a thorough mode coupling analysis 
of the effects of the dipolar interactions in the ferromagnetic phase is still a very challenging 
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theoretical problem and the topic of ongoing research 
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D. Renormalization Group Theory of Time Dependent Ginzburg Landau Models in 

the Ferromagnetic Phase 

The spin-wave theory, reviewed in the preceding section, is only a first step towards 
a more rigorous theory of the critical dynamics of dipolar ferromagnets below the Curie 
temperature. For a thorough understanding of the dynamics in the ferromagnetic phase a 
renormalization group theory or a mode coupling approach analogous to section III.B, which 
takes into account the effects of the critical fluctuations, would be necessary. A detailed 



analysis requires a treatment of a modified model J ||112|| appropriate for the dynamics of 
isotropic ferromagnets, where dipolar forces are included. 

Recently, the effects of the dipolar interaction on the critical dynamics of the n- 
component time-dependent Ginzburg-Landau models (model A |112| ) below the critical 
temperature have been studied within a generalized minimal subtraction scheme [|251|| . The 
corresponding Langevin equation of motion reads (Model A: a = 0; Model B: a = 2) 

^-^)« + «->< 

where the stochastic forces are characterized by a Gaussian probability distribution function 
with zero mean and variance 

(C a (x, t)C V, 0) = 2AA: B T(zV) a 5 (3) (x - x')5(t - t')8^ . (3.66) 

The Ginzburg-Landau effective free energy functional reads 

rmin(d,n) 

H[{S a }] = -j E [(r + q 2 )P^ + (r + g + q 2 )P^]S a (q)S^-q) 



a,/3=l 



E S a ( q )S a (-q) 

a=min(d,n)+l 



(3.67) 



where the general situation n ^ d is considered. Here P^p = 5 a @ —q a qP /q 2 and P^g = q a q /q 2 
denote the transverse and longitudinal projectio operator. 
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Those relaxational models neglect mode coupling terms resulting from the reversible 
motion of the spins in the local magnetic field. It is expected ||251|| , however, that most 



of the conclusions based on the relaxational models will also hold for models with mode 
coupling terms. 

As a consequence of the spontaneously broken symmetry there are n — 1 massless Gold- 
stone modes in the ordered phase of ideally isotropic systems. These massless modes lead 
to infrared singularities (coexistence singularities) in certain correlation functions for all 
temperatures below T c . Based on the analysis of the effects of the critical and Goldstone 
fluctuations for the 0(n)-symmetric time-dependent Ginzburg-Landau models ||250|| , it has 



been investigated ||251|| how the coexistence anomalies are modified when dipolar forces or 



weak anisotropies are included. For later reference the coexistence anomalies of the isotropic 
relaxational models are collected in Table VI. 



The analysis in Ref. ||251|| shows that the influence of the dipolar interaction on the 
coexistence singularities is quite subtle. Although the model explicitly breaks the 0{n)- 
symmetry, not all transverse modes lose their Goldstone character, but only their effective 
number is reduced by one. Hence, while for n = 2 a crossover to an asymptotically uncritical 
theory takes place, for n > 3 coexistence anomalies persist, governed by a dipolar coexistence 
fixed point ||251| . 



Below T c there are two preferred axes, the axis defined by the direction of the spontaneous 
magnetization and the axis defined by the wave vector q. This leads to a complex structure 
of the correlation functions already on the harmonic level. It is quite remarkable, however, 
that a one-loop theory for the two-point cummulants becomes an exact representation in 
the ordered phase in the coexistence limit (q — > and u — > 0) p51| . 



For n = d > 3 it is found ||251|| , that the power laws characteristic of the coexistence 
anomalies are not changed by the presence of the dipolar interaction. Hence the same power 
laws as in Table VI apply to the dipolar case also. It is also interesting to note that there is 
the following exact amplitude ratio of the longitudinal response function in the dipolar and 
the isotropic case: 
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XL(q,^)dipolar _ Tl-2 ^ 6g ^ 

XL (^Osteotropic n-1 

For n = 3 this universal amplitude ratio is 1/2, in accord with the results of Pokrovsky ||222| 
and Toh and Gehring ||252|| , obtained in the framework of a spin-wave theory. 

As is apparent from the value of the dipolar coexistence fixed point u* CD = 
6(4 — d)/(n — 2), the situation n = 2 requires a separate discussion. In this case there 
are no massless modes left, since the dipolar interaction reduces the number of Goldstone 
modes from n — 1 to n — 2. As a consequence the crossover below T c is from a critical theory 
at T c to a Gaussian theory, i.e. the fluctuations die out on leaving the critical region. A 



qualitative summary of the various crossover scenarios is given in Table VII [|251 
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TABLE IV. Coefficients in the Taylor expansion of for three-dimensional cubic lattices, 
taken from [Aharony73a, Aharony73b] . c denotes the coordination number, v a the volume of the 
primitive unit cell and ctj are lengths characterizing the dipolar interaction. 



Lattice sc bcc fee 

c 6 8 12 



v a \a 



31 



1 4/3^ 1/V2 



ai[a] 1.2755 1.7420 2.8313 

a 2 [a] 0.1649 -0.321 -0.335 

a 3 [a] 1.7700 0.8210 1.823 



TABLE V. Asymptotic behavior of the scaling functions for the longitudinal and transverse 
Lorentzian line width in the paramagnetic phase. The different regions dipolar critical (DC), 
isotropic critical (IC), dipolar hydrodynamic (DH), and isotropic hydrodynamic (IH) are defined 
by DC: y » 1, x < 1; IC: y < 1, x < 1; DH: y » x, x » 1; IH: y < x, x » 1. 



7 



7 



DC 
IC 
DH 
IH 



,V2 



1/2 



y 5/2 
1 

y 5/2 
x l/2 
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TABLE VI. Coexistence anomalies of the isotropic relaxational models for the longitudinal 
dynamic susceptibility xl(^^) and correlation function Gx(q, u). 



Model A Model B 



a 2 

XL(q,0) oc q d ~ 4 oc q d ~ 4 

Re [ X l(0, co/q a )} oc i^ 4 )/ 2 oc (uj/q 2 )^/ 2 

G L (0,Lo/q a ) ocoj d / 2 ~ 3 ot(uj/q 2 ) d - 4 



TABLE VII. The influence of the dipolar interaction on the coexistence anomalies. The ta- 
ble summarizes the various crossover scenarios possible if the number of components n and the 
dimensionality d of space is varied. 



d=2d=3 d=4 
n = 1 Crossover to a Gaussian theory as for g = 

n = 2 Crossover to a Gaussian theory no anomalies 

n > 3 u cd = 6 ^-2^ ~ * anomalies logarithmic corrections 
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Figures captions: 



Figure 3.1: Scaling function 7 for the transverse width in Lorentzian approximation 
versus r = l/g£(l + (g^) 2 ) 1 ^ 2 an d V 9 — arctan(g D £). 

Figure 3.2: Scaling function 7 L for the transverse width in Lorentzian approximation versus 
r = 1M(1 + (<? D 2 ) 1/2 and V = arctan(g D £)- 

Figure 3.3: Scaling function / y T for the transverse width in Lorentzian approximation 
versus l/q£ for values of (p — Nir/20 with N indicated in the graph. 

Figure 3.4: Scaling function 7 L for the longitudinal width in Lorentzian approximation 
versus for values of ip — iW/20 with N indicated in the graph. 

Figure 3.5: Scaling functions for the transverse (solid) and the longitudinal (point-dashed) 
widths versus q/q D in Lorentzian approximation at the critical temperature. Inset: Scaling 
functions for the transverse (solid) and the longitudinal (point-dashed) Onsager coefficients 
versus q/q D at the critical temperature. 

Figure 3.6: Scaling function of the transverse Kubo relaxation function <f) T (r, </?, 7r) at 
ip = 1.49 (close to the critical temperature) versus r T and r 

Figure 3.7: Scaling function of the longitudinal Kubo relaxation function <f) L (r,(p,T L ) at 
ip = 1.49 (close to the critical temperature) versus Tl and r = \J(q D /q) 2 + (l/q£) 2 . 

Figure 3.8: Scaling function of the longitudinal and transverse Kubo relaxation function 
4> a (r,ip,TT) versus r a for three different values oi ip (ip — 1.49, 1.294, 0.97). In each graph 
the scaling function is shown for r = (solid), r = 1 (dashed), and r = 10 (dot-dashed). 

Figure 3.9: Scaling function of the transverse Kubo relaxation function T (r, </?, tt) at 
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ip = 1.294 versus r T and r = ^ (q D /q) 2 + (l/g£) 2 . 

Figure 3.10: Scaling function of the longitudinal Kubo relaxation function (f) L (r,(p,r L ) at 
ip = 1.294 versus t l and r = <J (q D /q) 2 + (1/ ^C) 2 - 

Figure 3.11: Scaling function of the transverse Kubo relaxation function (f> T (r,ip,r T ) at 
ip = 0.97 versus r T and r = ^/ (q D /q) 2 + (l/q£) 2 . 

Figure 3.12: Scaling function of the longitudinal Kubo relaxation function (j) L (r,(p,r L ) at 
y? = 0.97 versus t l and r = J (q D /q) 2 + (l/g£) 2 . 
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Figure 3.1: Scaling function 7 for the transverse width in Lorentzian approximation 
versus r = l/g£(l + (f^O 2 ) 1 ^ 2 and (p = arctan(g D £). 
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Figure 3.2: Scaling function / y L for the transverse width in Lorentzian approximation versus 
r = 1/^(1 + (q D 2 ) 1/2 and up = arctan(g D 0- 
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Figure 3.4: Scaling function 7 L for the longitudinal width in Lorentzian approximation 
versus l/q£ for values of f = iW/20 with N indicated in the graph. 
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Figure 3.5: Scaling functions for the transverse (solid) and the longitudinal (point-dashed) 
widths versus q/q D in Lorentzian approximation at the critical temperature. Inset: Scaling 
functions for the transverse (solid) and the longitudinal (point-dashed) Onsager coefficients 
versus q/q D at the critical temperature. 
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Figure 3.6: Scaling function of the transverse Kubo relaxation function T (r, ip, vr) at 
V? = 1.49 (close to the critical temperature) versus r T and r 
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Figure 3.7: Scaling function of the longitudinal Kubo relaxation function L (r, ip, t£) at 
<f = 1.49 (close to the critical temperature) versus r L and r = ^ (q D /q) 2 + (l/g£) 2 . 
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Figure 3.8: Scaling function of the longitudinal and transverse Kubo relaxation function 
(f) a (r,(p,T T ) versus r Q for three different values of </?(</? = 1.49, 1.294, 0.97). In each graph 
the scaling function is shown for r = (solid), r = 1 (dashed), and r = 10 (dot-dashed). 
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Figure 3.9: Scaling function of the transverse Kubo relaxation function T (r, (p,Tr) at 
tp = 1.294 versus r T and r = J{q D /q) 2 + (l/g£) 2 . 
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Figure 3.10: Scaling function of the longitudinal Kubo relaxation function (f) L (r,(p,r L ) at 
tp = 1.294 versus t l and r = J{q D /q) 2 + (l/g£) 2 . 
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Figure 3.11: Scaling function of the transverse Kubo relaxation function (f) T (r,(p,T T ) at 
(p = 0.97 versus Tt and r = J (q D /q) 2 + (l/g£) 2 . 
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Figure 3.12: Scaling function of the longitudinal Kubo relaxation function (p L (r,ip,r L ) at 
= 0.97 versus r L and r = J (q D /q) 2 + {1/qC,) 2 - 
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IV. APPLICATION TO EXPERIMENTS 



In this chapter we review the results from the numerical solution of the mode coupling 
equations for dipolar ferromagnets above the transition temperature and compare them with 
experiments. Furthermore we give a short account of the experimental data recorded below 
T c and future theoretical developments. 

There are several experimental techniques like neutron scattering (NS), electron spin 
resonance (ESR), magnetic relaxation (MR), hyperfme interaction (HFI) and muon spin 
resonance (/xSR) with various complementary characteristics: different wave vector range, 
short-range point-like probes in real space (HFI, /iSR) vs. probes in reciprocal space (NS, 
ESR, MR). The critical dynamics of isotropic ferromagnets such as EuO, EuS, Fe, Ni and 
many other materials have been studied by one or several of the above methods. 

To emphasize the decisive role of the dipolar interaction we note that the experimental 
situation prior to the theoretical work in Refs. [79-81] was, however, puzzling in many ways. 
On the one hand, in hyperfme interaction experiments on Fe and Ni one observed a crossover 
in the dynamical critical exponent from z = 5/2 to z = 2 [228 , ^5 , 113 , 11^ 



i.e., a crossover 



to a dynamics with a non conserved order parameter. This was confirmed by electron spin 



resonance and magnetic relaxation experiments ||147| - |150|J6^ , |152|| , where a non vanishing 
Onsager coefficient at zero wave vector was found. These results indicated that the critical 
dynamics of isotropic ferromagnets cannot be explained solely on the basis of the short range 
exchange interaction, which would lead to an exponent z = 5/2 for the whole wave vector 
range. However, on the other hand a critical exponent z = 5/2 was deduced from the wave 
vector dependence of the line width observed in neutron scattering experiments right at the 
critical temperature by three different experimental groups in Refs. [I~8|,[5"7l, [|188H192| and 
26 l| , |28| -p0|] . Nevertheless and even more puzzling, the data for the line widths above the 



transition temperature [[188|J189| , |31| could not be described by the Resibois-Piette scaling 
function resulting from a mode coupling (MC) theory ||231|| and a renormalization group 
(RG) theory ||1 74| , |1 23f | , which take into account the short range exchange interaction only 
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(see chapter 0) . 

Those apparent experimental discrepancies could be resolved by the mode coupling the- 
ory |f75-8*I| described in chapter [TIT], which on top of the exchange interaction takes into 



account the dipole-dipole interaction present in all real ferromagnets. The results of this 
theory are reviewed in the following and compared with the experimental data. 

A. Neutron scattering 

Inelastic neutron scattering has been used since decades to investigate the spin dy- 
namics of magnetic systems providing information about the spin correlation function 
48| , |194|208| - |210| . [57|| . It has become one of the most important experimental techniques 



for studying the properties of condensed matter systems. 

The cross section for magnetic inelastic scattering of polarized neutrons (with polariza- 
tion P) is given by [ |18iyi66|| 



1 r°° 

x ■ 



die** (S^ ■ S^,(t)> + zP ■ (S^ x S^,(t)> , (4.1) 

-oo 



2%h 

where the spin fluctuations transverse to the neutron wave vector transfer (scattering vector) 
Q = k — k' are defined by S/~ d = Q x [S/^ x Q] /Q 2 . = 1 + d denotes the position vector 
of a nucleus in a rigid lattice, where d is the site of the nucleus (with a gyromagnetic ratio 
gd and a form factor -Fd(Q) within the unit cell 1). The quantity r = —0.5391 fm is a useful 
unit for the magnetic scattering length. 

In neutron scattering experiments using unpolarized neutrons one measures the trans- 
verse scattering function S T (Q,u) 



S T (Q, oc ( S°" - 9^0- ] (S^)S^(-co)) , (4.2) 



cap _ V V f.Q<*( f .\qP 



where q denotes the neutron wave vector transfer with respect to the nearest reciprocal 
lattice vector r, i.e. the scattering vector reads Q = q + r. Eq. (|4.2| ) implies that only spin 
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fluctuations transverse to the neutron wave vector transfer q can be measured in the forward 
direction, where r = 0. In order to measure the longitudinal fluctuations one has to go to a 
finite Bragg peak r with q perpendicular to Q ~ r. The separation of the longitudinal and 
transverse spin fluctuations can be achieved by using a polarized neutron beam. A typical 
scattering geometry is shown in Fig. 4.1, where the polarization P is defined by a small 
guide field H v . For more details on neutron scattering using polarized neutrons we refer the 
reader to the books by Marshall and Lovesey ||181|| , Williams [[26811 and Lovesey ||166 |. 



Recent advances in high resolution neutron scattering techniques, like e.g. the Neutron 
Spin Echo (NSE) technique ||187|| , have allowed to extend the early studies to critical fluctu- 



ations of substantially larger wavelengths and even made it possible to distinguish between 
different line shapes ||190 . 



In this subsection we compare the results of mode coupling theory with neutron scat- 
tering data on EuO, EuS and Fe. For a quantitative comparison we need several material 
parameters (non-universal quantities) like the magnitude of the Brillouin zone boundary, 
exchange constant, etc. (see Table |V111[ ). 

1. Shape crossover 

Due to the recent advances in high resolution neutron scattering techniques it became 
possible to resolve not only the width but also the shape of the spin correlation func- 
tions. This progress made it possible to test the predictions for the line shape from MC- 
theory ||262| , p.l8|| and RG-theory f6l|]22|| , which take into account solely the short-range 



exchange interaction. The latter theories for the critical dynamics were expected to fit the 
experimental data in an increasing quantitative way as one was able to measure the correla- 
tion functions at smaller and smaller wave vectors, since the influence of additional irrelevant 
interactions should diminish as one moves closer to criticality. It came as quite a surprise 
when Mezei found a nearly exponential decay of the Kubo relaxation function by spin echo 
experiments on EuO at q = 0.024A -1 and T = T c [|190||, which was in drastic disagreement 
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with the bell-like shape predicted |262| , |TT8i |6Tl,p2| . 



Actually - in contrast to the expectation - only at larger wave vectors the RG and 
MC-calculations, taking into account exchange interaction only, were in accord with the 
experiment. This was shown by a comparison [77[] of the shape of the RG-theory |6l| . |22| 



with the heuristic shape of Ref . |77| . Also, a discrete version of the MC-theory gives quite 
reasonable an agreement with experiment for large wave vectors close to the Brillouin zone 
boundary []53| , ^4| . 

In order to explain the anomalous exponential decay found by Mezei [ 190|1 , it was as- 



serted in Refs. ||167| , |i~5f that neither MC- nor RG-theory is valid in the time regime probed 



by Mezei's ||190|| experiments. To fit the experimental results a "hybrid theory" [H| was pro- 



posed, which is a phenomenological interpolation scheme between the short and long time 



limits. In contrast it was shown on the basis of a mode coupling theory |82| , |85| , which take 
into account the effects of the long-range dipolar interaction between the spins, that the 
experimental data can be explained quite naturally without any need of a special treatment 
of the short-time behavior. The results of this mode coupling analysis will be reviewed in 
this subsection. 

We have seen in chapter |TJ that the dynamic crossover in the line width sets in at 
wave vectors one order of magnitude smaller than found in the static susceptibility. Here 
we consider the shape crossover, i.e., changes in the functional form of the spin correlation 
function. From the solution of the dipolar MC-equations one finds the following scenario 
for the shape crossover. 

• In the immediate vicinity of the critical temperature the transverse relaxation function 
shows a nearly exponential decay in time for wave vectors smaller than the dipolar 
wave vector q D . Passing to wave vectors larger than g D there is a shape crossover to a 
Hubbard- Wegner shape [ |118|J262|| corresponding to the isotropic Heisenberg fixed point 
(see chapter [D]). This shape crossover takes place close to the dipolar wave vector q u 
in remarkable contrast to the crossover in the line width which sets in at wave vectors 
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almost one order of magnitude smaller. 

• For wave vectors smaller than q D and for temperatures in the intermediate vicinity of 
the critical temperature the longitudinal Kubo relaxation function exhibits a Gaussian 
decay at small times and damped oscillations for larger times. This shape is quite 
different from the exponential decay found for the transverse scaling function. Passing 
to larger wave vectors the shapes of the transverse and longitudinal relaxation functions 
become identical (Hubbard- Wegner shape). 

• For temperatures well separated from the critical temperature (i.e. for temperatures 
where the correlation length becomes less than the typical dipolar length scale q~ x : 
g D £ < 1) the shape of the transverse and longitudinal relaxation becomes identi- 
cal. Upon passing from small to large wave vectors both relaxation functions show 
a crossover from a hydrodynamic shape, which is nearly exponential, to the bell like 
Hubbard- Wegner shape. 

In order to be specific and in view of spin echo experiments on EuO above T c [ |192|jl93| | 
we have displayed the Kubo relaxation function versus time for a set of wave vectors (q = 
0.018A" 1 , 0.025A" 1 , 0.036A" 1 , 0.071A" 1 , 0.150A- 1 ) in Figs. 4.2-5. In Figs. 4.2 and 4.3 
these are plotted versus the scaled time variables r a for a temperature in the immediate 
vicinity of T c (T = T c + 0.25i^) and in Figs. 4.4 and 4.5 for a temperature T = T c + 8K well 
separated from T c . For wave vectors q <C q u = 0.147 A -1 one can infer from Fig. 4.2 that the 
transverse scaling function decays nearly exponentially in time. This implies a Lorentzian 
like shape for the frequency dependent relaxation function for q <C q D (solid line in Fig. 4.2). 
For larger wave vectors, q > q D , the curves look similar to Gaussians for small times and 
oscillate for larger times, > 3 (see point-dashed curve in Fig. 4.2). The oscillations of 
the longitudinal relaxation function (Fig. 4.3) leads to the side peaks in Fig. 4.9, where the 
longitudinal correlation function is plotted versus frequency for a fixed wave vector. 

Further away from T c the line shape crossover of the transverse relaxation function is 
much less pronounced and the shape resembles more a Gaussian even for wave vectors 
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much smaller than q u . In Figs. 4.4 and 4.5 the transverse and longitudinal Kubo functions 
are shown for the same wave vectors as in Figs. 4.2 and 4.3 at T = T c + 8K. For this 
temperature the shape of the longitudinal and transverse Kubo function is nearly the same 
(compare Figs. 4.4 and 4.5). This is precisely what one would have expected, because the 
influence of the dipolar forces decreases with separation from the critical point and then 
there is no difference any more between longitudinal and transverse modes. 

To single out the line shape crossover the appropriate time scales are r a of Eq. |3.39 



used in Figs. 4.2-5. On the other hand, for comparison with experiments it may be 
more convenient to present the transverse Kubo relaxation function versus the time scale 
t' = Ai or g 5 / 2 t. Such plots are exhibited in Fig. 4.6 for T = T c and the wave vec- 
tors q = 0.018A- 1 , 0.036A- 1 , 0.150A- 1 , 0.3A- 1 and in Fig. 4.7 for T = T c + 0.5K and 
q = 0.018A" 1 , 0.025A" 1 , 0.036A" 1 , 0.071A" 1 . 

At T = T c triple-axis spectrometer scans at q > 0.15A -1 ||193|| show non-Lorentzian line 



shapes in agreement with predictions based on models with short range exchange interaction 
only and in agreement with the above MC result. Neutron Spin Echo (NSE) studies at much 
smaller wave vectors |190| , |193| | lead to Lorentzian line shapes in agreement with the above 
crossover scenario for dipolar ferromagnets. In Fig. 4.8 data from NSE studies [19U| on EuO 



right at the critical temperature are shown for the transverse Kubo relaxation function at the 
wave vector q = 0.024A -1 (solid line) versus time in nsec. We have used the theoretical value 
for the non universal scale A = 7. 1/5. 1326 me V A 5 / 2 (see Table |VTTTD . There is an excellent 
agreement with the experimental data for t < 1 nsec. The experimental data are above the 
theoretical curve for t > 1 nsec. This may be due to finite collimation effects in this time 



domain, as noted by Mezei ||190|| . To substantiate this point we have also plotted in Fig. 4.8 
the relaxation function at q — 0.028A -1 (point-dashed curve), which is significantly higher 
than the curve for q = 0.024A -1 for t > 1 nsec. The fairly large difference of the curves with 
q = 0.024 A -1 and q = 0.028 A -1 is implied by the close proximity to the crossover region. 
In order to exhibit the difference from the MC theory including only short range exchange 



interaction the dashed curve in Fig. 4.8 shows the solution of Eqs. [3.40| and |3.41| for this 
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special case, i.e., y = 0, x = and p cut = q BZ /q with q = 0.024 A -1 . The result differs 
drastically from the actual line shape including the long-range dipolar interaction. It is 
important to realize, that the crossover in the line shape starts nearly at q D , whereas the 
line width still scales with the isotropic critical dynamic exponent z = 5/2 in this wave 
vector region. 

Very recently the theoretical predictions concerning the line shape of the longitudinal 
relaxation function [S2],j35|l have been verified experimentally |33],|93|. Fig. 4.9 shows the longi- 
tudinal Kubo relaxation function for EuS at T = 1.006T C and q = 0.19A -1 . The normalized 
experimental data de-convoluted by a maximum entropy method agree quite well with the 
theoretical prediction (solid line). Especially the double peak structure corresponding to the 
over-damped oscillations in time are observed. 



2. Line width crossover 



From the frequency and wave vector dependent Kubo relaxation function one can define 
a characteristic line width by the half width at half maximum (HWHM) 

= iite = 0)] , (4.3) 

which obeys the scaling law 

rH W (?)= A HW^7H W (^)- (4.4) 

Let us start by comparing the theoretical and experimental line widths for Fe, EuO and 
EuS precisely at the critical temperature. Fig. 4.10 shows the scaling functions for the 
transverse and longitudinal line width in the Lorentzian approximation, normalized such 
that the line width approaches 1 for q ^> q D . As already noted earlier, a remarkable result 
of mode coupling theory is that the dynamic dipolar crossover of the longitudinal line width 
starts near the dipolar wave vector q ~ q D , whereas for the transverse line width it is shifted 
to a wave vector which is about one order of magnitude smaller than q D . This explains 



Re 
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why no appreciable deviation from the exchange scaling prediction T T (g) oc q 5 / 2 in the wave 
vector range accessible up to now was found in Refs. |4^j5^], ||1 88| |1 92|| and [[264| ; ^8H 30| in 
neutron scattering experiments right at the critical temperature. The experimental data in 
Fig. 4.10 are collected from Refs. H|57||, p8| - M and |264|,|28|-g0|] and are normalized with 
respect to the theoretical value of the non-universal frequency scale Ai or = 5.1326 x A. For a 
comparison the solid line in Fig. 4.10 shows the scaling function of the transverse line width 



as it has been obtained from the solution of the complete mode coupling equations |82|,|85 



Whereas the general functional form is quite similar to the scaling function obtained in 
the Lorentzian approximation, their asymptotic values at large values differ by a factor of 
approximately 1.2, which leads to an even better quantitative agreement between theory 
and experiment. It still needs experiments at smaller momentum transfer to detect the 
increasing of T because of the crossover to the dipolar regime. 

In contrast to the data on the wave vector dependence right at the critical temperature, 
experiments on the scaling function above the critical temperature revealed huge deviations 
from the Resibois-Piette scaling function, as is evident from the neutron scattering data on 



Fe ||188|| shown in the scaling plot in Fig. 4.11, which become the more pronounced the 
closer to T c . 

This figure also shows a comparison of the HWHM (dashed-dotted line) of the transverse 
Kubo functions resulting from Eqs. |3.40| and |3.41| and the Lorentzian line width 7^. (solid 
line) determined by Eq. |3.49| . In the case of Fe curves a, b, c and d correspond to the 
temperatures T — T c = 1AK, 5.8K, 21.0K, and51.0i^. All curves are normalized with 
respect to their value at criticality. This implies that the non universal frequency scale is 
then found to be Ai or = 5.1326 A = 107.2 meV A 5 / 2 in the Lorentzian approximation and 
A HW = 1.37 x A lor = U7meVA 5 / 2 for the HWHM. The latter value is in quite reasonable 
agreement with the experimental value A cxp = 130 meV A 5 ^ 2 . The Lorentzian line width 
and the HWHM have nearly the same (g£) -1 -dependence. Especially for small wave vectors 
the difference comes about only because of the different non universal frequency scales. In 



summary, the line width is affected by the line shape crossover only to a minor extent, since 
the former is only a coarse feature of the relaxation function. At small wave vectors only 
the non-universal scale is modified. This allows one to use the Lorentzian approximation for 
some later applications of the theory to NS, ESR, MR, and HFI experiments. 

As can be inferred from Fig. 4.11 the complete solution of the mode coupling theory 
is in reasonable agreement with the experiment close to T c (T — T c = 1AK, 5.8K) and 
gives an improvement over the Lorentzian approximation. The minor differences may be 
due to the following reasons, (i) As mentioned above the measured scaling functions of the 
transverse relaxation function were fitted to an exponential line, which is not the correct 
shape, (ii) Because the dipolar crossover temperature T D of Fe is 8.6K, static crossover 
effects not taken into account in the experiments may cause some changes. Furthermore the 
non universal scale of the correlation length £ is affected by experimental uncertainties. A 
change in £ would lead to a horizontal shift of the data points in Fig. 4.11. 

Larger differences show up for temperatures further away from T c (T — T c = 21K, 51K), 
which cannot be accounted for by the shape crossover or static crossover effects. Here 
the measured line widths are larger than the theoretical. In order to explain this, it is 
necessary to take into account the van Hove terms and further relaxation mechanisms due to 
irrelevant interactions. Such irrelevant interactions are unimportant as concerns the critical 
behavior, but nevertheless may play an important role for temperatures well separated 



from T c (see Ref. [j35|j). Pseudo-dipolar forces have been studied in Refs. ]85| and Ref. [[|. 
Those additional interactions (crystal field interactions, the spin orbit interaction leading 
to pseudo-dipolar terms) are presumably less important in magnets with localized moments 
such as EuO and EuS. This can be seen from the comparison of the line width with the 



experimental data. Fig. 4.12 (Fig. 1 from Ref. ||193|| , and Fig. 1 from Ref. ||192|| ) shows a 
scaling plot of the transverse line width data for temperatures T — T c = 0.5K, 2K, and AK 
in EuO, which are in quantitative agreement with the mode coupling theory (solid lines). 

Recently the transverse and for the first time the longitudinal spin fluctuations have been 
measured on EuS using inelastic scattering of polarized neutrons p3jp3| . The observed re- 
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laxation rates follow rather precisely the transverse and longitudinal line width calculated 
from the dipolar mode coupling theory (see Figs. 4.13 and 4.14). Especially, these mea- 
surements confirm the behavior of the longitudinal line width predicted by mode coupling 
theory gggl . 



3. Constant energy scans 

In this subsection we return to an analysis of the line shape and review results obtained 
by scans at constant frequency as opposed to constant wave vector in the previous sections. 
As will become clear in the following, certain features of the line shape of the correlation 
function can be accentuated by constant energy scans for the scattering function S T (q,u). 
Characteristic quantities in these scans are the peak position q and the HWHM Aq. For 
an isotropic Heisenberg ferromagnet RG-theory predicts a flat curve for the reduced peak 
positions q (A/to) 2 ^ 5 plotted versus the scaled frequency to^ 5 ^ 2 /A ||123| , |78| . This theoretical 



result has been confirmed experimentally in certain energy and wave vector regions |29 
However, more recently Boni et al. |31[ have found pronounced departures from the isotropic 
scaling law in EuS in a region, where according to Refs. [^,^ll] the dipolar interaction should 
have a considerable effect on the dynamics. This is quite similar to the situation in constant 
wave vector scans, where one finds deviations from the Resibois-Piette scaling function 



188,189 



One of the most striking new features introduced by the dipolar interaction is the gen- 
eralized dynamical scaling Eq. ( |3.31| ) with the additional scaling variable y = q n /q. In 
conventional constant q scans plotted versus the scaling variable x = 1/q^ the line widths 
are not represented by a single scaling function but by a series of curves as exhibited in 
Figs. 3.1 and 3.2., where each curve corresponds to a fixed temperature. The failure of the 
isotropic scaling law in constant energy scans can also be attributed to the influence of the 



dipolar forces p2|,|85[] 



The generalized scaling law for the Kubo functions and characteristic frequencies leads 
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to the following scaling law for the peak position in constant energy scans |85 



, A .2/5 



UJ J 



{ \ 

UJ 

V> — s z/2 



(4.5) 



From the above generalized scaling laws it becomes obvious that one obtains a set of curves 
parameterized by the scaling variable <p> = arctan(<7 D £), if Eq. (|4.5|) containing the two 
parameter scaling function Q is plotted versus uj^ 5 ^ 2 /A. In Fig. 4.15 the scaled peak position 
g (Ai or /u;) 2//5 is plotted versus uj = u;£ 5 / 2 /A lor for the following set of scaling variables <p = 
a) 1.490, b) 1.294, c) 0.970 and d) <p = 0. Case d) corresponds to an isotropic ferromagnet, 
i.e., q D = 0. Ai or is related to A by A lor = 5.1326A. The above values for ip correspond 
in the case of EuS (q u = 0.27A -1 , £ = 1.8lA) to the temperatures T = a) 1.01T C , b) 
1.06T C , and c)1.21T c . Due to the generalized scaling law Eq. ( |4.5| ) the curves coincide with 



the isotropic theory for high frequency but deviate for small frequency. The frequency, 
where the deviation from the isotropic result sets in, increases upon approaching the critical 
temperature. 

The steep drop off of the scaled peak positions at particular ip dependent values of uj 
has been explained as follows | 85| . The dipolar forces imply that the order parameter is 



no more conserved and the uniform relaxation rate r o becomes finite. Hence the scattering 
function S T (q,u) remains finite for vanishing wave vector 

S T (q = 0,u;)^-^- for T >T C . (4.6) 

Since r o is proportional to Eq. ( |4.6|) reduces to 

S T (q = 0,uj) oc \ for T = T C . (4.7) 

Therefore, the maximum of the constant energy scan at q = increases strongly with 
decreasing frequency lu, whereas the local maximum at finite q is shifted to smaller q (see also 
the inset in Fig. 4.15). As a result of this competition only the maximum at q = survives 
for low frequencies. In order to substantiate this, typical constant energy scans are shown in 
the inset of Fig. 4.15, where S T (q, u)/S T (0, u) (in arbitrary units) is plotted versus \ jr for 
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if = 1.294 and for a set of scaled frequencies ( Cj = 10 L / 10 with L — 8, 10, 12, 14, 16 indicated 
in the graph). The corresponding scaled peak positions are indicated in the scaling plot by 
crosses. This behavior explains, why the scaled peak positions for dipolar ferromagnets show 
such a steep drop off at small frequencies. 

Finally, it is important to note that the characteristic deviations from the isotropic theory 
exhibited in Fig. 4.15 result from the crossover in the time scale and not so much from the 
crossover in the shape function. This can be inferred from Fig. 4.16 where the scaled peak 

/ .\ 2/5 

positions q ( — J in Lorentzian approximation is shown. Compared to Fig. 4. 15 the maxima 
are overemphasized in Fig. 4.16, but the portions with the steep slope are nearly identical. 
The differences of the exact mode coupling theory and the Lorentzian approximation can 
also be seen in the insets of Figs. 4.15 and 4.16. 



B. Electron spin resonance and magnetic relaxation 

In electron spin resonance (ESR) and magnetic relaxation (MR) experiments one mea- 
sures the electronic response function at zero wave vector and finite frequency 

r"(q,^) X a (q) , 

X(q=0,u;) = - — —— Hq-tf, 4.8 

iu + r Q (q, cu) 

where r a (q, u) is the frequency dependent relaxation function in Eq. ( |3.19| ). Therefrom one 



determines the kinetic coefficient 

L( W ) = r a (q, W )f(q) | q ^o (4.9) 

for the homogeneous magnetization dynamics. In Lorentzian approximation one finds for 
the kinetic coefficient at zero frequency 

T( c\ g2j2 ' u rr 1 jui.2 X L (k,q D )x T (k,q D ) 

Lo(%,t) = ^KTJ dkk TL{ktqu)+r r {ktqjs y ^ 

where the dependence of the Onsager coefficient on the dipolar wave vector and the tem- 
perature is now indicated explicitly. Note that the kinetic coefficient for the homogeneous 
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relaxation is the same for the longitudinal and transverse spin fluctuations. Using the gen- 
eralized dynamic scaling and introducing polar coordinates (r = a/x 2 + y 2 , y/x = tan</?) 
one finds |nj] 

L (q D ,O = BF(l/q D O, (4.11) 
with the universal crossover function 

Wfcfl = fl + 4.V 7/4 Lrr^ pl^fff* , (4.12) 
V J 7 L (^)+7 T (^) 

and the non-universal constant 

B = ^Aq^. (4.13) 

If there were no dipolar interaction [g = 0), one would simple find a vanishing Onsager 
coefficient due to the factor g 2 in equation ( [4.10|) . With regard to the crossover function F 
of Eq. (|4.12|) it is natural to define the reduced crossover temperature by 

r D = ^V^ = (& 2 ) V0 , (4-14) 

c 

where the crossover exponent equals the susceptibility exponent 7 for a ferromagnet with- 
out dipolar interaction ||. If one neglects the dipolar crossover of the correlation length, 
the scaling variable g 2 £ 2 can be written as g 2 £ 2 = (r cross /r) 7 and the crossover temperature 
in terms of the dipolar wave vector is given by 



The crossover temperatures T D resulting from Eq. ( |4.15| ) can be found in Table [IX . 



The Onsager coefficient L = T a x a does not depend on the sample shape ]7I] and is the 
same for the transverse and the longitudinal mode. The universal crossover function F(p) 
with p = l/q D £ is plotted in Fig. 4.17 in units of its value at criticality F(0). For tempera- 
tures larger than the dipolar crossover temperature F(p) shows a £ 7//2 power law behavior, 



which was first shown in Refs. |p.22j , p27 [. In the strong dipolar limit, i.e., very close to the 
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transition temperature, the Onsager coefficient approaches a constant |72| reflecting the non 
conserved nature of the order parameter due to the presence of the dipolar interaction. 

The kinetic Onsager coefficient at zero wave vector has been measured by magnetic re- 
laxation experiments ||147| . |9lB| . |i~4l||64T| and ESR experiments ||148| , |147| , |149| , |150|| for the Onsager 
coefficient in EuS, EuO, CdCr 2 S 4: , CdCr^Se^ and Ni | |236| , |238|| . In order to extract the 
contribution of the critical fluctuations one has to subtract a non-critical background Lb g 
which in the critical region of all ferromagnets (except Ni) is very small, L^ g <C L cr (see 



Table |LX|) . The critical part of the Onsager coefficient 



hL (q D ,O=L d ^ = BF(p) 



(4.16) 



is compared with the crossover function F in Fig. 4.18 (Fig. 1 of Ref. [152 1). The non- 
universal parameter hL& = BF(0) was fixed by the data for the kinetic coefficient in the 
center of the critical region. There is an excellent agreement (except Ni and EuO) of the 
observed crossover with the results of the mode coupling theory. Especially, the data for EuS 
follow the theoretical function quite closely. The ESR-data on Ni [|236|[, shown in Fig. 4.18, 



deviate from the mode coupling result. Sporel and Biller ||23£ ], however, find an increase in 
the electron paramagnetic resonance line width of Ni near T c , in accord with mode coupling 
theory but opposite to what was found in Ref. ||236|1 . It seems that the observation of the 



peaklike broadening depends critically on the quality of the sample. Recently Li et al. [165 



have demonstrated that it is possible to observe the critical broadening of the linewidth in 
ultrathin Ni films. As further possible sources for the deviations in Ni and EuO from the 
dipolar crossover sample imperfections (e.g. oxygen vacancies in EuO, and internal stress 
in Ni) have been suggested ||152|| . 

The values of the fit-parameter L<j are listed in Table [IX. But, even these values agree 
with the theoretical prediction as can be inferred from the inset in Fig. 4.18. The minor 



differences can be attributed to some uncertainties in the dipolar wave vector ||152 |. 

Because ESR experiments are performed in a magnetic field B there should be an effect 
on the relaxation rate starting at 7I? = T(q = £~ l ) [140]. Those data points are not given 
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in Fig. 4.18. The influence of the magnetic field on the relaxation rate at zero wave vector 



has recently been studied by Kotzler et al. |fL54j| . In analyzing the effect of the magnetic 
field in terms of the internal isothermal susceptibility Xmt = 5M/5H- mt (H- mt is the internal 
magnetic field related to the external magnetic field by demagnetization corrections) the 
data were found to collapse on a single curve. Hence, empirically, the field effect on the 
kinetic coefficient can be accounted for by using q^Xmt instead of q^ 2 as the scaling variable 
and the same scaling function F as for the case of zero magnetic field. 

Up to now we have studied the behavior of the uniform relaxation at zero frequency only. 
Actually, the kinetic coefficient is frequency dependent, which was observed first by De Haas 



and Verstelle p8| . Measuring the uniform dynamic susceptibility in EuO they have found 
a deviation from a Debye spectrum (Lorentzian form), which according to Eqs. (|4.8|) and 
(|4.9|) is equivalent to a frequency dependent kinetic coefficient. Similar deviation from a 



Lorentzian have been found quite recently in EuS [0,[59|]. In order to study the frequency 
dependence of the kinetic coefficient one would have to solve the full frequency dependent 
mode coupling equations and therefrom deduce L(q D ,£,u). As a first approximation one 
can use |T55]|59l 

which is obtained as the first iteration step in a self consistent determination of the frequency 
dependent kinetic coefficient, starting with the damping coefficients T L,T (q, q D ,£) from the 
Lorentzian approximation. A scaling analysis of the frequency dependent kinetic coefficient 
gives 

L(q ,£,u) = BF(l/q D £M z /A) , (4-18) 

with the scaling function 

F(l/g D £ |W /A£) = (l + J-) ' / drr^ 



1 \ V4 7^5/2 X L (r,<p)x T (r,<p) 



^d^ 2 / n 7 L ( r ) V) + 7 T ( r ) v) + i{u/kq z D )r z sin 2 <p 

(4.19) 
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The real and imaginary part of F is shown in Figs. 4.19a and 4.19b versus the scaled 
frequency Cj = uj / Kq z D for several values of <p = arctan(g D £), indicated in the graphs. For 

small and large values of Cj the real part of the scaling function behaves asymptotically as 

(l-(Cj/Cj cl (l/q D 0) 2 forc2>«l , 
HaF(l/ fo e,w/A^) « F(l/q D 0{ 1 iy (4.20) 

{{Cj/Cj c2 {l/q D 0)- l - l/z foro>»l . 

The corresponding scaling functions for the Onsager coefficient at zero frequency F(l/q D £) 
and the scaling functions u C 2(l/q D £) and uj c \{l/q D £) characterizing the large and low fre- 
quency behavior, respectively, are shown in Figs. 4.20a,b,c. As can be inferred from these 
figures, the scaling function u) c i(l/q D £) is nearly constant. This finding is in accord with 
experiments by Dombrowski et al. |59[]. The experiments have been analyzed assuming a 
Lorentzian shape for the kinetic coefficient 



/fef • (4-21) 
1 + iuj/uj c 

This shape differs from the theoretical results presented above, especially in the large fre- 
quency limit, i.e. for uj ^> Aq z D . At low frequencies the theoretical result has the same 
expansion as the Lorentzian approximation used in analyzing the experiments. It is found 
experimentally that the characteristic frequency u c is nearly independent on temperature. 
This is in accord with the above theoretical result, that the scaling function u c i(l/q D £) is 
nearly constant over the whole temperature range. For a more detailed comparison with 
the experiment further analysis of the data on the basis of the results of the mode coupling 
theory for the frequency dependence of the Onsager coefficient seems to be necessary. It 
would be interesting to see, whether the experiment confirms the predicted large frequency 
behavior L(u) oc uj^ 1 ^ 1 ^. 

Recently, the homogeneous magnetization dynamics has also been investigated in the 
ferromagnetic phase p9]| . It is found that the scaling function for the kinetic coefficient 
below T c agrees exactly with that observed earlier above T c . Since a complete theory of 
the critical dynamics below T c is still lacking, these results are not explained yet. However, 
those experimental findings are a clear indication of the importance of the dipolar interaction 
below the critical temperature. 
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C. Hyperfine Interactions 



There are several nuclear, i.e., hyperfine interaction (HFI), methods commonly used to 
study critical fluctuations in magnets. These are nuclear magnetic resonance (NMR) ||235 



Mossbauer effect (ME) P^M , perturbed angular correlations (PAC) [p28| j95|j4^JlTl of 
gamma rays, and muon spin rotation (/xSR) (see also section [IV (J| ). The application of 
nuclear techniques to study critical phenomena in magnets has recently been reviewed by 



Hohenemser et al. [l.T3|. For recent work on the application of /i-SR for the investigation 
of spin dynamics the reader may consult Refs. | |239| , |5!^|HI3| . All of the hyperfine interaction 
methods are local probes which are related to a wave vector integral of the spin correla- 
tion functions. As such they offer a complement to neutron scattering. Dynamic studies, 
using hyperfine interaction probes, utilize the process of nuclear relaxation produced by 
time-dependent hyperfine interaction fields, which reflect fluctuations of the surrounding 
electronic magnetic moments. 

The hyperfine interaction stems from the magnetic interaction of the electrons with the 
magnetic field produced by the nucleus. The hyperfine interaction of a nucleus with spin 
I, g-factor and mass itin with one of the surrounding electrons with spin S and orbital 
momentum L can be written in the form |20|J2^ 



_ Ze 2 g N 

^hyp 



ll.L + ^)(x)I.S-ll.S + 3(I " x)(S - x) 



y3 ^ 



(4.22) 



The first term represents the interaction of the orbital momentum of the electron with the 
nuclear magnetic moment of the nucleus. The second term is the Fermi contact interaction 
and the last two terms represent the dipolar interaction 0. The Fermi contact interaction 



8 Actually, also the Fermi contact interaction comes from the magnetic dipole interaction. In the 
evaluation of the matrix elements of the (residual) dipolar interaction, represented by the last two 
terms, an infinitesimal sphere around the origin has to be excluded. The singular part of the dipolar 
interaction coming from the interior of the sphere is represented by the Fermi contact term. 
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is finite only for electrons having a finite probability density at the nucleus, i.e. bound 
s-electrons or itinerant electrons. 

For a substitutional nucleus the main contribution to the hyperfine field comes from 



the bound s-electrons ||91|| , which are polarized by the magnetic ions. The large (and often 
negative) hyperfine field for transition-element ions comes from the polarization of the core 
s electrons by the spin density of the unpaired 3c? electrons, which then contribute to the 



hyperfine field through the Fermi contact interaction [261] . The residual dipolar interac 



tion is negligible in this case. Note also that this term is zero for cubic symmetry. The 
Hamiltonian, Eq. ([4.22 ), can also be used for the analysis of spin resonance experiments 



with muons (/xSR). However, these do not have bound electrons and hence the Fermi con- 
tact term involves only conduction electrons and is of the same order of magnitude as the 
(residual) dipolar interaction |56|J185|| (see also section |IVD| ). For instance in Pd, there is 



a dipolar contribution by the Pd /-electrons which also polarize the s-electrons which in 
turn contribute via the Fermi contact term to the hyperfine field at the muon site. 

After the above general remarks, let us now return to those hyperfine interaction probes 
(ME, PAC, NMR), where the Fermi contact term gives the dominant contribution to the 
hyperfine field at the nucleus. Then the corresponding interaction Hamiltonian reduces to 

H{t) = A 

contact 

I • S(t) . (4.23) 

If the spin auto-correlation time, r c , is much shorter than the Larmor period, 1/Vj,, and the 
nuclear lifetime, t^, (motional narrowing regime) the nuclear relaxation rate r^ 1 is directly 
proportional to the (averaged) spin autocorrelation time r c 

r" 1 = 0*7-0 (4.24) 



with a hyperfine coupling constant C n f | 113 |. The nuclear relaxation rate equals the spin- 
spin relaxation rate r^ 1 = T 2 _1 in the case of NMR, is proportional to the relaxation induced 
excess line width Ar for ME, and for PAC tr is given by the relaxation time of the pertur- 
bation factor. 
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The (averaged) spin autocorrelation time r c is denned by a time-integral 

+00 

= 2 / dt \ E GQQ ( r - °, *) = g E / q ^(q, w = 0) . (4.25) 

—00 

over the spin autocorrelation function G aQ (r, £) = ~({S a (r,t), S a (0,0)}). Hence hyperfine 
interaction methods provide an integral property of the spin-spin correlation function. Upon 
using the fluctuation dissipation theorem (FDT), which in the special case uo = reduces to 

G-(q,u = 0) = 2k B Tp^ y (4.26) 

one finds for the autocorrelation time 

_ KT f * l^X°(q,0) 



The q-integration extends over the Brillouin zone (BZ), the volume of which is V q . 

Important information about the behavior of the auto-correlation time can be gained 
from a scaling analysis. Upon using the static and dynamic scaling laws Eq. (|4.27 ) can be 
written as 

rcoe^/^V "^'^! - (4-28) 

If there were no dipolar interaction, one could extract the temperature dependence from the 
integral in Eq. (|4.28|) with the result r c oc This expression can be used to define an 

effective dynamical exponent z e g(r), which depends on the correlation length by 

'T - T„ N w 



r c cc e^" 1 oc • (4.29) 

In the presence of dipolar forces one finds after introducing polar coordinates as in sec- 



tion IV B 



where z = 5/2 and the non universal constant H is given by 



H = {KT)2 / 3 (q D af 2 . (4.31) 
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The lower cutoff ro is 



r 



1 + 4r, . ( 4 -32) 



where q BZ is the boundary of the Brillouin zone. In the critical region it can be disregarded 
and replaced by r = 0, since q BZ ^> q D and the integrand in Eq. Q4.30j) is proportional 
to y/r for small r. For very small £ (outside the critical region) the cutoff reduces the 
autocorrelation time with respect to the critical value. 

One should note, that the dominant wave vectors contributing to the relaxation time 
r c in Eq. fl4.27p are close to the zone center. Quantitative estimates for the wave vector 
range probed by yU-SR experiments on Fe and Ni, which also apply for other hyperfine 
experiments, are given in Ref. [230|. 



The averaged autocorrelation time is a sum of two parts 

r c = ^(r L + 2r T ), (4.33) 
which we call longitudinal and transverse relaxation times 

r a = HI a { V ) , (4.34) 

where we have defined 

( 1 V 3/4 7 Y a (r (o) 

IM = [ 1 + ^e) J o drV ~ r rkf)' a = L,T - (4 ' 35) 

These two relaxation times are shown in Fig. 4.21 as functions of the scaling variable l/q D £, 
where we have neglected the lower cutoff r . Since the static longitudinal susceptibility 
does not diverge one finds that the longitudinal relaxation time is non-critical, whereas the 
transverse relaxation time diverges like tt oc £. This corresponds to an effective dynamical 
exponent z e ^ = 2. When leaving the dipolar critical region there is a crossover to the 
isotropic Heisenberg region, where both curves join and the relaxation times tl = tt = r c 
are characterized by another power law r c oc £ 3 / 2 corresponding to an effective dynamical 
exponent z e s = 5/2. 
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Let us now compare with hyperfine experiments on Fe and Ni. Early PAC ||228| , |95| and 



ME ||143| , |144|| were performed well in the dipolar region and an effective exponent z = 2 was 
found. The crossover in the dynamic exponent from z = 2.5 to z = 2 was first observed by 
Chow et al. p4"| , |113|| . The autocorrelation time r c is shown in Fig. 4.22 in units of the non 
universal constant H (see Table [LX| ) versus the scaling variable 1/<? D £. The data points in 
Fig. 4.22 are results of HFI experiments on Fe and Ni [[228jj9^ , f44] , |113|1 for the autocorrelation 
time r c (in units of their non universal frequency scale H). As before there is no fit parameter 
for the scaling variable 1/<? D £. So we conclude that our MC-theory accounts well for the 
experimental data demonstrating the universal crossover behavior from z e ff = 5/2 to z e R = 2 
as the critical temperature is approached. 

Finally we would like to note that ME [[15|,f|(J and PAC j|9| experiments on the dynamics 
of Gd show neither z = 2.5 nor z = 2.0 but 1.74 < z < 1.82 depending on whether the system 
is assumed to have Heisenberg or Ising static exponents. Recent /x-SR experiments give 
similarly depressed values of the critical dynamic exponent z |[259| , [260]| (see section [IV D| ) . In 
Ref. | |229| the /i-SR data [ |259| , |105| ] have been compared with the results from MC-theory for 
isotropic dipolar ferromagnets [^T],^BJ, and it is found that the available data are found to be 
in agreement with the MC-theory. Gd is supposed to have an uniaxial exchange anisotropy 
along the c-axis. Hence one expects that in the asymptotic region there is one critical mode 
along the c-axis with z c g = 2 and two uncritical modes with z e g- = perpendicular to the 



c-axis ||72|| . Since for polycrystalline probes r c is an average over the relaxation rates of 
these modes, it is possible that the experimental data in the crossover region yield effective 
exponents less than z = 2. Here more experimental work with single crystals and a more 
detailed mode coupling analysis beyond the scaling analysis in Ref. JF2| would be needed to 
clarify the situation. 
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D. Muon Spin Relaxation (^SR) 



In muon spin relaxation (/iSR) experiments one observes the muon precession in an 
applied magnetic field or in a local magnetic field inside the sample^]. For recent work on 
the application of /zSR to the investigation of spin dynamics we refer the interested reader 
to Refs. pg^Pl. 

In analogy with other implanted probesQ like those discussed in section [IV C] , the muon 
spin will relax through interactions with fluctuating magnetic fields in its surroundings. 
However, in contrast the local magnetic field at the interstitial site of the muon has compa- 
rable contributions from the Fermi contact field and the dipolar field. In rare earth materials 
the residual dipolar field is high, and often the dominant field ||135 , |103|, 136 1 ■ It has a reduced 



symmetry as compared to the Fermi contact coupling. Whereas the Fermi contact coupling 
is isotropic in space, the residual dipolar field has an anisotropy with respect to the wave 
vector. In the case of a hyperfine field dominated by the residual dipolar field this implies 
that the spin fluctuations transverse and longitudinal with respect to the wave vector con- 
tribute with different weights to the /iSR relaxation rate. Hence, one has a situation quite 
different from e.g. PAC and ME measurements where the interaction of the substitutional 
probe with the host atoms is mainly a contact hyperfine coupling. This feature of /i-SR offers 
the possibility to design experiments which allow to distinguish between spin fluctuations 
longitudinal and transverse to the wave vector q [ g7l| . [27^ : |229| . 



Among the basic ferromagnets muon relaxation has been measured in Fe [|108|1 and 
Ni |[205|| , which have cubic crystal structure, and in Gd [|259| , |260|1 which crystallizes in a 



hexagonal lattice. The data for the muon relaxation rate in Fe were interpreted by a 



dynamic exponent z = 2.0 [108] which would indicate that the data are taken in the dipolar 



9 Muon spin relaxation experiments are often performed in zero applied external field 
10 The positive muon ends up in interstitial sites in solids, and in metals the actual site is in most 
cases identical to that preferred by hydrogen. 
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region. Actually, the data cover both the dipolar and the isotropic region - the crossover 
temperature for iron is T D — T c = 8.6K (see Table |V111| ) - with fairly large where the 
data error bars in the isotropic region, so that a crossover from z = 2.5 to z — 2.0 is 



not excluded by the data f69]| . The observations in Ni [|205|| seemed to be conflicting with 



PAC MM measurements, if one assumes that in both cases the contact field is the dominant 



contribution to the local magnetic field at the muon site. However, one has to take into 
account that the interaction of the muon with the surrounding is more complicated than for 
the other hyperfine probes, especially as already pointed out by Yushankhai [|273|1 one has 
to take into account the classical dipolar interaction between the muon magnetic moment 



and the lattice ion magnetic moments. This has recently been done in Ref. ||271|J272|| , the 
results of which we will review next. 

A zero-field /iSR measurement consists of measuring the depolarization function 

p z {t) m, 



P z (t) = ^Tr [p m a z a z (t)] 



(4.36) 



where p m is the density operator of the magnet and a z the projection of the Pauli spin 
operator of the muon spin on the 2-axis. TrLA] stands for the trace of A over the muon 
and magnet quantum states. Since the magnetic fluctuations are sufficiently rapid, the 
depolarization they induce can be treated in the motional narrowing limit. One finds [ [271 



P z {t) = exp [-\ g t] 



(4.37) 



with the damping rate 

vrX> f d 3 q 



A 5 



V J (2tt) ; 



{[G(q)C(q)G(-q)r + [G(q)C(q)G(-q)] w } 



(4.38) 



where 



C a/3 (q) 



\u=0 



(4.39) 



is the symmetrized correlation function at zero frequency, and 
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V 



An 



iI{9lHb)' 



(4.40) 



7 M = 851.6 Mrad s -1 T _1 is the gyro-magnetic ratio of the muon. The tensor G characterizes 
the type of interactions by which the muon spin is coupled to the magnetic moments of the 
material considered. It is given by 



G af5 {q) = r^H5 aP + D af3 (q) 



(4.41) 



The first term describes the Fermi contact coupling, where r M is the number of nearest- 
neighbor magnetic ions to the muon localization site and H measures the strength of the 
hyperfine coupling (see Table |VIII|) . The simple form for the contact coupling is valid only 
if the muon site is a center of symmetry. The second term characterizes the dipolar coupling 



and can be calculated via Ewald's method |271|| . To lowest order in the wave vector q one 
finds 



47T 



B a(3 {q = 0) 



q a qP 



(4.42) 



The values of the tensor components 5 Q/3 (q = 0) depend on the muon site. For muons in a 
tetrahedral or octahedral interstitial site in an fee lattice one finds B a/3 (q = 0) = |<5 a/3 [|156 
The situation is more complex for a bec lattice ||156 



For simplicity we confirm ourselves to the case of an fee lattice and refer the reader to 
Refs. ||271| , |272| , |229|| for a discussion of the more complicated cases of an bee cubic crystal or 
a hexagonal lattice structure. For an fee lattice the damping rate is found to be [[271|| 



A 3 



fio 2 k B T riBz 



3 4tt 



7, 



dqq 2 



2p 



x T {q) 
r T (g) 



2 V +{1 _ p yX L (q) 



(4.43) 



where we have defined 



P 



1 r.H 
- + J! — 
3 4vr 



(4.44) 



When the Fermi contact interaction is large compared to the classical dipolar interaction, 
i.e., \p\ ^> 1, Eq. |4.43j reduces to the result obtained in section |1V (J| 
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A' 



1 Mo,.. _ m2 fc s T /" 9SZ 



g X T (g) , x L (q) 

. ^(g) + r%) 



(4.45) 



67r 2 4vr v,/ ' /f 

On the other hand when the Fermi interaction is negligible, i.e. p ~ 1/3, we get from 
Eq. |Q3 

16^0 2 fc B T 



^dipola 



274vr 



7, 



x T (q) 2 x L (q) 



(4.46) 



While Eq. (|4.45| ) has been derived under the hypothesis that the muon spin interacts with 
the lattice spins through the isotropic hyperfine interaction, Eq. (|4.46|) has been obtained 
supposing that the coupling is only due to the classical dipolar interaction. One should 
note that the relative weight of the transverse and the longitudinal modes is reversed in 
going from a pure contact to a pure dipolar interaction of the muon with the lattice spins. 
Hence in combing PAC with /z-SR measurements it should be possible to distinguish between 
transverse and longitudinal relaxation rates. 

Upon introducing polar coordinates as in section |IV C| one can write the damping rate 
in scaling form 



A, 



2p 2 I T {v) + {l-p) 2 h{v) 



(4.47) 



where the non universal constant W is given by 

W 



^ 3/2 rwrlK' 2 



k B T c 



(4.48) 



3P V 47r g L jj, B 

The theoretical values of W for the four magnets considered primarily in this section are 
listed in Table [VIII| . The scaling functions I L,T (ip) are defined in section |IV C| . Fig. 4.21 



indicates that while It{^p) exhibits a strong temperature dependence in the whole temper- 
ature range of the critical region, Il(^>) is practically temperature independent for g D £ > 1. 
Therefore, the temperature dependence of the /xSR damping rate X z depends on the rela- 
tive weight of It{v) an d Iii'-f) which is controlled by the parameter p. As first noticed by 
Yushankhai ||273|| , the transverse fluctuations do not contribute to A 2 if r^H/Air = —1/3. In 
this case X z becomes temperature independent near T c . 
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In Figs. 4.23 and 4.24 experimental data for Ni [ |205|| and Fe ||108|| are compared with 
the results from MC-theory. One notices that most of the data have been recorded in the 
crossover temperature region where X z is not predominantly due to the transverse fluctu- 
ations. Again the MC-theory gives a quantitative explanation of the observed relaxation 
rates. 
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TABLE VIII. Material parameters for Ni, Fe, EuO and EuS. The parameter are defined in 
the main text. 





Ni 


EuO 




Fe 


a(A) 


3.52® 


5.12( 2 ) 


5.95® 


2.87® 


T c (K) 


627® 


69.1® 


16.6® 


1043® 




0.013® 


0.147® 4 ) 


0.245® 35 ) 


0.045® (0.033)® 


£o (A) 


1.23® 


1.57® 


1.81® 


0.95® (0.82) ® 


r^H/Air (T = K) 


-0.11® 


0(?) 


0(7) 


-0.51® 


T(A) 


0.00184 


0.725 


3.16 


0.0239 (0.0129) 


A ex p (meVfrl 2 ) 


350® 


8.7®(8.3)(® 


2.1®)(2.25)( 12 ) 


130(® 


A lor (meFl 5 / 2 ) 


241 


7.09 


2.08 


90.0 (123) 


W th (MHz) 


0.358 


5.35 (5.60) 


6.87 (6.41) 


2.98 (2.56) 


The values have been taken from the references (1) [|142fl 


, (2) gg, (3) m, (4) h, (5) g, 


(6) |p0l, (7) HH, 


(8) p|, (9) H 


I (io) @, (11) I 


3, (12) H, and (13) 


p8f. When no 


reference is indicated, the parameter 


has been computed from a formula given 


in the main text 



with parameters given in tables III. 1 and IV. 1. 
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TABLE IX. The dipolar wave vector q D , the Curie temperature T c , the experimental non 
universal constant for the Onsager coefficient Lj, L^ g and for the autocorrelation time H exp . 





EuS 


EuO 


CdCr 2 Si 


CdCr2Se4, 


Fe 


Ni 




0.245 


0.147 


0.058 


0.034 


0.045(0.033) 


0.013 


T C [K] 


16.6 


69.1 


84.4 


127.8 


1043 


627 




38 


24 


5.9 


4.4 




3.0 


hLb g (/j,eV) 


1.8 


0.64 


0.01 


0.01 






H cxp [W- 13 sec] 










6.0 


9.2 



The values are taken from Refs. [152] and the references in Table VII [ 
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Captions to the figures: 



Figure 4.1: Typical scattering geometry in neutron scattering experiments. The magnetic 
guide field Hy is used to define the beam polarization P. In order to measure the longi- 
tudinal spin fluctuations S L one has to go to a finite Bragg peak r with the wave vector 
q perpendicular to the scattering vector Q = r + q. For unpolarized neutrons only the 
transverse fluctuations can be measured in the forward direction, r = 0. 

Figure 4.2: Scaling function of the transverse Kubo relaxation function versus the scaled 
time variable r T at T = T c + 0.25K for q = 0.018A" 1 (solid), q = 0.0251" 1 (dotted), 
q = 0.036A- 1 (dashed), q = 0.071 A" 1 (long dashed) and q = 0.150A- 1 (dot-dashed). 

Figure 4.3: The same as in Fig. 4.2 for the longitudinal Kubo relaxation function. 

Figure 4.4: Scaling function of the transverse Kubo relaxation function versus the scaled 
time variable tt at T = T c + 8.0K for the same set of wave vectors q as in Figs. 4.2 and 4.3. 



Figure 4.5: The same as in Fig. 4.4 for the longitudinal Kubo relaxation function. 

Figure 4.6: Scaling function of the transverse Kubo relaxation function versus t' = A\ or q z t 
at T = T c for q = 0.025A" 1 (solid), q = 0.036A- 1 (dotted), q = 0.150A" 1 (dashed) and 
q = 0.300A- 1 (long dashed). 

Figure 4.7: Scaling function of the transverse Kubo relaxation function versus r' = K\ OT q z t 
at T = T c + 0.5K for q = 0.018A- 1 (solid), q = 0.025A- 1 (dotted), q = 0.036A" 1 (dashed) 
and q = 0.071 A" 1 (long dashed). 

Figure 4.8: Transverse Kubo relaxation function Q T (q,g,t) at q = 0.024A~ 1 (solid line) 
and q = 0.028A -1 (dashed-dotted line) for dipolar ferromagnets versus time in nsec. The 
dashed line is the Kubo relaxation function for short range exchange interaction only at 
q = 0.024A" 1 . Data points for q = 0.024A" 1 from Fig. 1 of Ref. £90 . 
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Figure 4.9: Normalized experimental data deconvoluted by a maximum entropy method. 
The solid line shows the spectral shape function predicted by mode coupling theory. The 
full circles (dashed line) indicate(s) the deconvoluted data. Figure taken from Ref. p3|. 



Figure 4.10: Scaling functions versus y^ 1 = q/q D at T c (in units of the theoretical non 
universal constant A) for (i) the HWHM of the complete solution of the MC equations 
for the transverse Kubo relaxation function (solid), (ii) transverse (point-dashed) and (hi) 
longitudinal line (dashed) width in Lorentzian approximation. Experimental data of the 
transverse line width for EuO ( U Ref. p§, Ref. |89| ) and Fe ( A Ref. |188| ). 



Figure 4.11: Scaling function of the transverse line width in Lorentzian approximation 
(solid lines) and of the HWHM of the complete solution (dashed lines) versus the scaling 
variable x = (<?£) 1 for a set of temperatures (T — T c = a) 1.4K, b) 5.8K, c) 21K, d) 51K). 
The Resibois-Piette function (solid) and the HWHM from the complete solution of the MC 
equations (dashed) without dipolar interaction is also plotted e). Experimental results for 
Fe from Refs. MM ( T — T c = (U) 1AK, (0) 5.8K, (A) 21.0K, (o) 51.0K ). 

Figure 4.12: Scaling plot of the temperature dependence of the relaxation rate of critical 
fluctuations in EuO above the Curie temperature T c = 69. 3K. The solid lines represent the 
results from mode coupling theory [80]. The data were taken between q = 0.018A -1 and 



q = 0.15A- 1 . Taken from Ref. p2|p3|1 (kj = 1/0 



Figure 4.13: Comparison of the longitudinal linewidth of Lorentzian fits normalized to 
A exp g 5 / 2 with A exp = 2.lmeVA 5 / 2 to the dipolar dynamic scaling function predicted by mode 
coupling theory for q D = 0.245A^ 1 for EuS. Figure taken from Ref. [|33|] (k = !/£)■ 



Figure 4.14: The same as Fig. 4.13 for the transverse widths. Figure taken from Ref. [33 



Figure 4.15: Scaled peak positions go(Ai or /^) 2 ^ 5 for constant energy scans of the scattering 
function for the complete solution of the MC equations versus the scaling variable cD for <p = 

110 



a) 1.490, b) 1.294, c) 0.970 and d) for the isotropic case, i.e., (p = 0. Inset: S T {q, u)/S T {0, u) 
in arbitrary units versus - for if = 1.294 for some typical values of the scaled frequency 

(u = 10 L / 10 with L = 8 (solid), 10 (dashed), 12 (point dashed), 14 ( ) and 16 

indicated in the graph by crosses). 

Figure 4.16: The same as in Fig. 4.15 for the scattering function resulting from the 
Lorentzian approximation. 

Figure 4.17: Universal crossover function F(l/q D £)/F(0) with F(0) = 0.1956 for the 
Onsager kinetic coefficient at zero frequency and zero wave vector versus the scaling variable 

Figure 4.18: Critical part of the Onsager kinetic coefficient for the homogeneous spin 
dynamics above T c of CdCr 2 Se 4 (Ref. p9fl ), CdCr 2 S 4 (Ref. |5(J), EuO (Ref. gig ), 
EuS (filled circles: Ref. open circles: Ref. M ), and M (Ref. |23§). The solid 

lines represent the result from mode coupling theory |^TJ. Inset: Nonuniversal ampli- 



tudes hLdq~ 3 ^ 2 for the kinetic coefficient and the relaxation rate at the critical point 
hL q q~ 1 q~ 1 / 2 = T T (q,T = T c )/q 5 ^ 2 fa 5.1326A for ferromagnets including Fe and Co. Figure 
taken from Ref. ||152|| . 



Figure 4.19: The real and imaginary parts of F are shown in a) and b) versus the 
scaled frequency Co = ui / Kq z D for several values of = arctan(g D £): cp\ = 0.99| (top curve), 
¥2 = 0.90| , ^3 = |(1 - 10- - 5 ), y 4 = f (1 - 10-°- 25 ), f 5 = f (1 - 10- - 1 ), = f (1 - 10-°- 025 ) 
(bottom curve). 

Figure 4.20: The scaling functions for the a) Onsager coefficient at zero frequency 
F(l/q D £) and b) the scaling functions uj C 2{l/q D ^) and c) io c \(l/q D £) characterizing the large 
and low frequency behaviour, respectively. 



Figure 4.21: Scaling functions Il,t{<p) f° r the transverse and longitudinal auto-correlation- 
time tl,t versus l/q D £- 
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Figure 4.22: Auto correlation time r c /H (in units of the non universal constant H) versus 
the scaling variable ^ (solid line). Experimental results for the auto correlation time in 



units of H exp for Fe (Ref. (TTJ: U) and Ni (Ref. M : 0, Ref. @: A, Ref. flIT3" | 



Figure 4.23: Temperature dependence of the /zSR damping rate for metallic Ni. The 
points are the experimental data of Nishiyama et al. ||205|| . The full line is the result of the 
model which takes the muon dipolar interaction into account. The dashed line gives the 
prediction when this latter interaction is neglected. 

Figure 4.24: Temperature dependence of the /iSR damping rate for metallic Fe. The 
points are the experimental data of Herlach et al. ||108|| . The curves are the predictions of 
mode-coupling theory for different sets of material parameters. In Fig. 4.24a (Fig. 4.24b) 
the muon is supposed to diffuse between tetrahedral (octahedral) sites. The full and dashed 
lines are present the results obtained with (q D ,£,o) equal to (0.033 A -1 , 0.82 Aj and to 
(0.045 A" 1 , 0.95 A) respectively; see Table [VITJ. 
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Figure 4.1: Typical scattering geometry in neutron scattering experiments. The mag- 
netic guide field Hy is used to define the beam polarization P. In order to measure the 
longitudinal spin fluctuations S L one has to go to a finite Bragg peak r with the wave vec- 
tor q perpendicular to the scattering vector Q = r + q. For unpolarized neutrons only the 
transverse fluctuations can be measured in the forward direction, r = 0. 
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Figure 4.2: Scaling function of the transverse Kubo relaxation function versus the 
scaled time variable r T at T = T c + 0.25K for q = 0.018A- 1 (solid), q = 0.025A- 1 (dotted), 
q = 0.0361" 1 (dashed), q = 0.071 A" 1 (long dashed) and q = 0.150A- 1 (dot-dashed). 
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Figure 4.4: Scaling function of the transverse Kubo relaxation function versus the 
scaled time variable r T at T = T c + 8.0K for the same set of wave vectors q as in Figs. 4.2 
and 4.3. 



116 




117 




Figure 4.6: Scaling function of the transverse Kubo relaxation function versus r' = 
Ai or <f t at T = T c for q = 0.025A' 1 (solid), q = 0.0361- 1 (dotted), q = 0.150A" 1 (dashed) 
and q = 0.300A _1 (long dashed). 
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Figure 4.7: Scaling function of the transverse Kubo relaxation function versus t' = 
A loi q z t at T = T c + 0.5K for q = 0.018A- 1 (solid), q = 0.025A- 1 (dotted), q = 0.036A- 1 
(dashed) and q = 0.07lA _1 (long dashed). 
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2.0 



Figure 4.8: Transverse Kubo relaxation function § T (q,g,t) at q — 0.024A -1 (solid line) 
and q = 0.028A -1 (dashed-dotted line) for dipolar ferromagnets versus time in nsec. The 
dashed line is the Kubo relaxation function for short range exchange interaction only at 
q = 0.024A" 1 . Data points q = 0.024A" 1 from Fig. 1 of Ref. ggD|. 
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Figure 4.9: Normalized experimental data deconvoluted by a maximum entropy 
method. The solid line shows the spectral shape function predicted by mode coupling the- 



ory The full circles (dashed line) indicate(s) the deconvoluted data. Taken from Ref. p3 
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Figure 4.10: Scaling functions versus y~ l = q/q D at T c (in units of the theoretical 
non universal constant A) for (i) the HWHM of the complete solution of the MC equations 
for the transverse Kubo relaxation function (solid), (ii) transverse (point-dashed) and (hi) 
longitudinal line (dashed) width in Lorentzian approximation. Experimental data of the 
transverse line width for EuO (U Ref. J28|, Ref. |THU[) and Fe (A Ref. pgfl). 
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Figure 4.11: Scaling function of the transverse line width in Lorentzian approximation 
(solid lines) and of the HWHM of the complete solution (dashed lines) versus the scaling 
variable x = (g0 _1 f° r a se ^ °f temperatures (T — T c = a) 1.4K, b) 5.8K, c) 21K, d) 51K). 
The Resibois-Piette function (solid) and the HWHM from the complete solution of the MC 
equations (dashed) without dipolar interaction is also plotted e). Experimental results for 
Fe from Refs. |TH||TH§ ( T - T c = (U) IAK, (0) 5.8K, (A) 21.0K, (o) 51.0K ). 
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Figure 4.12: Scaling plot of the temperature dependence of the relaxation rate of critical 
fluctuations in EuO above the Curie temperature T c = 69. 3K. The solid lines represent the 
results from mode coupling theory ||80|| . The data were taken between q = 0.018A" 1 and 
q = 0.15A- 1 . Taken from Ref. H§ Oi = 1/0- 
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Figure 4.13: Comparison of the longitudinal linewidth of Lorentzian fits normalized 
to A ex pg 5,/2 with A cxp = 2.1meVA 5 ^ 2 to the dipolar dynamic scaling function predicted by 
mode coupling theory for q D = 0.245A -1 for EuS. Figure taken from Ref. |33[ (k = !/£)• 
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Figure 4.15: Scaled peak positions q (Ai OT /uj) 2 ^ for constant energy scans of the scat- 
tering function for the complete solution of the MC equations versus the scaling variable 
Co for if = a) 1.490, b) 1.294, c) 0.970 and d) for the isotropic case, i.e., ip = 0. In- 
set: S T (q, lu)/S t (0, uj) in arbitrary units versus £ for if — 1.294 for some typical values of 
the scaled frequency (u = 10 L / 10 with L = 8 (solid), 10 (dashed), 12 (point dashed), 14 
(_.._.._) anc i ]_6 (_..._...—) indicated in the graph by crosses). 
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Figure 4.16: The same as in Fig. 4.15 for the scattering function resulting from the 
Lorentzian approximation. 
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Figure 4.17: Universal crossover function F(l/q u £)/F(0) with F(0) = 0.1956 for the 
Onsager kinetic coefficient at zero frequency and zero wave vector versus the scaling variable 



129 




Figure 4.18: Critical part of the Onsager kinetic coefficient for the homogeneous spin 
dynamics above T c of CdCr 2 Se 4 (Ref. p9fl ), CdCr 2 S 4 (Ref. pOfl ), £wO (Ref. pg ), 
£mS (filled circles: Ref. PI , open circles: Ref. M ), and iVi (Ref. |23§). The solid 
lines represent the result from mode coupling theory |JTJ. Inset: Nonuniversal ampli- 
tudes hLdq~ 3 / 2 for the kinetic coefficient and the relaxation rate at the critical point 
fiL q q~ l q~ l l 2 = T T (q, T = T c )/q 5 ^ 2 « 5.1326A for ferromagnets including Fe and Co. Figure 
taken from Ref. ||152| . 
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Figure 4.19: The real and imaginary parts of F are shown in a) and b) versus the 
scaled frequency Cj = u/Aq z D for several values of (p = arctan(g D £): ipi = 0.991 (top curve), 
y? 2 = 0.90f , cp 3 = f (1 - 1(T - 5 ), ¥» 4 = |(1 - 10-°- 25 ), V?5 = f (1 - 10- ai ), ¥>6 = |(1 - 10-°- 025 ) 
(bottom curve). 
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Figure 4.20: The scaling functions for the a) Onsager coefficient at zero frequency 
F(l/q D £) and b) the scaling functions u c2 (l/q D ^) and c) u c i(l/q D £) characterizing the large 
and low frequency behaviour, respectively. 
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Figure 4.21: Scaling functions Il,t{^p) f° r the transverse and longitudinal auto- 
correlation-time tl,t versus 1/<? D £. 
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Figure 4.22: Auto correlation time t c /H (in units of the non universal constant H) 



versus the scaling variable ^ (solid line). Experimental results for the auto correlation 



time in units of H exp for Fe (Ref. |TTB|: U) and Ni (Ref. 0, Ref. |§: A, Ref. |IT3 
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Figure 4.23: Temperature dependence of the /iSR damping rate for metallic Ni. The 
points are the experimental data of Nishiyama et al. ||205|| . The full line is the result of the 
model which takes the muon dipolar interaction into account. The dashed line gives the 
prediction when this latter interaction is neglected. 
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Figure 4.24: Temperature dependence of the /iSR damping rate for metallic Fe. The 



points are the experimental data of Herlach et al. ||108|| . The curves are the predictions of 
mode-coupling theory for different sets of material parameters. In Fig. 4.24a (Fig. 4.24b) 
the muon is supposed to diffuse between tetrahedral (octahedral) sites. The full and 
dashed lines present the results obtained with (g D ,£o) equal to (0.033 A -1 , 0.82 Aj and 
to (0.045 A" 1 , 0.95 A) respectively; see Table |VTTT| . 
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V. OTHER DIPOLAR SYSTEMS 



Dipolar interactions are of importance not only for the critical behavior of three- 
dimensional isotropic ferromagnets, but also for other magnetic materials such as antiferro- 
magnets, ferromagnets with an uniaxial or planar exchange anisotropy. Furthermore, dipolar 
interactions play an important role in two-dimensional systems, especially with regard to 
the existence of long-range order. 



A. Dipolar Antiferromagnets 

1. Hamiltonian and equation of Motion 

The order parameter of an antiferromagnet is the staggered magnetization, i.e., in the 
case of two sublattices the difference of the magnetization on these. Due to the alternating 
nature of the order parameter one supposes that the effect of the dipolar interaction on the 
long-range order parameter fluctuations is averaged out and thus one expects no influence 
on the static critical behavior. Indeed in Ref. || it was shown that in contrast to the ferro- 
magnetic case, the main parameters characterizing the dipolar interaction become irrelevant 
for the critical statics of antiferromagnets. Yet there is an important effect on the dynamics, 
since the dipolar forces lead to new damping processes in the long wave length dynamics 
of the magnetization, because they break conservation of the total magnetization. Since 
the magnetization and staggered magnetization are coupled dynamically, there should be a 
crossover in the dynamic critical exponents and related features in the dynamical scaling 
functions. This effect has been studied theoretically for temperatures above and at the Neel 
temperature T/v by means of mode-coupling theory |73| . 



The Hamiltonian of a dipolar antiferromagnet reads 

ff = -££(-W* + ^)sr^, (5.i) 

l^V a/3 

with spins S/ at lattice sites xj. The first term in brackets is the exchange interaction Ju> 
and the second the dipole-dipole interaction with the dipolar tensor given by 
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2 Vlx^-X;/^ |Xi-Xi/| 5 7 

We define Fourier transformed quantities by 

Sf = ^=£e iq - x S«, (5.3) 
4 = ^E^ (xrX,,) ^> (5-4) 

q 

< = ^E eiq ' (X!-Xi,) <> ( 5 - 5 ) 

q 

and get 

H = ~ E ( ^ + <*) ^q^-q ( 5 - 6 ) 

q 

for the Hamiltonian in terms of Fourier components. In order to study the long-wavelength 
behavior of the model one needs to know the behavior of the dipole-dipole interaction 
tensor at small wave vectors. The result depends on the lattice structure and the 
dimensionality of the system. For a three-dimensional system the expansion has been given 
in section III (see Eqs. ( |3.4| )). In two dimensions the dipole-interaction tensor has been 
analyzed by Maleev [Maleev76] . He finds for small wave vectors 

where q = {q x ,q y ,0) is an in-plane wave vector and z denotes the direction perpendicular 
to the plane r = (r x , r y , 0). The constants and A^ depend on the lattice structure 

^ (0) = |GW£4, (5.8) 

r^O 

A^ = -{g L ^ B )\ (5.9) 

where v 2 is the volume of the 2D unit cell. 

With the standard commutation relation for spin operators [S£, S@] = iUe^S^. one 



Af = A^ (^8^ - z a z^ + A™ q (z a zP - \ , (5.7) 



gets the following equations of motion ||217 
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d qx 



E 
k 



J kJ [ <Jk+q°-k "t" °k°q-k 



(^lf + ^k) (^k^q-k + ^q+k^-ky 

I A %y ( QX qz , qz qx \ 

"T^k \^ D k D q-k "T °q+k D -k J 

AXZ ( qx qy , Q2/ A 

/i k \^k q— k q+k — ky 

j_ av z ( qz qz I qz qy qy qy qy \ 

-|-^l k y q+k — k k q— k °k°q-k °q+k°-k,l 



(5.10) 



U. Critical Behavior of 3D Dipolar Antiferromagnets 



Now we turn to the critical dynamics of dipolar antiferromagnets. Specializing to the 
three-dimensional isotropic case and retaining only those terms which are dominant in the 



long wavelength limit the Hamiltonian for a simple cubic lattice reduces to |73 

m«mL 



H = J 



(go) 2 - 6) 5 al3 



q a q f} 



9- 2 



+jj (Q-(qa) 2 ) 6 a ?N«N^, 



(5.11) 



where M° = denotes the magnetization and iV? = the staggered magnetization, 

and we have introduced the abbreviation / q = J v a -£^k. The wave vector q = ^(1, 1, 1) 

characterizes the antiferromagnetic modulation. 

Starting from this Hamiltonian mode coupling equations for the correlation functions of 

the magnetization and staggered magnetization have been derived and analyzed in Ref. 

For instance for the longitudinal magnetization mode one finds 

i 



<1BZ 



^) 2 sin 2 ^(k + ,e,t)^ / (k,t) 



r^(q,t) = 4A 2 f J k 2 dk J d(costf) 
o -1 

(2q ' k + (sin 2 tfd>£(k + , t)^(k, t) + cos 2 ^(k + , t)*£(k, t 



+ 



( k 2 + e- 2 )(k 2 + +^ 

where we have used the abbreviations k + 



(5.12) 



k + q, A af = a 3 / 2 y/lffi and q BZ 
It is found that the memory kernels obey generalized scaling laws of the form 



2tt 3 /_3_ 
2a V 4?r ' 
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t m T n{^ 9, £ t) = A 2 ai q 2z ^ T N (x, y, r) , (5.13) 
^M T N(l:9A,t) ^ (f}^ N (x,y,T) , (5.14) 

where the scaling variables are defined by x — l/q£, y — qA/q, and r = A a {q z t with the 
dynamic critical exponent z = 3/2. The characteristic wave vector is related to the 
dipolar wave vector q D = yfg/a by 

<?A=(-J (q D af 3 q BZ . (5.15) 

The value 3/2 is the dynamical critical exponent for isotropic antiferromagnets and any 
crossover to dipolar behavior is contained in the dynamic scaling function. In evaluating 
the transport coefficients V^ N only two mode decay processes have been considered (see 
Eq. (4.12)). Without dipolar interaction, i.e., for the isotropic exchange antiferromagnet, 
the magnetization modes can decay into staggered modes only ||140| , |263| , |231| , |129 . 



The dipolar interaction leads in addition to a decay into two magnetization modes, 
dominating for small wave vectors. In the strong dipolar limit (T — ► Tjy and q — > 0) the 
mode coupling equations for the magnetization modes can be solved exactly with the result 



r M (9>£>*) = /1A a Qbz 



9 8 9/ 



S(t) (5.16) 



corresponding to a relaxation of the magnetization modes enforced by the dipolar interaction. 
The phase space for this decay is the full Brillouin zone, whereas the decay into staggered 
modes is weighted by the critical static susceptibilities. Therefore, with separation from the 
critical temperature and for larger wave vectors the decay of the magnetization mode into 
two staggered modes dominates the decay into two magnetization modes. 

The mode coupling equations for dipolar antiferromagnets have been solved in the 
Lorentzian approximation in Ref. ]73]| . The results for the transverse scaling function for 



the magnetization and staggered magnetization 7^- N (x, (p) versus x for different values of 
ip = Ntt/30 with N — 0, 1, 14 are displayed in Figs. 5.1 and 5.2. The corresponding plots 
for the longitudinal scaling functions are very similar J73|. The curves with N = corre- 
spond to the isotropic case. It is important to note, that the effect of the dipolar interaction 
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on the scaling functions of the magnetization and staggered magnetization modes is rather 
different. The scaled line width for the magnetization modes exceeds the isotropic curve by 
an amount which increases with T on approaching T c (i.e. larger values of <p = arctan(g D £)). 
This reflects the crossover from a diffusive behavior Tm oc q 2 q^ 2 to a relaxational behav- 
ior Tm oc q^ 2 • The line width of the staggered magnetization at fixed scaling variable x 
decreases with T approaching T c since the magnetization modes become uncritical. But, 
this change is much less pronounced than for the magnetization modes and more impor- 
tant, the asymptotic hydrodynamic dependence on the wave vector and correlation length 
is unmodified, i.e., the hydrodynamic behavior is always relaxational. 

In the dipolar (D) and isotropic (I) critical (C) and hydrodynamical (H) limiting regions 
the mode coupling equations can be solved analytically. These regions are defined by DC: 
y 3> l,x <C 1; IC: 1/ < 1, i < 1; DH: y ^> x, x ^> 1;IH: j/ < 1, 2; > 1. The results are 
summarized in Table |X|, where we note that these asymptotic power laws are the same for 
the longitudinal and transverse fluctuations. In the immediate vicinity of q = q A there is 
a dynamic crossover from z e $ = 3/2 to z c q = 2 for the staggered magnetization and from 
z e ff = 3/2 to an uncritical value z e g = for the magnetization modes. 

RbMnF 3 is one of the most thoroughly studied isotropic antiferromagnets [255|j5T|jr52] [ . 
Taking the values for the exchange coupling and the static susceptibility of the magnetization 
from Refs. ||255| , [270| , |121|| , the characteristic wave vector q A is about 0.02A" 1 . Unfortunately 
the neutron scattering data [255 , |51| , 162 for the line width of the staggered magnetization 
modes are limited to the isotropic region and to our knowledge there are no data for the line 
widths of the magnetization modes. 

The dipolar effects could be observed more readily in EuSe and EuTe because of their 
smaller transition temperatures 4.8K and 9.7 K implying larger q A . Experiments in the 
appropriate wave vector and temperature region are, however, still lacking. 
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B. Uniaxial Dipolar Ferromagnets 



We have mainly concentrated on the influence of dipolar forces on otherwise isotropic 
ferromagnets. But, even in systems which are anisotropic to start with such as uniaxial 
ferromagnets, the dipolar interaction has a significant effect, because of its long-range nature. 

The influence of the dipolar interactions on the critical statics of isotropic and uniaxial 
ferromagnets is quite different. Whereas for isotropic ferromagnets the dipolar interaction 



leads only to a slight modification of the critical exponents, Larkin and Khmelnitskii ||161 



have discovered that uniaxial dipolar ferromagnets show classical behavior with logarithmic 
corrections in three dimensions. This system was then studied by Aharony || by means of 
the renormalization group method. Aharony employed an isotropic elimination procedure 
and calculated the critical exponents to first order in (3 — d) as well as logarithmic correc- 
tions @. Note, that the dipolar interaction leads to a shift in the upper critical dimension 
from 4 to d c = 3. 

The existence of logarithmic corrections was verified experimentally for a number of 
uniaxial ferromagnetic substances |||90|]. However, these experiments were performed in 
regions of the reduced temperature, where departures from the asymptotic behavior are 
expected and are observed indeed. In particular, a maximum in the effective exponent of 
the susceptibility |5D| has been found. On the basis of a generalized minimal subtraction 
scheme the latter crossover from Ising behavior with non classical exponents to asymptotic 
uniaxial dipolar behavior, which is characterized by classical exponents with logarithmic 
corrections, has been analyzed in Ref. [ S3 ,[S8| . The theoretical results for the specific heat 
and susceptibility are in quantitative agreement with measurements on LiTbF^, where solely 
the strength of the dipolar coupling constant entered as an adjustable parameter |83|j88|| . 

The Landau-Ginzburg free energy functional for a n-component uniaxial spin system 
with an isotropic exchange coupling and dipolar interaction is given by ||161 ,B 

1 



H - 2 



r + k 2 + g 2 ^ 2 



S a (k)S a (-k) 
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J ^(k 1 )5 Q (k 2 )5 /3 (k 3 )5 /3 (-k 1 -k 2 -k 3 ) 



(5.17) 



Here ^"(q) are the n components of the spin variables. The (^-dimensional wave vector 
k = (P) <?) is decomposed into q, the component along the uniaxial direction, and p, the 
remaining (d-1) components. The bare reduced temperature is given by r = (T — T c °)/T c ° 
and g 2 is a measure of the relative strength of the dipolar interaction. The dipolar term 
g 2 $2 suppresses the fluctuations in the ^-direction. 



In Ref. |76j the dynamics of uniaxial dipolar magnets was described in terms of a simple 
Time Dependent Ginzburg Landau Model (TDGL) with the equation of motion 



l s ° (k - ( ) = - r »^7) + " (k - ( )' 



(5.18) 



Since the order parameter is a non-conserved quantity in the presence of the dipolar in- 
teraction, the kinetic coefficient To is wave vector independent. The random forces T](k, t) 
represent the uncritical degrees of freedom and are characterized by a Gaussian probability 
distribution. A renormalization group analysis for the above TDGL gives for the dynamic 



critical exponent [|7g] 



z = 2 + ct] , with c = 0.92 



(5.19) 



This may be compared with the value for the TDGL model for spin systems with short 
range exchange interaction | |100| , |101| , where c = 0.73 and the van Hove prediction c = — 1. 
Since for uniaxial dipolar magnets the fourth order coupling is marginal at d — 3, there 



are logarithmic corrections to the classical behavior. For the statics one finds ||161| ,p| for the 
susceptibility x an d the specific heatC 



X oc 



r^llnrl 1 / 3 , Coc|lnr| 1/3 , 



(5.20) 



where r = (T — T c )/T c is the reduced temperature. The effect of the logarithmic corrections 
on the dynamics susceptibility can be approximated by |76 



X{uj,p,q,T) 



^ 9 Q , N I 1 , / N I 1 It. 

-Tr+P 2 + g^ + m(o) | rrin/^o) r 1 / 3 



2/ 



(5.21) 
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where /xi(0) is proportional to the reduced temperature r, and Zo is some constant. In the 
hydrodynamic region q <^ip £ _1 the relaxation coefficient implied by the latter result is 

T oc r (| lnr | +const.) _1/3 . (5.22) 



This result has also been found by Maleev |177| ] using mode coupling arguments. A detailed 



theoretical study of the dynamics of uniaxial dipolar ferromagnets, which takes into account 
not only the relaxational dynamics on the basis of TDGL model, but also includes nonlinear 
terms resulting from the Larmor precession in the local magnetic field is still lacking. A 



mode coupling analysis of this problem is in progress ||107|| . 



C. Two dimensional Systems 

Two-dimensional systems are interesting, because matter behaves qualitatively different 
as compared to three dimensions. We recall that the Bloch argument |2l!| , |211| , |158|| , which is 



proven on a rigorous basis by the Hohenberg-Mermin- Wagner theorem [ [L11|J186|| , excludes 
conventional long-range order in isotropic systems with a continuous symmetry and short- 
ranged interaction such as the Heisenberg ferromagnet, superfluid Helium films, and two- 



dimensional crystals [|202| . Any hypothetical broken continuous magnetic or translational 



symmetry would be overwhelmed by long wavelength spin or phonon excitations, since the 
phase space for these Goldstone-like modes is enhanced in two dimensions. Two dimensions 
seem to be the borderline dimension where thermal fluctuations are just strong enough to 
prevent the appearance of a finite order parameter. 

However, since the dipolar forces reduce the fluctuations longitudinal to the wave vector, 
related to the anisotropy of the dipolar interaction with respect to the separation of two spins 
one expects that a finite order parameter could exist in two-dimensional isotropic systems. 
In fact we will review in the following spin-wave theories, which show that the dipolar 
interaction leads to the existence of long-range order in two-dimensional ferromagnets and 
antiferromagnets and other systems which belong to the same universality class. 
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1. Ferromagnets 



The possibility of a finite order parameter in two-dimensional ferromagnets was shown 
first by Maleev [178], who evaluated the low temperature properties, in particular the 



magnon excitation spectrum for Heisenberg ferromagnets in the presence of dipolar forces. 



Linear spin-wave theory [|1 78|[ results in a magnon dispersion relation 



£ k = y/(Dk 2 + n a) (Dk 2 + tt ka sin p k ) , (5.23) 

where Dk 2 = 2S(Jq — Jk) at ka < 1, flo = 2nS(gLfiB) 2 /v2a <C 2SJo, a = SA Q /Q ~ 1, and 
(fk is the angle between the wave vector k and the magnetization M. At sufficiently small 
k (and for k not parallel to the magnetization) the energy of the spin-waves has the form 

oc f2 | sin (p^ \Vk. (5.24) 

An analysis of the relative deviations of the magnetization from saturation (in the framework 
of linear spin-wave theory) results in an order of magnitude estimate for the transition 
temperature (for small dipolar couplings) 

It was also shown |109|j63| , [T6|] that an exchange anisotropy leads to a suppression of thermal 



interaction and in turn to the existence of long-range order two-dimensional ferromagnets. 

In two-dimensional ferromagnets, the dipolar interaction results in an easy-plane 
anisotropy ||178| , [223|| . This is due to the fact that the dipolar energy is minimized when all 



spins are aligned in plane. Thus the problem of a three-component Heisenberg ferromagnet 
with dipolar interaction is reduced to the problem of a two-component (planar) ferromagnet 
with renormalized values of the temperature and of the dipole-dipole interaction constant. 
The magnetic dipolar forces change fundamentally the nature of the low-temperature phase. 
This can be seen as follows. In the absence of dipolar interaction any long-range order is 
impossible since it gets destroyed by thermal fluctuations; the mean square of the spin fluc- 
tuations is given by (| SSq | 2 ) oc k B T/q 2 . A planar (pure) exchange ferromagnet would 
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undergo a Kosterlitz-Thouless ||145| , |146|| transition. The Kosterlitz-Thouless phase is char- 
acterized by local but no long-range order and by bound vortex pairs. There is a divergent 
susceptibility throughout this phase. However, due to the dipolar interaction the spectrum 
of the spin fluctuations at small wave numbers is changed. Including the dipolar interaction, 
Eq. ( ]5.7p , one obtains the following reduced Hamiltonian 



n 



H 



JS 2 



v 2 a 



d 2 q 



k B T k B T ^ J (2tt) 2 
in terms of the azimuthal angle 0, where we have assumed a fixed length S 2 



(5.26) 



1 and 

an ordering along the x-direction. The relative strength of the dipolar interaction in two 
dimensions is characterized by the parameter g 2 = A^ / Ja 2 . In the derivation of Eq. (|5.26|) 
the in-plane components of the dipolar interaction tensor have been used ||214|| = 
A^gap _ A (i)££ + Q(q 2 ) ( A and A 1 see Eqs. (pj) and (p|)). Neglecting all but the 
terms quadratic in q , the mean square fluctuations about the local ordering are given 
by 1213 

k R T f d 2 q 1 



J J (2tt) 2 q 2 + g 2 \q\ sin ip c 



(5.27) 



where ip^ is the angle between q and the total magnetization. The above integral no longer 
diverges due to the presence of the dipolar interaction. Hence the dipolar interaction stabi- 
lizes ferromagnetic long-range order in more than one dimension ||221|,|178 |. Besides stabi- 



lizing the long-range order, the dipolar forces alter the vortex behavior ||214|| . Ref. |[214]| also 
contains an analysis of the static critical behavior. More recently the critical behavior has 
been reanalyzed in Ref. |18| . In passing we note that there are two reviews on excitations 



in low-dimensional systems by Pokrovsky et al. |[222| , [224 



Actually the situation is more complicated in thin magnetic films. Mono- and two-layer 
iron orders perpendicularly to the planes for T = OK, only for layer numbers above 3 the 



expected parallel orientation is observed |TTJ. This comes from the spin-orbit interaction, 
which effectively leads to a perpendicular anisotropy, which overcomes the dipolar inter- 
action for the very thinnest layers. Upon raising the temperature a reorientation phase 
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transition ||206|| from a perpendicular to an in-plane orientation of the magnetization has 



been observed experimentally. It was argued by Pescia et al. ||212| , [213[| that this transition 
is driven by the competition between spin-orbit and dipolar interaction. The (classical) 
transition temperature has been estimated ||225| , [213|| to be of the order Tr ~ (A — fl) 7 ^^2 , 
where c is the number of nearest neighbors, Q = 2ir(gLiiBS) 2 /a 3 characterizes the strength 
of the dipolar interaction, and A is the single-ion anisotropy constant favoring perpendicu- 
lar orientation. Since the single-ion and the dipolar anisotropy parameters scale the same 
under a renormalization group transformation [|213| , [225|1 , it depends on the relative values 
of the single-ion anisotropy and the dipolar parameters, whether the transition tempera- 
ture for the reorientation phase transition is smaller or larger than the Curie temperature. 
However, since the dipolar interaction parameter increases with the number of layers, the 
value of Tr is reduced for thicker films, as indeed observed experimentally |[206|| . Some of 



the above theoretical work is still rather controversial ||164|| and alternative mechanisms have 
been proposed for the experimentally observed reorientation phase transition ||163|| . Due to 
the combination of various anisotropics, dipole-dipole interaction, reduced dimensionality 
and enhanced importance of thermal fluctuations the phase diagram of thin magnetic films 
shows a variety of new and interesting phases | |206| , |iT||rT] with properties still awaiting a 
theoretical explanation. 

2. Antiferromagnets 



Recently it has been shown [|215|| that long-range order is also possible in two-dimensional 



antiferromagnets due to the anisotropy of the dipolar interaction. The existence of the long- 
range order is a consequence of a subtle interplay between exchange and dipolar interaction. 
The classical ground state of an isotropic pure exchange antiferromagnet has a continuous 
degeneracy and hence no long-range order. The dipolar interaction lifts the continuous 
degeneracy of the ground state such that a spin alignment perpendicular to the plane is 
energetically favored (see Fig. 5.3). In other words, whereas the exchange interaction imposes 
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the antiferromagnetic order, the dipolar interaction prevents thermal fluctuations from its 
destruction. This can be seen most easily by considering the equations of motion at zero 
wave vector. Upon approximating the longitudinal spin components by their equilibrium 
expectation value Sf ~ S , e* qo ' X! the linearized form of Eqs. ( |5.10| ) (see Section |5TT| ) reads 

d 



S x = 2S (A£ - A") fl* , (5.28) 
( J qo - Jo) + (AZ - A™)] S$ , (5.29) 



dt 

s s « = 2S 



and an analogous set of equations for Sq and S*. The wave vector q = ^(1, 1,0) charac- 
terizes the antiferromagnetic, staggered modulation of the ground state. The coefficient on 
the right hand side of Eq. (|5.28|) is nonzero due to the anisotropy of the dipolar interaction 



in 2D. The resulting energy gap at zero wave vector then is 



E = 2SJA" - A% v/(J qo - Jo) - (Ar - Ag) • (5.30) 



A more detailed information on the spin-wave spectrum can be obtained by linear spin-wave 
theory ||215| . Introducing Bose operators by a Holstein-Primakoff transformation [141], given 



here only up to harmonic terms, 

Sf = \/f(^W), S? = TiM(ai-aS) , Sf = ±(S - aM) , (5.31) 

where the upper (lower) sign refers to the first (second) sublattice, the Hamiltonian in the 
harmonic approximation is given by 

H = Eq {A} Oqflq + T^q («q «-q + ftq^-q) + 

C q a q Ct-q-q + Cq a q a -q-qo + C q a q a q+qo + C* CLq+q^q} (5.32) 

with the coefficients 

A q = S(2J CLo — J q — J q+qo ) + S(2A z ^ o — A x * — A q + qo ) (5.33a) 
fi q = 5(J q+qo - J q ) + S« +qo - A™) (5.33b) 
C q = iSA x v. (5.33c) 
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In this description the primitive cell is the crystallographic, which is half the magnetic. 
Diagonalization of the Hamiltonian, achieved by a Bogoliubov transformation, results in a 
spin-wave spectrum with two branches 

4 = ^(fixifia) (5.34) 

where 
and 

^2 = (^q ~ — ^4 q + qo + -Sq+q ) 2 + 16 [C q+qo ( A q+q() — -B q+qo ) — C q (Aq — Bq)] 



><[Cq (^.q+qo + -^q+qo) ~~ Cq+q (^q + -^q)]- 

In Fig. 5.4 the dispersion relation is shown for three values for the ratio of dipolar and 
exchange energy k = with isotropic nearest-neighbor exchange interaction (J < 0). 

The two branches can be resolved only for large values of k. 

In particular the gap is proportional to the square root of the difference of the magne- 
tostatic energy between the configurations of in-plane and out-of-plane magnetization. In 
a three-dimensional simple cubic lattice the first root in Eq. (|5.30|) vanishes because of the 
symmetry, but in two-dimensional systems there is a finite gap for perpendicular antiferro- 
magnetic order. We note that for sufficiently large exchange energy the gap is the geometric 
mean of dipole and exchange energy, which in turn implies that the gap is much larger than 
the dipolar energy for k < 1. 

The Neel temperature Tjv has been determined via a high-temperature expansion, by 



employing Callen's method ||2X5| — 12 17| | . If the dipolar interaction is weak in comparison with 
the exchange interaction one can give an order of magnitude estimate for the transition 
temperature in the framework of linear spin-wave theory based on the Holstein-Primakoff 
transformation 
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J 




In 


OS) 



T N oc — ttttv • (5-35) 



This, however, gives an overestimate for the transition temperature, because it uses a 
temperature-independent dispersion relation. Actually, the magnon frequency softens with 



increasing temperature. Those effects have been accounted for in Ref. ||216|| by an extension 
of the Tyablikov decoupling scheme due to Callen In essence this leads to a replacement 
of the saturation magnetization of the spins S by the temperature-dependent order param- 
eter a. The resulting transition temperature is lowered with respect to the estimate from 
linear spin-wave theory. Eq has been obtained by linear spin-wave theory and Eq and T^ 1 
by the method of Callen. The results are summarized in Table |X| and show a quite satisfac- 
tory agreement with experimental data. This theory has been extended to antiferromagnets 



on a honeycomb [218] and several other lattices [219 
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Tables: 

TABLE X. Behavior of the scaling functions of the line widths of an isotropic dipolar antifer- 
romagnet in asymptotic regions. 



IC 



DC 



IH 



DH 



L,T 
L,T 

In 



,3/2 



,-1/2 



3-1/2 
x 3/2 



r 3/2 
r 3/2 



TABLE XL Exchange energy |J|, lattice constant a, energy gap Eq, spin-flop field H s {, Neel 
temperature Tjv and zero temperature order parameter gq. 





\J\ 


a 


^0 


^0 


Eq 


rpexp 
N 




0"o 




[K] 


[A] 


[K] 


[K] 


[K] 


[K] 


[K] 




K 2 MnF A 


8.5 a 


4.17 a 


7.4 b 


7.1 


7.6 


42 a 


41 


2.33 


Rb 2 MnF 4 


7.4 C 


4.20 g 


7.3 b 


6.5 


7.0 


38§ 


36 


2.33 


Rb 2 MnCk 


11. 2 f 


5.05 e 


7.5 f 


6.1 


6.6 


56 f 


48 


2.32 


(CH 3 NH 3 ) 2 MnCl 4: 


9.0 f 


5.13 e 




5.3 


5.7 


45 f 


39 


2.32 



a Reference p4]| , b Reference [EBS], c Reference [266] , d Reference ||267|| , e Reference [ |16C 
Reference 1124311 , 9 Reference fE3l . 
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Captions to the figures: 



Figure 5.1: Scaling function for the line width of the transverse magnetization above T/v 
versus x — l/g£ for different values of if = Nn/30 (N=0,l,2,...,14). 

Figure 5.2: Scaling functions for the line width of the transverse and longitudinal 
staggered magnetization above Tn versus x = l/q£, for different values of <p — Nir/ 30 
(N=0,l,2,...,14). 

Figure 5.3: Classical ground state of a two-dimensional dipolar antiferromagnet with a 
dipolar interaction much smaller than the exchange interaction on a quadratic lattice. Taken 



from Ref. 217 . 



Figure 5.4: The spin-wave dispersion relation (Eq. ( |5.34| )) of pure exchange antiferromag- 
nets on a quadratic lattice with nearest-neighbor interaction (solid line) and with additional 
dipolar interaction (S = 1/2), for the ratios of dipolar energy to exchange energy k = 
along the -[£,£> 0] direction: k = 0.1 (long dashed), k = 0.01 (dashed) and k = 0.001 
(dotted). The splitting of the two magnon branches is visible only for k = 0.1. 
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N=0 



10" 2 10" 1 10° 10 1 

x=lA£ 

Figure 5.1: Scaling function for the line width of the transverse magnetization above 
T N versus x = l/q£ for different values of tp = Nir/30 (N=0,l,2,...,14). 
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Figure 5.2: Scaling functions for the line width of the transverse and longitudinal 
staggered magnetization above T N versus x = l/q£ for different values of ip = Nn/30 
(N=0,l,2,..,14). 
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Figure 5.3: Classical ground state of a two-dimensional dipolar antiferromagnet with 
a dipolar interaction much smaller than the exchange interaction on a quadratic lattice. 



Taken from Ref. 217 
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0.5 



Figure 5.4: The spin-wave dispersion relation (Eq. (|5.34 )) of pure exchange anti- 
ferromagnets on a quadratic lattice with nearest-neighbor interaction (solid line) and with 
additional dipolar interaction (S = 1/2), for the ratios of dipolar energy to exchange energy 



K 



4\J\a 3 



along the -[£, £, 0] direction: k = 0.1 (long dashed), k = 0.01 (dashed) and 



k = 0.001 (dotted). The splitting of the two magnon branches is visible only for k = 0.1. 
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VI. STOCHASTIC THEORY 



A. Classical field theory and dynamic functional 

If one tries to describe the critical dynamics of ferromagnets, one is in general not in- 
terested in the complete, highly complicated, microscopic time evolution. Rather, one is 
usually interested in the dynamics on time scales characteristic of its slowly varying dynam- 
ical quantities, such as those associated with hydrodynamic, order parameter or Goldstone 
modes. 

A method to extract from the microscopic dynamics the equations of motion appropriate 
to these slowly varying quantities is the projection operator formalism of Zwanzig | |276|| and 



Mori |]1971| . Since the semi-phenomenological equations derived from this formalism form 
a common basis of mode coupling and dynamic renormalization group theories we briefly 
review the main ideas and results of this approach. 

The idea is to eliminate the fast variables by introducing a projection operator V, which 
projects onto the subspace of slowly varying modes. Descriptions of this formalism can be 
found in Refs. ||139| , |198|| . A short account is given in Appendix B. After having chosen an 



appropriate set of slow variables - which is crucial for the validity of the dynamics - one can 
derive formally exact non linear equations of motion for these modes (here these modes are 
denoted by {S a (t)}) 

j t s a (t) = v a (is(t)}) 

+E - T)} f_M^mgmw^m + ut) , (6 . 1} 



where 



P eq ({S}) = exp 



1 -M{S}) 



k B T 

is the equilibrium probability distribution function. The first term in these equations. 



(6.2) 



v a ({S(t)})=<tCS a ,{S(t)}>, (6.3) 
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is called the streaming term and describes the systematic part of the driving forces. This 
term can be expressed in terms of Poisson brackets. (The conditional average < X; {a} > 
is defined by < X; {a} >=< X5(S — a) > /P eq ({S(t)} where the average is taken in the 
microcanonical ensemble ||140|| .) The second term contains the memory kernel 



M a(3 (t; {S(t)}) =< Ut)Q(0); {S(t)} > , (6.4) 

and characterizes the damping of this adiabatic motion by frictional effects arising from 
the "random" forces. The "random" forces are given in terms of the projection operator 
Q — 1 — V and the Liouville operator C 

( a (t)=exp[itQ£]iQ£S. (6.5) 

Whether this force can be regarded as random or not depends crucially on the choice of the 
projection operator, i.e., whether V contains all the slowly varying quantities of the system 
under consideration. 

The semi-phenomenological equations of motion are obtained from the above exact equa- 
tions by making the following basic, and at least plausible assumptions. First one assumes 
that one can neglect memory effects in M a p and makes the Markovian approximation 

M Ql3 (t; {S(t)}) « 2L Qf3 ({S(t)})5(t) . (6.6) 

Second one takes the kinetic coefficients to be independent of the slow variables 

L af3 ({a}) « L af3 , (6.7) 

and third one assumes that all slowly varying variables are contained in the choice of the 
projection operator so that the forces f a are really random and have a Gaussian probability 
distribution 

1 /•*! 



w({C}|io<t<ti) ~exp 



dtUt)L- a lc {t) 

to 



(6.8) 



4 

Then the exact equations reduce to the following semi-phenomenological equations of the 
Langevin form 
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j t S a {t) = v a ({S(t)}) - £ L% 6H ^ }) + Ut) , (6.9) 



where the streaming term, 

Va({S(t)}) = -A/£ \^-Q aP ({S}) - Qa P ({S}) 5 -^^- , (0.10) 

can be expressed in terms of Poisson brackets 

QMS}) = {S a , Sp} P = -Q Pa ({S}) . (6.11) 

The equilibrium distribution does not determine the driving force in the Langevin equa- 
tions uniquely. Only the dissipative couplings which are generated by the derivative of the 
effective Hamiltonian 7i are related to the static properties of the system. As a direct conse- 
quence with each static universality class are associated several dynamic universality classes 
depending on the structure of the Poisson brackets. 

After these general preliminaries we can now start to set up a stochastic equation of 
motion for the dynamics of dipolar ferromagnets. If the interaction between the spins in 
an isotropic Heisenberg ferromagnet is solely given by the short range exchange interaction 
the dynamics are described by precession of the spins in a local magnetic field generated by 
the surrounding spins. The rotation invariance of the Heisenberg Hamiltonian implies that 
the total spin is conserved. As shown explicitly in chapter A| this is no longer the case if 



one takes into account the long-ranged dipole-dipole interaction. Hence the spin dynamics 
of a isotropic dipolar ferromagnet are described by the following semi-phenomenological 
Langevin equations 

as(x,t) gg[gj) f mm , „ , fi19 . 

= A/S(x ' t] x " L w^t) + c(x ' t] ' (6 - 12) 

where S(x, t) is the local spin density. The local magnetic field is determined by the func- 
tional derivative of the Landau-Ginzburg functional T~i[S\ with respect to the spin density. 
For the present case the Landau-Ginzburg functional is given by 



nS} = -\j (r + q 2 )S^ + g^ S"(q)tf»(-q) 

q L ^ - 

~T\ III sa ^) sa ^) sP ^) S ^-^ " * - <b) • (6-13) 



qi 92 93 
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S a (q) (a=l,2,...,n) are the components of the bare spin variable with n equal to the space 

dimensionality d. We further used the abbreviation J = J d d q/(2ir) d , r is the reduced 

g 

temperature-deviation and g denotes the relative strength of the dipolar interaction. 

The dipolar interaction in Eq.( |6.13D breaks the symmetry of the spin fluctuations trans- 



verse and longitudinal to the wave vector q, which is reflected in the free propagator 

<r \ <r J 

where G L (r, g,q) = (r + g + q 2 )' 1 and G T (r, q) — (r + q 2 )~ l ■ 

The time scale is characterized by the quantity A. The mode coupling coefficient / 
determines the strength of the coupling of the spin density to the local magnetic field. Since 
the order parameter is a non conserved quantity the leading term in the Onsager operator 

L = A( 7 -V 2 ) (6.15) 

is wave vector independent, A7, and not — AV 2 as for a conserved order parameter. The 
quantity 7 characterizes - like g in statics - the relative strength of the dipolar interaction. 
The random forces C( x >^0 have a Gaussian probability distribution 

u> f ([CiCM)] I t < t < h) ~ exp (^--^J^ dt J d rf <i(x,*)£ _1 CtCM)) > (6-16) 

which is uniquely determined by the lowest moments 

(C(x,t))=0, (6.17a) 
(C»(x, t)0(x, 0) = 2L5(x - xT)5(jt - if)S i5 . (6.17b) 

Note that we have set /3 = 1/ k B T = 1 . 

In particular for the implementation of the dynamic renormalization group theory in a 
way analogous to static critical phenomena it is convenient to introduce a functional which 
generates the perturbation expansion for the frequency dependent correlation and response 
functions, which is equivalent to the equations of motion. We shall use the functional integral 
formulation [|127|Jl7| , which converts the Langevin equations into a dynamic functional with 
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one additional field [182]. The idea is that, instead of solving the Langevin equations (|6.12|) 
for the stochastic spin fields S(k, t) in terms of the random forces £(k, t) and then averaging 
over the Gaussian weight, one can eliminate the random forces in favor of the spin variables 
by introducing a path probability density W^dS 1 }) via 



W{{S})V[S] = w({(})V[C} . 



(6.18) 



Furthermore, it is convenient to perform a Gaussian transformation in order to "linearize" 
the dynamic functional. This is accomplished by introducing response fields S | 127 , 182|1 by 



W({S}) = J V[iS) exp { J[S, S 



(6.19) 



For more details on the general formalism we refer the reader to Refs. [|127] , |6ll , |l"7| , |182|| . The 
generating functional for the correlation functions is 



Z[h,h] = jf I V[S}V[iS}exp 



Jll [S, S] + £ / d d x f 1 dt (hiSi + KS,, 

4 J Jto 



(6.20) 



where the normalization factor M is chosen such that Z[h — 0,h — 0] = 1. The symbols 
T>[S] and T>[iS\ denote functional measures of the path integral. The Janssen-Dominicis 
functional is given by 



J£[s,§\ = E/ dt dd% 



to 



where the functional of the forces is 



\ dt 11 2 6Si 



(6.21) 



(6.22) 



<SS(x,t) 5S(x,t)' 

In formulating the perturbation theory one splits the functional, Eqs. ( |6.21| ) and ( |6.22 ), 
into a harmonic part 



' k Jld 



\(-f + k z )S a (k,u)S a (-k,-Lu) + 



lit) 



5 aP + X(-f + k 2 )G a ^(r,g,k)- 1 



S a (k,cj) 



(6.23) 



which is bilinear in the stochastic spin fields S and S and a part J- xrA [S, S] containing the 
interaction terms. Here G af3 (r, g,k) denotes the static propagator, Eq. ( p.!4| ). Further- 
more we have introduced the short hand notations J k = J d d k/(2Tr) d , § w = j duj/2ir. The 
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Fourier transform of the spin density is defined by S(x, t) = J k j' w S(k, a;)e^ k ' x wt ) . Hence 
the generating functional can be written as 



Z[h,h] = ^=exp 



Jmt \Sh'5h / 



x / V[S]V[iS}exp 



JtoUS, S]+J2j ddx £ dt faSt + hSi) 



(6.24) 



We start with the discussion of the harmonic part, Eq. ( p. 23| ). It can be written in matrix 
notation as 



-2L(k) iu - A a ?{k) 

-iuj - A a P(k) 



A af3 (k) = -L(Jfe)[G a/3 (r>ff> k )]~ 1 > 



Jharm[5, S] = ~\ J J S"(k, uj) A^(k, u) (§ P {-k, -(j), S^-k, -u)f . (6.25) 

The matrix A a/3 (k, uj) is given by 

A«*(k,u;) = | ' I. . | , (6.26) 
where we have defined 



(6.27) 



L(k) = A( 7 + fc 2 ). 



The harmonic part of the partition function 



Z [h,h] = / V[S\V[iS\exp 



J^rm[S, S] + £ / d d X dt (hiSi + kS, 



can be calculated explicitly with the result 



Z [h, h] = exp 



1 

2 Jk 



~h a {k, u), h a (k, u)) (A a/3 (k, lu))- 1 (hP{-k, -u), h p {-k, 



(6.28) 



(6.29) 



-uj 



T 



(6.30) 



Hence the free propagators are found to be 
(S a (k, u)SP(k>, uj'))o (S a (k, u)SP(U, uj')) ' 
(S a (k, u)S^k', */))„ (S<*(k, u)S?(k', «0>o , 



6(oo + u')5(k + k') (A a/3 (-k, 



-uj 



(6.31) 
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where 

(A^(k,^)) _1 = A- L \k,u;)Pf (k) + A T 1 (k,^)P^(k) . (6.32) 



Hence the longitudinal and transverse propagators are given by 

/ - A a {k)\ \ 

A; 1 (k,^)= (6.33) 
\l/[iu-A a (k)} 2L(k)/[co 2 + A a (k) 2 ]) 

where 

!r + k 2 for a = T, 

(6.34) 
r + g + k 2 for a = L. 

In summary, one gets for the response propagators 

R ap {k, u) = (S Q (k, u)S p (-1t, -u)) = R T (k, uj)Pf(k) + R L (k, u)Pf{k) , (6.35) 

where the transverse and the longitudinal parts are given by 

R L > T (k, uo) = — \ — , (6.36) 

-ico - A L , T {k) 

and the projection operators are defined by P£ — <5 a/3 — q a q f3 /q 2 and P£ = q a q /3 /q 2 . For 
the correlation propagators one obtains 

C a/3 (k, u) = {S a (k 1 uj)S^(-k, -w)) = C T (k, oj)Pf{k) + C L (k, oj)Pf{k) , (6.37) 

where 

CW {k , u) = J 2 y ) m2 . (6.38) 

The diagrammatic representations are depicted in Figs. 6.1. a. 
The interaction part of the dynamic functional 

Jte[S, S] = Juc[S, S] + J RE [S, S] , (6.39) 

consists of a mode coupling vertex J7mc [S, S] , originating from the Larmor term S x ^ in 
Eq. Q6.12|) , and relaxation vertices Jke[S, S], originating from the nonlinear terms of the 



static Hamiltonian, L^. The mode coupling vertex is given by 
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Juc[S,S] = Xfe 



a(3-y 



k J ui J p J u 



xS a (k, u)S p (p -)S^(-p - -, -v - -) . 



(6.40) 



It consists of an isotropic part 



J M ci[S,S] = \fe a ^ [ [ [ /(p-k)^(k,a;)^(p_^_)^(p + ^ + ), (6.41) 

J k Jul Jp J v 

and a dipolar part 



JMcn[S,S]=\fge aPj [ [ [ [ ^tf±^(k,u,)^(p_, i/_)S"(p + , u + ) , (6.42) 

Jk Jui Jp Jv p+ 

whose graphical representations are shown in Figs.6.1.b-c. We have used the short hand 
notations p± = =F (p ± k/2) and v± = =F iy ± uj/2). 

The dipolar mode coupling vertex can be symmetrized by the substitution p — > — p in 
Eq. ( |6.41 ) leading to 



Jucd[S, S] — -Xfge a/ 3-y 

Z JkJuiJpJv 



^S a (k, w)^(p_, i/_)5"(p + , + 
pi 



(6.43) 



The relaxation vertex is given by 



Jbe[S,% = -±F°™ f If I I I 

o! J ki J K2 J k:j J u>i J u>2 •> u>c 



A 



7+(E k * 



x 



k-2 Jkz Jul 1 Ju!2 JU>3 

S a C£ ki, X)wi)^(-ki, -w 1 )5^(-k 2> -cu 2 )S 5 (-k 3 , -w 3 ) . 



(6.44) 



It contains a diffusive part 



u 



OI J ki J R2 J k$ J u>i J u)2 J U13 



X, 



:5 Q (E ki, E^)^("ki, -^)^(-k 2 , - W2 )5 5 (-k 3 , -w 3 ) 

i i 

and a relaxational part 



(6.45) 



J R 4S,S] = -^-XF a ^ 5 f f 

01 Jk\ Jk2 



7 



'&2 Jkz Ju>\ Ju)2 JUI3 

xS a (J2 k i; E^)^("ki, -ui)£T(-k 2 , -u 2 )S 5 (-k 3 , - W3 ) 



(6.46) 
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whose diagrammatic representations are given in Figs. 6.1.d-e. This decomposition is not 
important from a technical point of view, since both terms have the same tensorial structure. 
But, it emphasizes that the relaxational vertices due to the exchange interaction and dipolar 
interaction are diffusive and relaxational, respectively. 

With the above perturbation theory at hand, we can now calculate the Greens functions 
defined by 

„ai,...,a N ,f3i,...,f3r, , , 1 n 

"•jy jy w l) •• ) k N, ^N] KiV+1, ^JV+li k N+N> ^N+ff) ~ 

SjN+N 

d/i ai (ki,a;i) ■ ... ■ bh p N[k N+ jf,uj N+ i!) 
It is convenient to consider the vertex functions 

N ' l ' ^(kij^i; ...; kjv,a?jv; \s.n+i,^n+i', k N+N' ^n+n) ' (6.48) 



which can be obtained from the cummulants by a Legendre transformation [17 . 



B. Self consistent one loop theory 

In the preceding sections we have derived a path integral representation of the dynamics 
starting from semi-phenomenological Langevin equations. Within this method a perturba- 
tion theory for the correlation and response functions could be formulated which is similar to 
the usual field theoretic procedure for static critical phenomena. In each finite order of per- 
turbation theory infrared divergences arise reflecting the infrared divergences in the second 
derivatives of the free energy close to a critical point. In order to remove these divergences 
(which appear combined with the ultraviolet divergencies at the upper critical dimension 
of the model) and to have a perturbation theory with a small parameter one can use the 
concepts of renormalized field theories, i.e., upon introducing renormalization factors and 
expanding around the upper critical dimension. 

Instead we follow here a alternative route, which consists in the resummation of certain 
classes of diagrams to infinite order in perturbation theory. Such methods are frequently 
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used in condensed matter physics quite successfully. These methods are characterized by 
the fact that they lead to self consistent equations for the correlation functions. 

In this section we formulate such a self consistent procedure for critical dynamics. It will 
turn out that the resulting equations are equivalent to mode coupling theory. 



1. Self consistent determination of the line width in Lorentzian approximation 



As we have already seen in the discussion of the mode coupling theory one can obtain 
quite reasonable results for the line width by assuming that the line shape is given by a 
Lorentzian. In the above formulation in terms of a dynamical functional the line width is 
given by £L7| 

F?i(?,o)] 2 



1 lor — Z - 



(6.49) 



r?(g)r&(g,o) 

where T%(q) is the inverse static susceptibility (a = L,T). Upon using the Fluctuation 
Dissipation Theorem (FDT) 



rSL(g,0)= A( 7 + g 2 ) + A/r^ 10 (g,0) r»(5) 



(6.50) 



one obtains to one loop order — 2r" 1 (g,0) = I^ (q)TQ 2 (q, 0) and consequently r^ r = 
r?i(g,0) = -§rf(g)rg2(g,0). To zeroth order we have T% 2 {q,0) = -2A(7 + g 2 ). The 
one loop contributions to (g, 0) are shown in Fig. 6.2. The transverse and longitudinal 
Lorentzian line width is found to be 



Tf or (q)x T (q) = Kl + q 2 ) + ^ff 



p JU> 



+Vtl(9, P, q)C T (P-, uj)C L (p+, u)+ 
+vl L {g, p, q)C L (p_, u)C L {p + , u) 



(6.51) 



and 



Tt{q)X L {q) = Kl + q 2 ) + 2{\fY 



p Jul 



vtt(9, p, q)<^ T (p-, ^)c T (p+, w) + 

+Vtl(9, P, q)C T (P-, lj)C l {p + , u) 



(6.52) 
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where we have introduced the notation p± = p ± |. The vertex functions v°Jg, p,q) are 
given by 

1 



^tt(#>p> q) 



Vtl(9,P, q) 



d-1 



d-1 



cT - 4d + 6 



(p+-p-) 2 (p+q) 2 

2 9 



p%p_ 



ip-qj 



rf-3 + 



P+ • P 



+ 



ip+ ■ qj 



p\q 2 



yP ■ q + 



^IiO,p,q) 



^tt(#>p> q) 



i - 



d-3 + 2 



IP+ • P- 

p+p 2 . 



,p+ • qj 

p\q 2 



Vtl(9, p>q) = 4 



1 - 



(p+ • q) 

p\q 2 



ip-qj 



(p-q) 



(p-q + p 2 



(6.53a) 



(6.53b) 



(6.53c) 



(6.53d) 



(6.53e) 



Upon substituting in the integrals ( |6.51| ) and ( |6.52| ) p + = p + | — > p the above vertex 
functions are identical to the vertex functions in section III, Eqs. (|3.20|) - (|3.24|) . The frequency 
integration can readily be done leading to 



/ C a (p_,^( P+ ,^) 



x°(p-)x /3 (p 



(6.54) 



A [(7 + P 2 -)/x a (p~) + (7 + vl)lx^P + )] ' 
i.e., the product of the static susceptibilities divided by the sum of the line width to zeroth 
order. This structure corresponds to the first step in an iteration procedure and suggests 
to define a self-consistent approximation by replacing the bare line width by the full line 
width. 



x a (p-)x^(p 



X a {P-)x P (PA 



(6.55) 



A [(7 + p-)/x a (p-) + (7 + p+)/x /3 (p+)] [i? or ( P _) + rf or ( P+ 

This self consistent approach corresponds to a partial resummation of the perturbation 
series as shown in Fig. 6.3. The resulting equations are identically with the mode coupling 
equations in Lorentzian approximation found by the conventional derivation in section III. 
The factorization approximation (two mode approximation) in the conventional derivation 
corresponds to the structure of the one loop diagrams here. 
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2. Self consistent equation for the Kubo relaxation function 



In the preceding section we derived for the sake of simplicity the mode coupling equations 
in the so called Lorentzian approximation for the line shape. Here we go beyond this 
approximation and derive a self consistent theory for the full wave vector and frequency 
dependent relaxation functions. The dynamic susceptibility can be written as |l7j (a = L, T) 

\(l + q 2 ) + \fT«(q,cu) 



The X-insertion is defined by (17 



r?i(-g, -w) 



(6.56) 



X a (q, to) = e a/37 / / S p (k - q/2, v - uj/2)S\-k - q/2, -v - cu/2) . (6.57) 

J k J v 

Upon using the (explicit) structure of the dynamic propagators C a (q, t)Q(t) = X a (q)R a (q^ t) 
one can show that the relation 



r?i(-9, -o>) 



-ILO 



\(l + q 2 ) + XfTZ 10 (q,u) T a 2 {q) 



(6.58) 



holds to one loop order. This implies for the Kubo relaxation function Q a (q, to), which is 
related to the dynamic susceptibility by 

1 



^(q,u) 



■luJ 



bc a (q,u) - x(q)] 



(6.59) 



a structure analogous to Eq. (|3.25|) 



X a (q) 



ru-q,-u) [A(7 + g 2 ) + A/r^ 10 (g,^)]r^( g ) 

To one loop order one obtains for 10 (q, t) 



(6.60) 



T xw (q,t) = 2Xf [ Q(t) 

Jv 



Vtt(9> P, q)C T (p_, t)C T (p + , t) + 
+v^ L (g, P, q)C T (p_,t)C L (p + ,t) + 

+i$ L {g,P,q)C L (P-,t)C L (jP + ,t) , 



(6.61a) 



T XA0 (q,t) = 2Xf [ 0(t) 



Vtt(9,P, q)C T {p„,t)C T {p + ,t) + 

+vh(g,p,<i)c T (p-,t)c L (p + ,t) , 



(6.61b) 
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where the vertex functions are identical to Eqs. ( |6.53a| - |6~53e|) . We have C a (q, t)Q(t) = 
X a {q)R a {q,t). Eqs. ( |6.60[ - p3TrjD can be written in a self-consistent form by replacing the 
propagator C a (q,t) by the full Kubo relaxation function $ a (g,i) on the right hand side of 
Eqs. ( |6 . 6 1 a| J6T 6 1 b| ) . This corresponds, as in the preceding section, to a partial resummation 
of the perturbation series. The resulting equations are identical with the mode coupling 
equations Eqs. ( |3.19| )- (|3.25|) in section III . Hence we have shown that mode coupling the- 
ory is equivalent to a self consistent one loop theory. The factorization approximation in 
the conventional derivation of mode coupling theories is here a direct consequence of the 
structure of the one loop theory. Furthermore the present procedure has the advantage to 
be extendible to a self consistent formulation of higher loop order. In order to justify the 
validity of the mode coupling approach one has to ask in what sense higher order terms are 
small compared with the self consistent one loop theory. This question will be addressed 
in Appendix C, where we restrict ourselves to the isotropic case (in order to simplify the 
discussion) . 
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Captions to the figures: 



Figure 6.1: Basic elements of the dynamical perturbation theory, a) Correlation and 
response propagators, b) isotropic mode coupling vertex, c) dipolar mode coupling vertex, 
d) isotropic relaxational vertex, and e) dipolar relaxational vertex. 

Figure 6.2: One-loop diagrams for the vertex function (q, u). 

Figure 6.3: Partial summation of the perturbation theory for r^^q, uj), where vertex 
corrections are neglected. This resummation leads to the mode coupling equations. 
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Figure 6.1: Basic elements of the dynamical perturbation theory, a) Correlation and 
response propagators, b) isotropic mode coupling vertex, c) dipolar mode coupling vertex, 
d) isotropic relaxational vertex, and e) dipolar relaxational vertex. 
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Figure 6.3: Partial summation of the perturbation theory for r^f (q, co), where vertex 
corrections are neglected. This resummation leads to the mode coupling equations. 
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VII. CONCLUSIONS AND OUTLOOK 



The past decade has witnessed substantial progress in understanding the critical dynam- 
ics of real ferromagnets, both on the experimental and the theoretical side. In this review 
we have concentrated on the interplay between exchange and dipole-dipole interaction. We 
have presented theoretical and experimental evidence that many aspects of the critical dy- 
namics of ferromagnets such as Fe, Ni, EuO, EuS, and many other magnetic materials 
are now fairly well understood on the basis of a mode coupling theory which takes into ac- 
count exchange as well as the dipole-dipole interaction. In this final section we conclude our 
overview by summarizing some of the main theoretical and experimental achievements and 
by briefly discussing a number of open problems in the field of critical dynamics of magnets. 

One of the most important theoretical advances towards our understanding of the critical 
dynamics of real ferromagnets was the realization that additional interactions such as the 
dipole-dipole interaction can lead to a qualitatively new behavior of the frequency and wave 



vector dependence of the spin-spin correlation functions [fT9| , pO[ . A description of a mode 
coupling theory, which on top of the exchange interaction takes into account the dipole- 
dipole interaction, has been given in section III. The main results of this mode coupling 
analysis can be summarized as follows. The dipolar interaction leads to an anisotropy of 
the spin fluctuations with respect to the direction of the wave vector and introduces a 
second length scale q~ x besides the correlation length £. The presence of a second length 
scale leads to generalized dynamic scaling laws, containing two scaling variables, x = 1/qS, 
and y = q D /q for the length scale and one scaling variable r a = Aq z Q a (x,y)t for the 
time scale. Due to the dipolar anisotropy the characteristic time scales l/[Aq z Q a (x, y)} are 
different for the longitudinal and transverse modes (a = L, T). This is mainly due to the 
non critical longitudinal static susceptibility implying that the longitudinal characteristic 
frequency Aq z Q L (x,y) shows no critical slowing down asymptotically. Furthermore, the 
dipolar interaction leads to a quite interesting crossover in the line width and line shape of 
the spin-spin correlation functions. 
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Line width crossover. Precisely at the critical temperature the dynamic critical expo- 
nent for the transverse line width undergoes a crossover from z = 5/2 to z = 2, which is 
displaced with respect to the static crossover to wave vectors smaller by almost one order 
of magnitude. This explains why up to now this crossover escaped detection by neutron 
scattering experiments right at T c . For the longitudinal line width the crossover is predicted 
to be from z = 5/2 to an uncritical behavior z = and, in contrast to the transverse width, 
it occurs in the immediate vicinity of the static crossover, characterized by q D . If the line 
shapes are approximated by Lorentzians the concomitant Lorentzian line widths obey the 
scaling law T a (q, £, g) = Aq Zr y a (x, y). For vanishing dipolar coupling g, the scaling functions 
coincide with the Resibois-Piette scaling function. If the strength of the dipolar interaction 
g is finite, the curves approach the Resibois-Piette scaling function for small values of the 
scaling variable x and deviate therefrom with increasing x. Since the dipolar interaction 
leads to a non conserved order parameter, the line width in the hydrodynamic limit is given 
by T oc q° instead of T oc q 2 as for an isotropic exchange ferromagnet. 

Line shape crossover. Close to T c the line shapes of the longitudinal and transverse 
relaxation function coincide in the isotropic Heisenberg limit, i.e., for large values of the wave 
vector q (q ^> q D ). In this limit the dipolar interaction becomes negligible and the shape is of 
the Hubbard- Wegner type as discussed in section II. Upon increasing the value of the scaling 
variable r = ^ (l/g£) 2 + {Qd/iYi ^ ne ^ me shapes of the transverse and longitudinal relaxation 
function become drastically different. Whereas the transverse relaxation function shows a 
nearly exponential temporal decay, over-damped oscillations show up for the longitudinal 
relaxation function. The line shape crossover sets in in the vicinity of the dipolar wave 
vector q D in contrast to the line width crossover, which starts at a wave vector almost one 
order of magnitude smaller. The situation is quite different for temperatures well separated 
from the critical temperature. Then, the difference in the shape crossover of the longitudinal 
and transverse relaxation function diminishes with decreasing q D ^. For q D ^ <C 1 the shape 
crossover as a function of r corresponds to the crossover from the critical (Hubbard- Wegner) 
shape to the hydrodynamic shape as discussed in section II. 
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On the experimental side, an important development began with the observation that the 



data for the line widths above the transition temperature | |188| , |189|J3T[1 in magnetic materials 
such as EuOsaid Fe could not be described by the Resibois-Piette scaling function resulting 
from a mode coupling theory ||231|| as well as a renormalization group theory ||1 74| . |1 23 



which take into account the short range exchange interaction only. Even more, the data 
could not be collapsed onto any other single scaling function. Those experimental results 
were explained quantitatively on the basis of the mode coupling theory |79 HSl|l described 
in section III. Clear indications of the importance of the dipolar interaction have also been 
observed in hyperfine interaction experiments on Fe and Ni, where one found a crossover in 



the dynamical critical exponent from ,2 = 5/2 to £ = 2 |[228|J^|iT3| , |114|| , and with electron 
spin resonance and magnetic relaxation experiments ||147H149| [|15(|[)4],|152| , where a non 
vanishing Onsager coefficient at zero wave vector was found. For all these data now there 
exists a quantitative theoretical description (see section IV). 

With advances in the neutron scattering technique subsequent experiments investigated 
the line shape of the spin-spin correlation function. This progress made it possible to test 
the predictions for the line shape from MC-theory | |262| , |118| | and RG-theory |5"T| . |2"2] ] , which 
take into account solely the short-range exchange interaction. Originally those theories for 
the critical dynamics were expected to fit the experimental data in an increasing quan- 
titative way as one was able to measure the correlation functions at smaller and smaller 
wave vectors, since the influence of additional irrelevant interactions should diminish as 
one moves closer to criticality. It came as quite a surprise when Mezei found a nearly 
exponential decay of the Kubo relaxation function by spin echo experiments on EuO at 
q = 0.024A" 1 and T = T c ||190||, which was in drastic disagreement with the bell-like 



shape predicted |262| . |118| . pT| . [22]| . The anomalous exponential decay found by Mezei ||190 



was shown again to be a dipolar effect |8^J8^]. Recently, significant experiments have been 
performed using polarized neutrons which allowed for the first time to measure the 

longitudinal and transverse spin-spin correlation functions separately. The results confirmed 



the theoretical predictions concerning the line width |80(] as well as the line shape of the Ion 
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gitudinal relaxation function f82| , |85|| . The experimental data, de-convoluted by a maximum 
entropy method, agree quite well with the theoretical predictions, especially the double peak 
structure corresponding to the over-damped oscillations in time have been observed. This 
can be regarded as a success of the mode coupling theory 

For the critical dynamics below T c experimental investigations are scare. In contrast to 
the very satisfactory situation above T c , where quantitative agreement between experiment 
and theory has been achieved, the situation in the ferromagnetic phase is by far less clear. 
Theoretical investigations so far have concentrated on the dynamics of isotropic exchange 
ferromagnets, neglecting dipolar effects |8lJ240|| . In this context it has been shown within the 



framework of mode coupling theory that the amplitude of the scaling function for the spin- 
wave frequency of the isotropic Heisenberg Hamiltonian is universal ||240|| . A comparison 



of the theory f81fl , taking into account the so determined universal amplitude ||240|| , with 



recent measurements of the longitudinal line width on Ni |34| below and not too close to 
the critical point gives quantitative agreement. 

But also in the ferromagnetic phase profound effects of the dipolar interaction on the 
spin dynamics have been revealed in recent experiments. Using polarized neutron scattering 
the spin excitations in EuS and Pd 2 MnSn have been investigated below T c by Boni et 



al. |37|j38[| . These investigations indicate that there are deviations from the theoretical 
correlation functions of isotropic exchange ferromagnets. Especially, an anisotropy of the 
spin fluctuations with respect to the polarization relative to the wave vector q has been 
observed. Still further experimental and theoretical progress is needed in order to reach 
a thorough understanding of the spin dynamics in the ferromagnetic phase. The most 
convincing experimental indication for the importance of the dipolar interaction below T c has 
been obtained recently in a measurement of the homogeneous magnetization dynamics []59 



It is found that the scaling function for the kinetic coefficient below T c coincides with the 
one observed earlier above T c . Up to now, there is no theoretical explanation for these 
interesting experimental findings, but work in the framework of mode coupling theory is in 
progress [|241| . 
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Besides three-dimensional ferromagnetic materials with isotropic exchange interaction 
there are other magnetic systems, where the dipole-dipole interaction influences the criti- 
cal behavior, such as antiferromagnets or ferromagnets with a uniaxial or planar exchange 
anisotropy. The critical dynamics of isotropic dipolar antiferromagnets is affected by the 
dipolar interaction featuring a crossover from diffusive to relaxational dynamics of the mag- 
netization and a crossover in the dynamic exponent ||. The effect on the staggered mag- 
netization is much less pronounced than for the corresponding order parameter modes in 
the ferromagnetic case. For ferromagnets with a uniaxial exchange anisotropy the dipolar 
interaction leads first of all to a reduction of the upper critical dimension from d c — 4 to 
d c = 3. As a consequence the critical statics is described by mean field theory with loga- 
rithmic corrections [|161|| . A time dependent Ginzburg Landau model |76j for the dynamics 



results in logarithmic corrections to the conventional van Hove theory of critical dynamics. 

We conclude with some comments on two-dimensional (2D) magnetic systems. These 
are interesting because the Hohenberg-Mermin- Wagner theorem ||lliyi86|| excludes conven- 



tional long-range order in isotropic systems with a continuous symmetry and short-ranged 
interaction. Two dimensions seem to be the borderline dimension where thermal fluctua- 
tions are just strong enough to prevent the appearance of a finite order parameter. However, 
the anisotropy of the dipolar interaction with respect to the separation of two spins leads to 
a suppression of the longitudinal fluctuations and thus a finite order parameter may exist 
in two dimensions (see section V). The dipolar interaction is responsible for the existence 
of long range order in both ferromagnetic ||177|| and antiferromagnetic thin films ||215| , |216 



There is quite an interesting phase diagram with various spin flop and intermediate phases 
for two-dimensional dipolar antiferromagnets in a finite external magnetic field ||217| . 



Two-dimensional ferromagnets are interesting for practical reasons as magnetic storage 
devices as well as from a fundamental point of view. With the recent advances in thin 
film technology it becomes possible to study the critical behavior of two-dimensional mag- 
netic systems [nj. Due to the combined effect of anisotropy, dipole-dipole interaction, and 



increased importance of thermal fluctuations in two dimensions, the phase diagram of fer- 
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romagnetic thin films shows a variety of new phases [ 
in this field are still rather controversial p!2| , pi3|J164 . 
question is the critical dynamics of all those phases. 
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[131 



lTIITTf . The theoretical predictions 
. A still open and very interesting 
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APPENDIX A: FLUCTUATION-DISSIPATION RELATIONS 



In this appendix we collect several fluctuation dissipation theorems, which are of impor- 
tance for the dynamics of systems described by nonlinear Langevin equations. 

Let T be the time reversal operation: t — ► —t. Then detailed balance (time reversal 
symmetry) implies ||128|| that 



T exp \Jll - H tl \ = exp \JZll - H^ t2 \ , (Al) 

where H t = TC(S(t)) is the stationary probability distribution function and J^[S,S] the 
dynamic functional. For simplicity we have assumed a one-component field S. Without loss 
of generality one can assume that the field S(x,t) is even or odd under time reversal, i.e. 

TS = eS, e = ±1 . (A2) 

In the case of spins S is odd under time reversal. The stationary distribution P^S] = e~ n ^ 
is characterized by the "free energy" TifS 1 ]. The time reversal symmetry implies that 

^>H*H>-«1). 

Now one uses the causality property of the response functions 

< S(t 1 )S(t 2 )...S(t k )S(t 1 )S(t 2 )--S{t k ) >= 0, if one tj > all U . (A4) 
Then for example 

< S(t)S(0) >= for t < . (A5) 
With the time reversal operation and Eq. (|A3| ) it then follows from Eq. ( |A5|) for t < 

Upon redefining t = —t one obtains for t > 

< s(t)s(o) >= e(t) < s(t) s -^M > (A?) 
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The same arguments can be repeated for < S(ti)S(t2)---S(tk)S(ti) > with ti > allt,-. The 
result is 



< s(ti)S(t 2 )...s(t k )s(i 1 ) >= e(ti,{tj}) < s(t l )s(t 2 )...s(t k ) s ^.fl )] > 



(A8) 



where 0(ti, {tj}) is an obvious generalization of the O-function. Note that these generalized 
FDT's hold for the cummulants and not for the vertex functions. In particular we get 



vk 2 

G u (k,t) = e(t)—G 02 (k,t). 



(A9) 



In a completely analogous way one can derive the following identities 1 17,128 



Xa p{x -x',t- = -Q(t - t')j t < S a (x, t)S%x', t') > 



(AlOa) 



UJ 



(AlOb) 



- e(t-t')— < S a (x,t)S^(x',t') >= L lP < S a (x,t)S\x',t') > -\f < S a (x,t)XP(x',t') > 

(AlOc) 

with X 13 = €fs^ u S fJ, S u . Therefrom one can deduce the following identities for the vertex 
functions 



rff (q, w = 0) = \q a 5 a ^ + \fT% 10 (q, u, = 0) Vf (q) , 



(Alia) 



- (q, w) = M a Tn p (g, u) - T"?(-q, -u 



a/3, 



7/3 



(Allb) 



where Tf{q) =< S a {-q)S f3 {q) > 1 is the static two-point vertex function, and T 



X,NN 



denotes a vertex function with one X-insertion. If the Hamiltonian is quadratic in the fields 
{S}, i.e., Ti = \ J d d xc a pS a S 13 then one has a further identity 



< S a {t)S^(t') >= Q{t - t')c aP < S a {t)S /3 (t') > . 



(A12) 
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APPENDIX B: DERIVATION OF NON LINEAR LANGEVIN EQUATIONS 



In this Appendix we give a short derivation of non-linear Langevin equations based on 
the projector formalism of Mori [197 and Zwanzig [275]. In our presentation we follow 



closely the papers by Kawasaki [|139|1 and Mori et al. ||198|,199 l. 

As a first step in deriving non-linear Langevin equations for a particular physical sys- 
tem, one has to choose the relevant slow variables. This choice is in general dictated by 
the conservation laws and spontaneously broken symmetries in the system under consid- 
eration. Let us assume that we have found an appropriate set of slowly varying variables 
{S} = (Si, S2, • • • , S n ). Then, the next step consists in separating an arbitrary dynamical 
variable X into a part that is associated with the slow variables {S} and the rest. This is 
accomplished by introducing an projection operator by 

VX(t) = ^2(X,* j )<S> j ({S}), (Bl) 

j 

where one defines the inner product as the value of the Kubo relaxation function at t — 

(A, B) = $ AB (t = 0) = 1 lim / dre- eT < [A, B^} > \ t=0 , (B2) 

which is identical with the static susceptibility. The functions <3> 3 -({S'}) and ^/jdS"}) are two 
suitable sets of functions which are orthonormal with respect to the inner product defined 
in Eq. flB2l) 



(B3) 



The microscopic dynamics of an arbitrary dynamic variable X following from the Heisenberg 
equations of motion can formally be presented by the Liouville equation 

d 



dt 



X(t) = iCX(t) . 



(B4) 



It was shown in Ref. ||139|| that with the aid of the operator identity, 

-e ltc = e itc iC + [ l dse i{t - s)c i£oe it{c - Co) i(£ - C ) 
/ Jo 



d 

dt' 



+e 



i(C - C c 



(B5) 
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one can derive a formally exact non-linear Langevin equation 

j t x(t) = ^(zCx^,)M{s(t)}) 

-t 



"E / dr (f x (r), f p ) § p {{S{t - r}) + f x (t) , (B6) 
p J0 

with the "random" forces 

fx(t) = e? tQC iQCX , (B7a) 

fa=Qi£*a({S(t)}), (B7b) 

where Q = 1 — V. Note that the adjoint operators of C and V are denoted by £ and V. 

The first term in Eq. ( |B6| ) is the so called "adiabatic" term, and it describes the time 
variation of the dynamic variable X which follows adiabatically the changes in the slow 
variables {S(t)}. The second term represents the damping of this adiabatic motion. The 
last term is the stochastic force acting upon the dynamical variable X. Note, however, that 
this interpretation depends on the appropriate choice of the slow variables and the sets of 
functions $ a ({5(t)}) and ty a ({S(t)}). Upon choosing (we assume that the statics is diagonal 
in the variables S a ) 

®a({S(t)}) = X - 1/2 S a , (B8a) 
y a ({S(t)}) = X - 1/2 S a , (B8b) 



with (S a ,Sp) = S a/ 3Xa one obtains Mori's generalized Langevin equations [[L97|| 
d 

where 



j t S a (t) = "£iuj aP S p (t) -Y,J*dr (f a (r), f fi ) S p {t - r) + f a (t) , (B9) 



iu a p = (i£S a , S^Xa 1 , (BlOa) 
f a (t)=e itQC Qi£S a . (BlOb) 

The latter equations are the basis for the mode coupling theory discussed in sections II 
and III. However, with the above choice of a linear projection operator the "random" 
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forces are orthogonal with respect to the slow variables S a only. It has been shown in 
Refs. [ |138j , |140| , p77|| that the "random" forces may not be really random since they contain 



products of the slow variables S a . Another choice for the projection operator is to include 
all the suitable symmetrized polynomials of {S} among the sets of functions {^/} and {$}. 
Now we restrict ourselves to the case of classical mechanics, where one has = {$}• 
This choice corresponds to the projection operator introduced by Zwanzig [|275| , [276|1 . The 
completeness relation for the Zwanzig projection operator reads 

E <M{s})^(m) = (an) 

with the equilibrium probability distribution function 

1 



P eq ({S}) = exp 



H({S}) 



(B12) 



k B T 

The resulting generalized Langevin equation takes the form 

jjS„(t) = »„({$(«)}) 

+£ [ * W(t - T)] f^W(t-mPMS(t-r)}) + ut) (B13) 

with the so called mode coupling term 

v a ({S(t)})=<i£S a ;{a}>, (B14a) 

and the memory kernel 

M a(3 (t- {a}) =< f a (t)t p (0); {a} > , (B14b) 

where we have introduced the notation < X; {a} >=< X5(S — a) > /P eq ({a}) for the 
conditional average. 

The semi-phenomenological equations of motion are obtained from the above exact equa- 
tions by making three basic, plausible assumptions, (i) First one makes a Markovian ap- 
proximation for the kinetic coefficients. This is justified by the fact that one has included 
all suitable symmetrized polynomials of the slow variables in the projection operator. 
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M af} (t; {a}) « 2L a(3 ({a})5(t) . 



(B15a) 



(ii) Second one assumes that all kinetic coefficients are independent of the slow variables 



L a p({a}) m L al3 . 



(B15b) 



(hi) Third one takes the random forces ( as Gaussian white noise 

1 /•*! 



w({(}\to < t < h) ~ exp 



dtCa(t)L-fa(t) 

4 Jto 



Then the generalized Langevin equations reduce to 



where 



-^(t)=MW)})-£L 



<w({s(o» 



+ Ca(t) , 



(B15c) 



(B16) 



H({S}) = -k B Tln(P eq ({S})) , 



(B17a) 



^«(W)}) = -a/E 



The Poisson brackets are defined by 



-QMS}) - Qapds})^^- 



(B17b) 



<2^({s}) = s p }p = -QMS}) ■ 



(B18) 



We close this appendix by presenting the generalized Langevin equation for the example of 
an isotropic ferromagnet 

H = \j d d x[rS 2 (x,t) + (VS(x,t)) 2 ] + »f d d x(S 2 (x,t)) 2 . (B19) 

The Poisson brackets are given by 



Q a/3 (k,k') = e a/37 ^(k + k'). 



(B20) 



Since the order parameter is conserved the kinetic coefficient is 



L(k) = Xk 2 , 
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(B21) 



(in general L(k) = \k a ). Hence one obtains the following generalized Langevin equation for 
isotropic ferromagnets 

Jt s ^ " = A/s x mhj ~ A < iV > to + « x ' ( » • < B22 > 

Note that the conventional (van Hove) theory of critical dynamics is obtained by making the 
following additional assumptions: (i) The Onsager coefficient L remains finite at the critical 
point, (ii) The mode coupling term is ignored: / = 0. (iii) A quadratic approximation is 
made for the free energy functional TC = \ J k x _1 (k)S(k, £)S(— k, t), where x(k) is the static 
susceptibility. 

One should also note that in mode coupling theory solely the latter approximation is 
made (and further additional approximations like two-mode approximation and neglecting 
the vertex corrections). 

APPENDIX C: VALIDITY OF MODE COUPLING THEORY, HIGHER ORDERS 

IN PERTURBATION THEORY 

In this section we analyze how higher orders in perturbation theory can be incorporated 
in the self consistent approach. For the sake of simplicity and clarity we restrict ourselves to 
the case of isotropic ferromagnets with exchange interaction only. The extension to dipolar 
ferromagnets is straightforward. 

First we consider the corrections from two-loop contributions to the vertex functions r 02 
and r n shown in Figs. A.l and A. 2. Due to the tensorial structure of the relaxation and 
mode coupling vertices they can not be contracted (F a ^ 5 e a 'js = 0). 

The 2-loop diagrams of r n and r 02 can be divided into two categories. The first category 
consists of "true" 2-loop diagrams connecting two relaxation vertices. The second type 
consists of diagrams whose structure is - up to vertex corrections - identical to the one-loop 
diagrams. Upon defining a renormalized mode coupling vertex, as shown in Fig. A. 3, the 2- 
loop contributions to the vertex functions Tn and To2 are given by the diagrams in Figs. A. 4 
and A. 5. 

186 



The modifications of the self consistent equations resulting from these types of two-loop 
diagrams are now discussed separately. 

(I) The vertex corrections in Fig. A. 3 correspond to the vertex function Ti2(q,p, v = 0,u). 
Since all diagrams have a bare mode coupling vertex with an external response line (see e.g. 

the one-loop contribution to T 12 (q,p, v — 0,co) is proportional to p • q, i.e., it has the 
same wave vector dependence as the original bare mode coupling vertex. The remaining 
frequency integral is a function f(uj) of the external frequency u. Because only the long 
time behavior is of importance for the critical dynamics the limit u — > merely leads to a 
renormalization of the amplitude, Xf, of the mode coupling vertex and not to a change in its 
wave vector dependence. (But, see the discussion of higher order vertex corrections later) 

(II) The two-loop diagram connecting two relaxation vertices is proportional to 

TqT ~ AY / / / C(p, u)C(k, c)C(q - p - k, u) . (CI) 

Jp J k Juj 

Upon introducing the dimensionless wave vector variables p = p/q and k = k/g, carrying 
out the frequency integrals, and replacing the correlation propagators according to Eqs. 
( p.55| ) one finds 

It* ~ q±~* . (C2) 

In comparing this result with the wave vector dependence, q 3 ~ z , of the one-loop diagram, 
connecting two mode coupling vertices, one recognizes that the contribution of Tq 2 r can be 
neglected in the long wave length limit. Hence this two- loop diagram does not change the 
self-consistent equations as well. 

The above discussion shows that in a self consistent theory the two-loop diagrams are 
small compared to the one- loop diagrams in the long time and long wave length limit. 

Now we discuss higher order vertex corrections, where an internal line of a one-loop vertex 
correction is decorated with a static insertion. It can be shown that the renormalization 
of the mode coupling vertex is of purely static origin ||87|| . This is due to an exact relation 
between the renormalization factor of the mode coupling vertex and the time scale and field 
renormalization. To lowest order the mode coupling vertex is given by 
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ri 2 (p, q) = f 3 (p, q) = A /oea/3 7 P • q. 



(C3) 



This expression can also be written in terms of the bare static susceptibilities as 



^s(p,q) 



(C4) 



showing that the wave vector dependence of the mode coupling vertex is determined by the 
static susceptibilities. In order to take into account the static IR divergences correctly, one 
way to proceed is to replace the bare static susceptibilities by the fully renormalized static 
susceptibilities 



in the mode coupling vertex. This replacement is suggested by the exact relation between 
the renormalization factors. It is also equivalent to replacing the nonlinear Landau- Ginzburg 
functional by 



where x _1 (k) is the fully renormalized static susceptibility. This approximation is often 
used in mode coupling theories without any justification. Here it is a consequence of the 
exact relation between the renormalization factors. 

Besides these static vertex renormalizations (due to the non linearity in the Landau- 
Ginzburg functional) there could also be dynamic vertex renormalizations. But this is not 
the case here as we will show next. This can be seen from the diagrammatic representation 
of the dynamic vertex renormalization (see Fig. A. 3. Any diagram contributing to the vertex 
corrections starts with a bare mode coupling vertex. Therefore to leading order in the wave 
number one can replace the remaining wave vectors (in the integrals of the corresponding 
diagrams) by their values at p = q = 0. This leads then only to changes of the amplitude. 
Hence we have 




(C5) 




(C6) 




(C7) 
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This gives a precise specification of the approximations involved by neglecting the vertex 
corrections. 

As we have seen above, the effects of the four point coupling can be taken into account 
by a static renormalization of the mode coupling vertex, Eq. ( |C7[ ). Then one is left with 
a dynamic theory with a harmonic (effective) Landau- Ginzburg functional and a mode 
coupling vertex. In this case the fluctuation-dissipation relation ( |A12| ) applies (see Appendix 
A), which in the case of a purely quadratic Landau- Ginzburg functional reduces to 

G n (k,t) = e(t) X -\k)G 02 (k,t) (C8) 

This implies for the Fourier transform Gu(w) + Gn(-w) = x^Gmiyi) and consequently 

2Re[G u (u;)]=x- 1 Go2(uj) (C9) 

By using the generalized FDT's in Appendix A we find 

iu -iu + [Xq 2 + A/r x ,io(q, w)]/x(q) 

This is exactly the same structure as the non linear Langevin equations (with the linear 
projection operator) (see section III) with the memory kernel 

x(q) 

Neglecting the mode coupling contribution (/ = 0) leads immediately to the conventional 
theory (van Hove theory) of critical dynamics. The mode coupling contribution to the 
memory function is given by 

1 



r Xil0 (q,i) = A/ / 



X (q + k/2) x(q-k/2) 



-$(q + k/2,t)$(q-k/2,t). 



x(q + k/2) 

(C12) 



In the strong coupling limit, where the van Hove term Xq 2 (in general there will be a 
crossover from van Hove to strong coupling) can be neglected, a scaling analysis of the 
above self consistent equations shows that the dynamic exponent is given by 
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z = i(d + 2-r ] ) (C13) 

i.e., with this generalized mode coupling theory we find the correct dynamic exponent. This 
resolves a long standing problem with the conventional derivation of the mode coupling 
theory, which gives a dynamic exponent with the wrong sign of rj. The correct expression, 
Eq. flC13|) , is a consequence of the fact that we have correctly taken into account the vertex 



renormalizations of the mode coupling vertex by static decorations of internal lines. 

The above analysis is not a rigorous derivation of a self consistent theory. It shows, 
however, that the neglected terms are small in the long time and long wave length limit. 
Furthermore, taking into account static corrections of the mode coupling vertex - suggested 
by an exact relation between the renormalization factors of the mode coupling vertex and 
the static field renormalization - we could show that the self consistent theory (generalized 
mode coupling theory) is capable of giving the exact result for the dynamic critical exponent 
with the correct sign of rj. 
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Captions to the figures: 

Figure A.l: One-loop and two-loop diagrams contributing to the vertex function r n . 

Figure A. 2: One-loop and two-loop diagrams contributing to the vertex function r 02 
(the self energy). 

Figure A. 3: Diagrams contributing to the one-loop vertex correction. 

Figure A. 4: One-loop and two-loop diagrams contributing to the vertex function Tn, 
drawn in terms of the renormalized vertex. 

Figure A. 5: One-loop and two-loop diagrams contributing to the vertex function r 2, 
drawn in terms of the renormalized vertex. 
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Figure A. 2: One-loop and two-loop diagrams contributing to the vertex function r 02 
(the self energy). 
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Figure A. 3: Diagrams contributing to the one-loop vertex correction. 
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Figure A. 4: One-loop and two-loop diagrams contributing to the vertex function r n , 
drawn in terms of the renormalized vertex. 
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Figure A. 5: One-loop and two-loop diagrams contributing to the vertex function T, 
drawn in terms of the renormalized vertex. 
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